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SUMMARY 


This report is the first supplement to the primary HILT OP program 
document published in January, 1973, and describes the modifications and im- 
provements made to the HILTOP electric propulsion trajectory optimization 
computer program up through the end of 1974 * 

New program features include the simulation of power degradation, 
housekeeping power, launch asymptote declination optimization, and powered 
and unpowered ballistic multiple swingby missions with an optional deep space 
burn,, 

The report contains the new analysis describing these features, a com- 
plete description of program input quantities, and sample cases of computer 
output illustrating the new program capabilities . 
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I. 


INTRODUCTION 


This document describes the modifications and improvements made to 
the HILTOP electric propulsion trajectory optimization computer program [l] 
during calendar years 1973 and 1974, New program features include the simu- 
lation of power degradation, housekeeping power, launch asymptote declination 
optimization, and powered and unpowered ballistic multiple swingby missions 
with an optional deep space burn. 

The power degradation model has been hypothesized by the authors in 
earlier publications [2, 3l . The model allows a single parameter (denoted 
"characteristic degradation time") to describe the power degradation behavior 
of an electric propulsion spacecraft to a degree which fundamentally affects 
the solution to the trajectory optimization problem. 

The option of simulating spacecraft housekeeping power applies to solar 
electric propulsion with specified reference power. The housekeeping power 
is a specified constant power generated by the solar arrays and shunted away 
from the thruster power-conditioners and directly to the spacecraft payload 
for ,r housekeeping" purposes. 

The launch asymptote declination optimization model was first hypothesized 
by the authors in the appendix of [l] , and later a more thorough treatment of the 
subject was put forth in [4] . A solution to the problem of optimizing electric 
propulsion heliocentric trajectories, including the effects of geocentric launch 
asymptote declination on launch vehicle performance capability, is developed 
using variational calculus techniques. The model of the launch vehicle per- 
formance includes a penalty associated with a non-easterly launch plus another 
penalty arising from a non-coplanar launch from the parking orbit. Provisions 
for range safety constraints are included. Optimal trajectories wall generally 
have the launch excess velocity offset from the initial primer vector. 
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Extra-ecliptic mission simulations now involve launches from the Earth 
in which the Earth’s ephemeris is generated by the program’s analytic ephemeris 
capability; previously, extra-ecliptic missions were generated simply by starting 
the trajectory on the x-axis at one astronomical unit from the sum This improved 
extra-ecliptic capability allows the launch date to be optimized together with the 
launch asymptote declination. Also, an additional set of boundary conditions has 
been added to the program for simulating extra-ecliptic missions. 

The program's ability to simulate all-ballistic missions has been ex- 
panded to include powered and unpowered multiple swingby missions with an 
optional deep space burn. This capability renders HILTOP a powerful tool for 
ballistic mission design and optimization, with tremendous flexibility for creating 
imaginative multi-target mission profiles. 

This report contains the analysis describing these new program features, 
a complete description of program input quantities, and sample cases of computer 
output illustrating the new program, capabilities. 
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II. FORMULATION 

A. SPACECRAFT AND TRAJECTORY MODELS 

The assumed spacecraft and trajectory models are as described in Ll] 
and are not repeated here; the nomenclature used in the following analysis is 
also described in [l] , except for the introduction of new symbols which are 
described in the text as they appear, 

l a Power Degradation, Historically, this electric propulsion power 
degradation model first appeared in the literature in [2] , and then soon after- 
ward an improved discussion appeared in [3] , in which several of the ramifi- 
cations and consequences of the theory are also discussed. For deeper insight 
into the analysis, the reader is therefore referred to [3] , from which much of 
the analysis below is extracted. 

The model discussed here is concerned with the manner in which the 
performance of a solar array degrades due to high energy particle damage. 

It is assumed that one can define a damage factor q which has a value of unity 
at the initial time and which decreases in value with time during the course of 
the mission such that the power output, p, of the arrays at any time may be 
written 

p = q yp ref 

where p „ is the reference power and y is the power factor which is a 
ref 

function of solar distance and array orientation relative to the sun. The damage 
factor q may also be thought of as a time -dependent efficiency. The derivative 
q is negative, and is assumed to be linearly proportional to q and to the density 
of particles impinging on the solar cells. For simplicity, it is assumed that the 
particle density d is of the same form as the density of photons striking the 
surface of the array, i.e., 

e 0 n 

a = -4- (i) 

r 
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where e is a unit vector along the sun-spacecraft line, n is a unit vector 
r 

normal to the arrays such that e • n > 0 implies the front of the panels faces 
the sun, and r is the solar distance of the spacecraft,, In the discussion to 
follow, an upper case symbol will denote a vector, a lower case symbol will 
denote a scalar, and a lower case symbol with a bar will denote a unit vector. 
Then 



q = “ kqd 

(2) 

and 

t 



-k f ddt* 



q = e J 0 

(3) 


where k is the constant of proportionality. It is convenient to introduce a 
parameter s, called degradation time, 

t 

s = J ddt' (4) 

0 

or 

s =d (5) 

Thus, under the above assumptions the degradation time accumulates at a faster 
rate when the spacecraft is nearer the sun, which is a characteristic one might 
expect. Note that when the spacecraft is maintained at 1 AU with the panels 
normal to the sun line, s = 1 and the degradation time is equal to the flight 
time. Then k may conveniently be thought of as the inverse of a reference time, 
called the characteristic degradation time denoted by T , and the degradation 
factor q becomes 

-s/t, 

q 53 e d . (6) 

Actually, there is little reason to allow the degradation time s to continue to 
accumulate during coast phases since the arrays may be turned edgewise to the 
sun and the degradation process may be effectively halted. Therefore, we ad- 
just equation (5) to read 

s = h d (7) 

0 v 
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where h is a step function equal to one during thrust phases and equal to zero 

O’ 

during coast phases,, 

The characteristic degradation time T is an engineering parameter that 

must be determined experimentally. For example, by exposing a solar cell to 

the particle emission of a solar simulator and measuring the performance of the 

cell over a period of time, one should be able to estimate a reasonable value of 

T , Another source of information would be measurements from actual space - 
d 

craft which employ solar cells for power supply. 

The assumed exponential form of the degradation factor, although intended 

for use with SEP systems, is applicable for NEP systems as well. The principal 

difference is in the definition of s ; for example, s* = h^,. The exponential form 

also permits one to evaluate radio-isotope systems by defining s = 1 and letting 

T represent the time for the radioactivity to dissipate to l/e of its initial 
d 

level . 

In the development which follows, the formulation applicable to SEP is 
used exclusively. The equations of motion and adjoint equations are given in 
[l] ; the equations affected by the inclusion of power degradation in the model 
are given below. 

Power degradation affects the problem in a very fundamental sense, be- 
ginning with the rocket-thrust term in the equations of motion: 

V - R = a e t - ^“R 

r (8) 

R = V 

where R is the position vector, r is the magnitude of R, V is the velocity 
vector, fi is the gravitational constant of the sun, a is the magnitude of the 
thrust acceleration and e is a unit vector in the direction of thrust. The 
thrust acceleration a is a function of several variables and may be written as 


follows : 



a ~ h 

a 


gyq 

V 


where g is a reference thrust acceleration evaluated under a prescribed set of 
conditions, V is the ratio of current to initial mass, and q is the degradation 
factor defined above . The power factor y is assumed to be a function of the 
density, d, of photons incident on the arrays, where d is as written in equation 
(1)* The general form of y is 
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i=0 



( 9 ) 


The coefficients a. are chosen so that this equation will adequately describe the 
power output of a given array. The only restriction on the is that their sum 
is equal to one. Then at r = 1 AU, with the arrays normal to the sun line, 
y — d = 1. The summation term in (9) represents the temperature effect. 

The mass ratio satisfies the differential equation 

• , K 7q 

v = - h , (10) 

<y c 


using v - 1 as an initial condition, where c is the jet exhaust speed which is 
assumed to be constant over the trajectory. 

The variational Hamiltonian becomes 

h = A * R - A • R + X v + \ r + \ s+Xg+Ac 
v v r s g 5 c 

= h aR^< A - V^V + X s d+X r]- ^< A ” R >- A ‘R W 

r 


and the adjoint equations are 


A=--^ A + ^(R. A)R+h o [^(A-S t -^ X ) + A ] 


ad 

Sr 




( 12 ) 
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( 12 ) 

(eont. 


X- =-h -23- (A *5 - * A> 

g (7 P t C V 


X = -h ^A 
c a 2 v 
c 


where 


gyq 




A' = h (A - e - - A ) 

s cr i/t, t c v 

d 


y* = 3y/3d = l a.(l + J-)d l/4 
i=0 


dd _ 1 

dR ~ 3 
r 


n - 3(e • n) e J 


The control variables are the thrust direction e^_ , the switch step 

function h , and, providing the array orientation is not constrained to yield 

a 

maximum power, the normal direction n. According to the Maximum Principle, 
these controls are chosen to maximize the variational Hamiltonian (11). The 
maximum of h^ with respect to h^ is seen to depend totally on the sign of the 
term in square brackets. That is, denoting 


a = ^ (A 


e - — A ) + A d + A , 
t c v } s T 


then choose 


f 1 if a > 0 
to if a < 0 


(13) 


(14) 


Maximizing h with respect to e, is also accomplished by inspection. Since 
v t 

appears only in the dot product with A and since the coefficient of that dot 
product is non-negative, i„ e 0 , 


h y gyq/l' a 0 


then h is maximized with respect to e. by making A" e, as large as possible, 
v t t 


7 



which is the same result as when there is no degradation., The control vector 

n appears explicitly in (11) through the density d as given by equation (1), In 

fact, since n only appears in h through d, n only affects h^ through its dot 

product with e , Letting the angle between n and e be denoted X such that 
1 r r 

e 8 n - cos X > then it is clear that thei’e will be a ,f best n angle X between n 
r 

and e to maximize h , but that n may lie along any element of a right circu- 
r v 

lar cone of half angle X about e . For the moment, we will put aside the question 
of the explicit direction of n and concentrate on defining the optimum X or » 
alternatively, the optimum d, The optimum value of d is determined by maxi- 
mizing h with respect to a, i. e. , by solving for the root of the equation 


Bh - 

v , oCr 
__ = h — = o 
od cr od 


(15) 


Performing the indicated differentiation yields 


or 


£2^3- a + X =o 

v r s 


y * 


X v 
s 


( 16 ) 

( 1 ?) 


where, using (9), 


and 


4 

y* = ) a. a + T) 

L 4 
i-0 



a = A # e, 
r t 



(18) 


i/4 

Because of the form of (18) equation (16) is a quartic in the variable d , and 
is solved by iteration in the program, A more detailed discussion of the solution 
to (16) is given in reference [3] for a specific set of coefficients a . For now, 
assume that the optimum value of d is found from (16), Then the optimum angle 
X is immediately obtained 


8 



( 19 ) 


cos X = d r 

Of course, equation (16) does not take into consideration the fact that d can 
never exceed the inverse square of the solar distance. Consequently, the right 
hand side of (19) may exceed unity, under which condition the program sets 
cos X ^ 1 (i«e., X“0) implying that n is directed along e^. 

If the d that represents the solution to (16) exceeds the upper limit 

allowed for d whether that limit is imposed by the problem or by nature (i.e. , 

, 2 

1/r ), the correct choice for d is that upper-limiting value. Likewise, on the 
lower side, d is physically limited to be non -negative. Thei'efore, a negative 
solution to (16) is disregarded, and d is set to zero which corresponds to 
X = ff/2 (panels oriented edgewise to the sun), and the engines are shut down. 

The precise definition of n has no bearing on the solution of the problem, 
except as it affects d as defined in (1). The appearance of n in the state and 
adjoint equations and the variational Hamiltonian is solely through the density d 
except in the equation for A where n appears explicitly as part of the partial 
Sd/SR, defined following equations (12). Actually, this partial is valid only if 
d is permitted to vary with R. That is, if either d = 0 or d = constant is im- 
posed then dd/SR becomes the null vector and the entire term drops from the 
equation for A. Furthermore, if d is the solution to (16), then the last term 
of A in (12) again drops out because the term in square brackets is the left side 
of (16). Therefore, the only time the term in question remains in the equation 
for A is when d = l/r which corresponds to cos X = e * n = 1 and implies 

n = e . Under this condition 
r 

cUl _ 2_ - 

SR " 3 6 r 

r 

The boundary condition pertaining to the initial degradation time is 

s (t ) = 0 

o 
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Although transversali/ty conditions are discussed in the next section, the trans- 
versality condition associated with degradation time is stated here for the sake 
of completenesso If the final degradation time, s^, is unspecified, 

i 

A = 6 . 
s f 


The initial value of A g is unknown and therefore becomes one of the independent 
parameters of the boundary value problem. From equations (12), it follows that 


X 

s 



X 


g 


and, therefore, 





(X - A 
g f g 


Using the boundary conditions 


A 



s 

o 



it follows that 


Thus, if an approximate value of A is available from a trajectory similar to 

gf 

the one of interest, a reasonable guess of A is readily available, 

s 

o 

2, Housekeeping Power, This program option applies only to solar 
electric propulsion and is currently not integrated with the power degradation 
option. The program therefore does not allow the use of both options simulta- 
neously. The new program input-quantity controlling the housekeeping power 
simulation is DPOW, which is the ratio of housekeeping power to reference power 
and is given the symbol Ap: 


Ap= p h /p 


ref 


The reference power p is defined as the power input to the power conditioners 

YG I 

under reference conditions, (For solar electric propulsion, reference conditions 
are equivalent to having the total spacecraft in a perfect solar- simulator such 


* 


*When the thrusters are at full -throttle. 
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that the simulated solar distance is one astronomical unit, the solar arrays are 
oriented normal to the source of energy, and no radiation damage or other de- 
gradation effects have occurred to any components). 

This spacecraft model assumes that the powers required by the operating 
components at each instant of time exactly match the power developed by the solar 
arrays: 


where p is the power developed by the arrays, p is the instantaneous power 
a 

delivered to the power conditioners, and p is the housekeeping power, which 
is constant with time. Currently, all trajectories generated by the program 
using the housekeeping power simulation must satisfy the condition that p > p^> 
so that p remains positive; in other words, engine shutdown when p -* 0 is 
not coded into the program. 

The power factor y is defined as the ratio of instantaneous power to 
reference power, 

y = p/P ref = ( 1 + ^p) l a t d 


where the summation quantity L is the power factor when p^ - 0. 

The spacecraft model has been expanded by deleting the old, total pro- 
pulsion s}rstem specific mass Qi and replacing it with a , the specific mass 
of the solar arrays, and a , the specific mass of the thruster subsystem. The 
quantity previously called the "propulsion system mass” is now expressed, 


m =ol p 
ps a ao t 


P ref 


<“t + « a (1 + Ap)) p ref 


where p is p evaluated under reference conditions (i. e. , p -p + p ), 
ao a ao rei u 

and the total specific mass ol has been replaced by (a + ol (1 + Ap)), also a 

t 3, 

constant since the housekeeping power option is intended to be used only in 
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simulations where the reference power is specified. The remainder of the 
analysis describing the optimization problem remains unchanged. 
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B. 


BOUNDARY AND TR A NS VERSA LIT Y CONDITIONS 


The basic boundary conditions and transversality conditions are described 
in [l] ; the nomenclature used in the following analysis is also described in [l] , 
except for the introduction of new symbols which are described in the text as they 
appear. 

In summary, when no intermediate -targets are present, individual trans- 
versality conditions are derived from the general equation 

• n 

d77+ [A* dV- A * dR + A dt^ + A ds+A dr + A dg+A dc-h dtj =0 

PsTg cvo 

where 77 denotes the performance index which is assumed to be of the form 


77“ 77 (v 


»o 


V “n’ V g > C> * 0 ’ V 


such that 


where 


d77~77 dv +77 
v ®o v 
«o “n 


77 denotes B77/3 x, 


dv +77 d^ + 77 dff + 77 dc + 77 dt + 77 dt 
a> n V ng c totn 

n on 

The individual transversality conditions are obtained 


x 

by equating to zero the coefficients of all independent differentials. 


X. Deep Space Burn , In simulations of trajectories which are all-ballistic, 
the program is capable of simulating a single deep-space burn, or impulsive 
velocity-change, at any point prior to arrival at the primary target. The three 
components of the incremental velocity AV are independent variables of the 
boundary value problem, such that, at a specified time, the spacecraft velocity 
is incremented: 

0 *> 

R = R + AV 


The use of this program option is described in the Sample Problems and Results 
section under sample Case C, Multiple Ballistic Swingby Mission, 
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2 0 V,, Optimization in LVI Mode, The optimization of the launch excess 
velocity V when using the Launch Vehicle Independent (LVI) mode of simulation 

CO 

is accomplished when the initial primer vector is forced to vanish: 


A - Oo 

o 

This is accomplished by setting the values of these three independent variables of 
the boundary value problem to zero and turning their triggers off; the three com- 
ponents of the departure heliocentric velocity become independent variables in- 
stead. The program contains special logic to circumvent the numerical singularity 
associated with the null primer vector. 

3. Launch Asymptote Declination Optimization . This section contains the 
equations employed in the optimization of the launch asjunptote declination. A more 
thorough exposition of this subject is found in [ 4 ] , in which some of the ramifica- 
tions of the theory are discussed. Historically, an incomplete form of this analysis 
was first published in the appendix of [l] . 

The assumed spacecraft and propulsion system models are as described in 
[l] and will not be repeated here. The launch vehicle payload capability is 
assumed to follow the simple exponential law 

m =b e - ( VV-h , (20) 

o 1 3 


where b , b . and b are pre-determined 
12 3 

v^ is a characteristic speed representative 
specific escape trajectory. 


constants for each launch vehicle and 
of the energy required to achieve a 


The velocity penalty incurred with non-due-East launches from the ETR 
is shown graphical^ in Reference [5] as a function of the parking orbit inclina- 
tion. This velocity penalty Aw is adequately approximated with a quadratic 
curve fit of the form 


Av. = c 1 
l 1 


.2 


+ C 2 i + C 3 


( 21 ) 
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Given a reference orbit inclination i and a circular orbit speed , 

the velocit 5 ' penalty Av associated with a non-coplanar departure from this 

circular orbit to the desired hyperbolic excess velocity at a declination <5 is 

defined as follows . Assuming the line of nodes of this reference orbit is an open 

variable, one may choose this variable to minimize the angle between the excess 

r 6 j 

velocity and the orbital plane. This minimum angle is 6 - i. Gunther has 
shown that the minimum incremental velocity required to achieve a given v^ 
along an asymptote not lying in the orbital plane from a specified circular orbit 
is obtained from the solution to a quartic equation in the sine of the out-of-plane 
angle . Defining 

s = sin(<5-i); p-v /v ; 

2 2 

P = s (p +4) ; 

2 _ 2 2 
q - s CL-s ) p ; 

1 1 

x=[y<q/2) 2 + (p/3) 3 + q/ 2] 0 - [M) 2 + (p/ 3) 3 - q/ 2] 0 ; 


y = ,/p 2 /4-x ; 


( 22 ) 


w=|[p/2+y +.,/(P/2 + y) 2 + 4(x/2+ yx 2 /4 + s 2 ) ] , 

then Gunther's solution for the magnitude of the minimum velocity impulse re- 
quired to accomplish the maneuver is 


v = v 
g o 



+ 3-2 



+ p w-w z )(2 +Pw) 


(23) 


and the penalty Av is the difference between 

o 

quired if the out-of-plane angle were zero, i. e. 


v 

g 


and the velocity increment re- 


Av = v 
g g 


-(f 


2 „ 2 
v + 2v 

00 o 


- v 


)• 
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Thus, the definition of the characteristic speed for those cases in which the 
asymptote declination lies outside the interval C-i, i] is 


v 

c 


J + 2v iJ + Av + Av 

V co o i g 


-v + v + Av.. 

O g L 


(24) 


The transversality conditions which optimize the launch phase are deve- 
loped as follows . The state and adjoint equations for the problem under considera- 
tion are precisely as formulated in [lh The only difference in the optimality 
conditions is the format and content of certain of the transversality conditions. 
Specifically, these differences are due solely to the new definition of v^ which 
is now a function of the direction of as well as its magnitude. Whereas be- 
fore the differential of v was, simply 


dv 

c 



it follows from (24) that the equivalent formula now is 


dv =(av /dv 
e g 


)dv + (dv /d 6 ) d 6 

> 00 g 


+ (dv /d i+ d Av /d i) d i 
g i 


where, from (21) 

d Av. /d i = 2 c, i + c_, 

i 12 

and, from (22) and (23), 

d v /d i = - dv /d 6 . 
g g 

The derivation of the partial derivatives dv /dv^ and dv /o6 is straight - 

g g 

forward although somewhat cumbersome. The equations for these partial 


16 



derivatives are as follows: 


2 2 2 2 
3v v w(3 + 2pw-w )(3p/Bv )+(3p + 2p w-3pw -4w) (°w/^v o ) 

p -IB 

^ V °° V g ^ V °° 2 / (1 + pw - w^) (2 + pw) 


2 9 2 

v ( 3 p+ 2 p" 1 w - 3 pw - 4 w) 
o 

2 v J( 1 +pw-w ^)(2 + pw) 
g' 


Bw 
36 ’ 


Bp/3v = l/v } 
» 0 


3w = 1 (I _3£ + Ji + fA 

5 v 2 ^ 2 d v B L \ 


CO CO 


, X ) -*x + ( P )( 1 Sp 

/ 2,: ”2 ' 3v o= 2 ' 2 3v 

2 .J x /4 + s “ 


) *] / ( 2 w- p /2 - y)J , 


f3x/36) + 4s(3s/56) ./ £ . w \ 3* 


2,/ x 2 /4 + s 2 


( | + y) ||]/( 2 w~p/ 2 -y)} , 


9 s/5 <5 = cos (6 - i). 


By/Bv^ | (p/ 2 ) dp/Bv^ - Bx/Bv^] / 2 y , 


Sy/3 6 = - (3 x/3 6)/2y, 


x = l r(q/ 2 )( 3 q/ 3 u) + (p/ 3 ) (B p/Buj + jq, T I {q/2) 2 + fa/a> 3 +q/2 T 
U ° L ,/(q/ 2) 2 + ( p/ 3) 3 5UXV J 


ui M 


(q/ 2) 2 + (p/ 3) 3 


-f^I/rf^ 3 - q/2 ] 


17 



w ith u - v or 6 , 

CO 

dq/3v w - 2ps 2 (1 -s Kdp/dv^) , 
dq/d 6 ~ 2p 2 s(l - 2s" 1 ) (3 s/d 5 ) , 
dp/dv co -2ps 2 (dp/d Voo ), 
dp/d 5 = 2s(p 2 +4)(S s/d 6). 

This concludes the computation of these two partial derivatives. Some symbols 
used in this section of the document apply to this section only; for example, p 
and q here are intermediate quantities in the analysis, and have no relation to 
power or degradation factor discussed elsewhere. 

The differential of g, the reference thrust acceleration, is 


d g “ - (g/ m Q ) d m 

where 

dm =-(bVb_)e ^c^^dv 
o '12' c 


and the differential of may be written 


dV = 

CO 


(V /v 

' CO CO 


) dv „ + < n p x v co )da+[(V co x lyxvylv 


x n I ] d5 , 

CO p 1 


where n is a unit vector along the North Pole and a is the geocentric right 

ascension of V , and therefore the transversality conditions associated with 
00 

V are as follows: 

CO 


For optimum launch parking orbit inclination: 


f(d Av. /d i 


d v /d 6 ) = 0. 
g 


( 25 ) 
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For optimum launch excess speed: 



For optimum launch asymptote declination: 

f(&v /d6)-A e [(V xn )xV / 1 V xn 

' g ' O 03 P co ' co p 

For optimum launch asymptote right ascension: 


] = 0 . 


( 26 ) 


(27) 


- A • (n x V ) - 0 . (28) 

o p “ 

In the above equations f is given by 

f=[k +k -(1+kJP -gX /m ] dm /dv 
s t r f go oc 

Due to the complexity of the equations defining ^v^/d^, the solution of equation 
(25) for the optimum value of i must be obtained using an iterative technique. 

The approach to the solution of the problem as formulated above differs 
in three basic respects from that of the problem where asymptote declination is 
ignored: (X) the condition (25) must be solved for the optimum parking orbit in- 
clination, given values of v^ and 6 ; (2) the asymptote declination 6 must be 
introduced as an independent parameter and (27) added as an end condition of the 
problem; and (3) the evaluation of V a becomes somewhat more involved. The 
computation of given v^, and 6, proceeds as follows: Denoteas C 

the obliquity of the ecliptic such that the matrix 




10 0 
0 cos € -sin € , 

0 sin € cos 6 
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operating on a vector expressed in ecliptic Cartesian coordinates yields the same 
vector in Earth equatorial coordinates. Then the right ascension of the initial 
pl'imer may be written 

a x = tan ^ F(X cos <r - X sin e )/X , 

X L v yo zo xo J 

where X , X , X are the given ecliptic coordinates of A • Then, the right 
xo yo zo u 

ascension of the asymptote is set 


a - a 


X 


or 


oc = + 


TT 


and V is evaluated 

CO 


V 

CO 



cos cl cos 6 
sin OL cos 6 
_ sin 6 _ 


4 0 Extra-Ecliptic Missions. An additional set of boundary conditions has 
been added to the program for simulating extra-ecliptic missions. 

The original set of boundary conditions is given on pp. 25-26 of reference 
[l] , which is now invoked by employing the program input quantity IOUT = 1 and 
corresponds to i, e, and r^ specified, ^ and CO optimized and f ~ 0, where 
i is inclination to the ecliptic, e is eccentricity, r^ is perihelion distance, S7 
is ascending node angle, CO is argument of perihelion, and f is true anomaly, 
all evaluated at the final time. 

The new set of end conditions for extra-ecliptic missions is invoked by 
IOUT = 2 and is similar to the original set except that the final true anomaly is 
optimized rather than specified as zero, which corresponds to having the space- 
craft arrive at the optimum point of the extra-ecliptic orbit rather than being con- 
strained to arrive at perihelion. The specified final conditions are i, e, and a, 

where a is the final semi-major axis (in place of r )„ r2 and co are optimized 

P 
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as before, using the transversality conditions 


T(ft) = k • C = 0 

(29) 

T(w) = h • C = 0 

(30) 


which are identical to equations (47) and (48) of reference [l] . C=RxA-R x A 
the vector constant of the motion on a given trajectory segment, h is a unit vector 
along the angular momentum of the final orbit, and k is the z-axis unit vector = 
Optimizing f is accomplished by satisfying the transversality condition asso- 
ciated with f; 



(A • K) + r (A'R) = 0 


(31) 


The program attempts to satisfy equations (29), (30), and (31) directly, since no 
closed-form solution has yet been found for casting these end conditions in the 
form of AR and AR similar to equations (54) and (55) of reference [l] . 
These end conditions possess poor convergence properties, a fact which the 
authors have experienced directly in attempting to obtain converged solutions 
using IOUT - 2„ 

5c Multiple Target Missions - Errata , Two equations in the section, 
Extension of Formulation for Multiple -Target Missions of reference [l] are 
in error o Equation (116) of reference [l] should read 


77 
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samp i 



and equation (117) should read 
B 77 
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c. 


AUXILIARY COMPUTATIONS 


This section presents equations employed in computations which are 
made after the iteration sequence involving the primary target is completed. 

1. Additional Block Print Variables. A standard print block is employed 
for printing information at various points along a trajectory. Each standard 
block contains a total of forty parameters, which are described in [l] . 

The standard block may be augmented in two ways. When power degra- 
dation, as indicated by the input variable TPOWER, is simulated, a single line 
of information is automatically added to each block, as displayed in the Sample 
Problems and Results section. Case B, Comet Rendezvous Mission. When the 
input variable MPRINT is 2 or 3, three extra lines of information are generated 
per block, as displayed in the Sample Problems and Results section, Case C, 
Multiple Ballistic Swingby Mission, These two types of additionally printed lines 
may appear simultaneously. 

(a). Power Degradation. The single line of power degradation 
information contains eight parameters as follows: 

Degradation time, s, since departure, in days, 
given by expression (4). 

Degradation time adjoint variable, 

Density parameter, d, used in expression (9) for 
the power ratio y, and given by expression (1), 
in AU’ 2 . 

3 y/3 r, in AU ^ . 

3 y/dd. 

DEGRAD The degradation factor, q, given by expression (6). 

CHI Solar array orientation angle X as given by, for 

example, expression (19), in degrees. 


S 

LS 

DENSITY 

DPOWR 

DPOWD 
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CHI REF 


Solar array orientation angle which the arrays 
would have if oriented for maximum power, in 
degrees. 


(b). Target-Relative Coordinates and Comet Magnitudes. The 
three extra lines which may appear via using MPRINT contain the following in- 
formation: 


R1 REL 
R2 REL 
R3 REL 


VI REL 
V2 REL 
V3 REL 


Cartesian components of current spacecraft 
position vector, with respect to the next astro- 
nomical body to be encountered along the 
trajectory in a moving coordinate system generated 
by that body, with the x-axis pointing outwai-d along 
the body f s heliocentric radius vector, the y-axis in 
the body T s orbit plane in the direction of the body f s 
motion, and the z-axis completing the right-handed 
orthogonal system, in kilometers, with the origin 
of coordinates at the body. 

Cartesian components of current spacecraft velocity 
vector, in kilometers/second, in the target-relative 
coordinate system described direct!}' above (see 
R1 REL), 


RMAG REL 


Magnitude of R1 REL, R2 REL, R3 REL, in kilo- 
meters . 


VMAGREL 


Magnitude of VI REL, V2 REL, V3 REL, in kilo- 
meters/second. 


S/C NUC MAG 


Nuclear magnitude (of comet) of the next astro- 
nomical body to be encountered along the trajectory, 
as seen by the spacecraft. 


M = M + M 
N o 1 




+ M 

2 




+ o 03 cos 


-1 


R, * (R, - R) 

targ targ | 


R 


targ 


|r 


targ 


-R 


where M q} M^, and M 9 are. magnitude constants 
( Continued on next page) 
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S/C NUC MAG 
(continued) 

S/C TOT MAG 


GEO NUC MAG 
GEO TOT MAG 
ANG(V, R) 

ANG(V, XY) 


R1 REL ECL 
R2 REL ECL 
R3 REL ECL 

VI REL ECL 
V2 REL ECL 
V3 REL ECL 

R3VIAG ECL 


VMAG ECL 


associated with the target, and C° is the radians- 
to-degrees conversion factor* The arc-cosine term 
is the phase angle. 


Total magnitude of the next astronomical body to 
be encountered along the trajectory, as seen by the 
spacecraft 


M T = M 3 + V° g 10 


|E - R | 
targ 


5 ^10 


R. 


targ 


where M Q , M., and M 5 are magnitude constants 
associated with the target. 


Shine as S/C NUC MAG, except as seen by the Earth. 
Same as S/C TOT MAG, except as seen by the Earth. 


Angle which (VI REL, V2 REL, V3 REL) makes with 
the positive x-axis in the target-relative coordinate 
system described under R1 REL, in degrees, 

ANG(V, R) - cos" 1 (VI REL/VMAG REL) 

Angle which (VI REL, V2 REL, V3 REL) makes 
with the xy plane in the tar get -relative coordinate 
system described under R1 REL, in degrees. 

ANG(V, XY) = sin’* 1 (V3 REL/VMAG REL) 


Same as R1 REL, R2 REL, R3 REL except ex- 
pressed in the ecliptic coordinate system of 
date. 


Same as VI REL, V2 REL, V3 REL except ex- 
pressed in the ecliptic coordinate system of 
date. 


Magnitude of R1 REL ECL, R2 REL ECL, 
R3 REL ECL, in kilometers. 


Magnitude of VI REL ECL, V2 REL ECL, 
V3 REL ECL, in kilometers/second. 
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2, Swingby Continuation Analysis, Auxiliary computations are optionally 
provided, invoked by the NAMELIST input vector MOPT4, whereby ballistic 
swingbys past the primary target may be simulated. 

In one mode of program operation, invoked by MOPT4(l) > 0, single 
swingbys past the primary target may be simulated to up to ten post-swingby 
targets per case. 

In another mode of program operation, invoked by MOPT4(l) < 0, multiple 
swingbys along a single trajectory may be simulated, first swinging past the 
primary target and then subsequently swinging past more targets downstream 
along the trajectory. One multiple swingby trajectory may be simulated per 
case. 

In either mode of operation, the following basic assumptions are made. 

The swingby continuation computations are independent of the trajectory leg 
leading up to the swingby target, which may consist of an optimized electric 
propulsion trajectory segment (if the swingby target is the primary target), 
except that the arrival V m and arrival time at the swingby target are used in the 
determination of the swingby passage conditions. Each swingby maneuver is 
calculated under the assumption of the patehed-conic approximation, and the 
swingby planet's sphere -of- influence is assumed to have zero radius as seen 
from interplanetary space and infinite radius as seen from the planetary vantage 
point. The passage time in the swingby planet's sphere -of- influence is neglected 
(taken to be zero in the heliocentric frame). 

Each swingby maneuver may be either unpowered or powered, and these 
two cases are discussed in the following sections. Since the unpowered swingby 
solutions are embedded in the wider class of powered- swingby solutions, tending 
to appear in pairs which are separated by a region of braking powered swingbys, 
the more general case of powered swingbys is discussed first. 
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(a). Powered Swingbys. This type of swingby maneuver is re- 
stricted to occur at the mutual perifoci of the approach and departure hyperbolic 
arcs; the powered phase is impulsive and the thrust is colinear (pro or con) to 
the velocity at closest approach. Whether the swingby is powered or unpowered, 
the trajectory segment leading up to the swingby planet has been pre-determined, 
this being the method by which the program has been designed to obtain swingby 
solutions. Therefore the swingby time and the arrival hyperbolic excess velocity 
V are known. In the following analysis, subscript A pertains to arrival at 
the swingby planet and subscript D pertains to departure. 

A basic assumption of the powered swingby problem posed here is that 
the flight time from the swingby planet to the next target is specified. This being 
so, the program is able to converge, by iteration, on some ballistic trajectory 
from the swingby planet to the next target having the specified transfer time, im- 
plying that the departure hyperbolic excess velocity V at the swingby planet 
is thereby determined. Therefore, the heliocentric trajectory before and after 
the swingby planet is determined, and it then remains to perform the required 
computations pertaining to the hyperbolic arcs within the sw r ingby planet's sphere 
of influence. 


and 


The closest approach distance is found by iteration as follows. Let 


2 

r v * 
„ , p "A 

O' = 1 + — L 

A n 



1 + 


r v 


2 

<»D 




where v^ - | V c&A [ , v^ = (v^l, is the (unknown) passage distance, 
and (j, is the swingby planet's gravitational parameter. Then the approach and 
departure hyperbolic bend angles are given by 
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A -1 . -1 . , 

•— = cosec a =sm (1/a.) 
2 A A 


D -1 . -1 „ , . 

— = cosec a = sin (1/a ) 

la iJ X.J 


and these must sum up to the total bend angle, which is specified in terms of 


V . and V 
50 A 00 D 


6 6 

* A D -1 

6 = - = cos 

T 2 2 


- V co A - V ~D 1 
U V V _ J 
00 A 


Therefore, using r as the independent variable, the zero of the quantity 


-1 -1 -1 

F = sin (l/a ) + sin (l/a ) - cos 
a n 


_ V 


A - V „ 
»A “P 

v v ^ 
“A 06 D 




is obtained by Newton's iteration, using the derivative 

aF _ / -IN r “A x A 

3r \ u / L r~z 


/ 


2 / 

Vco T) a T> 


a 2 - 1 
D 


When the iteration is converged, the passage distance r is in hand, and the 
impulsive velocity increment is computed. 



where the square-root -quantities are the hyperbolic speeds at closest approach., 
The remaining parameters defining the planetocentric transfer are computed as 
follows. The inclination of the swingby orbit plane to the planet's equator is given 
by 

i = cos ^ (h - n ) 

P 


where h is the unit vector along the angular momentum of the hyperbolic passage 
trajectory and n is a unit vector pointing toward the swingby planet's north pole. 
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The ascending node angle of the swingby orbit plane is computed as 

Q - tan * (-h /h ) 
x y 

and is placed in the proper quadrant by using the system library routine DATAN2. 
The argument of perifocus is given by 

— X A A 

to = cos (r ° r ) 
p n 

A A 

where r is the unit vector pointing toward the closest approach point and r^ 
is the unit vector lying along the line of nodes and pointing toward the ascending 
node. This is adjusted for the proper quadrant by the test, 


A A 

Ifh (r xr) <0, oo -* 2 77 - co 
z n p z 

In the right-handed planetary reference frame, the z-axis is toward the planet's 
north and the x-axis points toward the ascending node of the planet T s equator on 
the ecliptic plane. 


(b). Unpowered Swingbys. This type of swingby maneuver is 

considered to be a powered swingby having Av = 0. The program adjusts the 

post-swingby heliocentric trajectory segment, by iteration, until the swingby 

departure V magnitude equals the given arrival V magnitude. The primary 

independent variable in this iteration is the post-swingby transfer time to the 

specified target, which was held constant in the powered swingby case. Thus 

v s= y = v , and the swingby passage distance is obtained from the 
cod "A 00 

formula, 


r = 
P 



2 v 


Iv - V 

1 “A “D 


- 1 


) 


The other orbital parameters are obtained from the same relations given above 
in the section, Powered Swingbys. 
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The program can generate multiple-revolution ballistic arcs, and a 
particular solution obtained by the program may not be unique, even for the 
same transfer time, All solutions are reachable, however, by means of in- 
putting an appropriate initial velocity guess for the trajectory segment in 
question. 


30 



III. 


PROGRAM INPUT 


The following consists of a complete description of program inputs. 
With respect to the basic HILTOP report [l j , several input quantities have 
been added, some have been deleted, and some definitions have been modified. 
The input quantities which have been deleted are CALPHA, DIGIT, HALT, and 
WEIGHT, and so these names must no longer appear in NAMELIST MINPUT 
input data sets. CALPHA has been replaced by ALPIIAA and ALPHAT. Also, 
Yi(4) and Yi(5) for i = 1, 2, 3, ...» 70 have been deleted (See the description 
of "dependent parameters" in the section just ahead, Definitions of Input 
Parameters) and should therefore also be removed from input data sets. 
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A. NAMELIST 

Inputs to HILT OP are given through the NAMELIST feature of the IBM 
Fortran IV programming language. The input NAMELIST is named MINPUT, and 
every input required or used in the program is declared by name in the list. The 
general form for assigning an input value to a quantity is, simply, 

NAME - VALUE 

where NAME is the name assigned to the variable and is included in the NAME- 
LIST, and VALUE is a numerical or logical quantity consistent in form (i.e. , 
logical, integer, or real) with NAME. Unless otherwise specified, all MINPUT 
names commencing with one of the letters I through N represent integers, 
whereas all names commencing with one of the letters A through H or O 
through Z are double precision floating point numbers. Each NAMELIST case 
must begin with the characters 

& MINPUT 

commencing in card column 2 and followed by at least one blank, and end with 
the characters 


&END 

preceded by at least one blank. Card column 1 is ignored on all NAMELIST 
input cards. Multiple data assignments on a single card are permissible if 
separated by commas. Blanks in the variable field, VALUE, are taken as 
zeroes. A comma following the last VALUE on a card is optional on the IBM 
system. The order of the input data assignments is arbitrary; i.e D , they need 
not be in the same order as listed in the NAMELIST. In fact, there is no re- 
quirement that any specific input parameter be represented in the input data 
set. If no value ls included in the inputs for a particular parameter, the de- 
fault value is used (see Default Values). For other details regarding the 
NAMELIST feature, the reader is referred to the IBM System 360/Fortran IV 
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Language manual. NAMELIST cases may be stacked back-to-back indefinite^, 
A single NAMELIST error majr wipe-out the remaining NAMELIST inputs. 



B. DEFINITIONS OF INPUT PARAMETERS 

Specific examples of the program inputs are given In the Sample Problems 
and Results section. Default-values of inputs are given in the next section. 

The program inputs, in alphabetical order, are: 


AAI Desired final extra-ecliptic inclination. Related to 

AE, AR, and IOUT. [deg] 

AE Desired final extra-ecliptic eccentricity. Related to 

AAI, AR, and IOUT. 

ALPHA A Specific mass of solar ar raj's, a . Ckg/kwJ 

a 

ALPHAT Specific mass of thruster subsystem, a . Ckg/kw] 

ALTITU This input variable is associated with program logic 

which has not been kept up-to-date, specifically, 
logic pertaining to optimum departure of a NERVA- 
type rocket from Earth orbit. This variable should 
be ignored. 

AN Trajectory- integration exponent in expression (37) of 

reference [lj. 


AR Desired final extra-ecliptic perihelion distance. Re- 

lated to AAI, AE, and IOUT. [AU] 


ASOL 


BI 

B1 ' 
B2 > 
B3 


Array of five elements consisting of the solar power 
law r coefficients in expression (17) of reference 
[l] . ASOL(l) > 0 tells the program to use the input 
coefficients rather than the internal coefficients. The 
coefficients are normalized internally, and the program 
executes the iterations to produce the required remark- 
able points of the power curve (which are printed). 

Efficiency coefficient b in expression (15) of reference 
[l] . Related to DI and El. 

Launch vehicle coefficients b-j, b 9 , and bg in 
expression (2) of reference Cl j . Used only if MBOOST 
is negative, [kg, m/sec, kg] 
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CNI 

CNEX 

CSTR 

CTANK 

CTRET 

DI 

DMRETR 

DPOW 

eci 

EC IX 
El 


Inclination to ecliptic of primary-target orbit. Input 
only when MOPT3 = 11. Related to ECI, OMI, SAI, 

SOI, TPI, EMUODD, and RADODD. [deg] 

Array of five elements, the first three of which may 
be currently used. Inclinations to ecliptic of inter- 
mediate-target orbits. Input CNIX(i) only when 
MOPTX(i) = 11. Related to ECIX, OMIX, SAIX, 

SOIX, TPEX, EMUODX, and RADODX. [deg] 

Structural factor, k , in expression (7) of reference 

[i]. s 

Propellant tankage factor, k^_, in expression (6) of 
reference [l]. 

Retro tankage factor, k , in expression (10) of reference 

Ci]. rt 

Efficiency coefficient d in expression (15) of reference 
[l] . Related to BI and El. [km/sec] 

Retro engine mass, m rs , in expression (10) of reference 

[i3 . Ck g ] 

Ratio of housekeeping power p to reference power P re f» 

The power transmitted to the propulsion system is that 

generated by the arrays less housekeeping power which is 

constant along the trajectory. The power output of the 

arrays normal to the sun at 1 AU is p „ + p, . 

ref h 

Eccentricity of primary-target orbit. Must be less than 
unity. Input only when MOPT3 = 11. Related to CNI, OMI, 
SAI, SOI, TPI, EMUODD, and RADODD. 

Array of five elements, the first three of which maybe 
currently used. Eccentricities of intermediate-target 
orbits. Input ECIX(i) only when MOPTX(i) =11. Related 
to CNDC, OMIX, SAIX, SOIX, TPDC, EMUODX, and 
RADODX. 

Efficiency coefficient e in expression (15) of reference 
[l] . Related to El and DI. 
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EMUODD 


EMUODX 

GAMMAX 

GAP 


HOUR 

IBAL 

INTPR 

iout 


Gravitational constant of primary -target. Input only 
when MOPT3 = 11. Related to ECI CNI, OMI, SAI, 

SOI, TPI, andRADODD. [m 3 /sec 2 ] 

Array of five elements pertaining to the gravitational 
constants of intermediate -targets. These inputs must 
be ignored at present. 

Maximum permissible value of the power function y 
when MODE = 5. At solar distances less than the value 
for which y~ GAMMAX, the solar arrays are assumed 
to be tilted such that y is maintained at the limiting 
value. 

Propulsion-corner proximity tolerance -interval, Act. 

See the discussion in the section Avoiding Corners in 
the Propulsion-time Function in reference [l ] . When- 
ever the thrust switch function 0“ grazes the zero-axis 
within the tolerance |Acr | on any trajectory, an internal 
counter is incremented, and the trajectory is considered 
to be in the neighborhood of a propuls ion -time corner. 
Positive value of GAP causes forced-thrusting case to 
be inserted, negative value causes bypass to next case, 
whenever the internal counter reaches the related input 
variable NHUNG. 

Hour-of-day of reference date (e.g. , 17.342D0). 

Related to MYEAR, MONTH, and MDAY. 

Ballistic option indicator. Setting IBAL 0 invokes 
option 1 discussed in the section Ballistic Trajectory 
Option of reference [l] . 

Indicator which specifies print-length when the iteration 
in subroutine INTERP fails. Value of 0 causes short- 
print and 1 causes detailed-print. 

Extra-ecliptic mission indicator. IOUT - 1 or 2 
indicates that extra -ecliptic target conditions are 
desired, in which the iterator dependent variable 
triggers Yl(2) through Y6(2) are set equal to 1, 
and for which the input LAUNCH (which see) 
should probably be set to 1, and parameters related 
to LAUNCH also set appropriately. Ordinarily MOPT2 = 3. 
No retro stage may be employed. (Continued on next page). 
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IOUT 

(cont) 


ERK 

ERL 


moT 


ISPIN 

ITF 


IT PR NT 
JPP 


= 1 i, e, specified; f = 0. 

= 2 i, e, a specified; f optimized,, 

In the above, i = final extra-ecliptic inclination, 
e = final eccentricity, = final perihelion distance, 
a = final semi-major axis, and f = true anomaly at the 
final time. Final Q and (X> are optimized in both 
cases. Related to AE, AR, and AAI. 

Numerical integration option (currently not used). 

Primer-origin-proximity step-size-control indicator. 
Value of zero causes bj^pass of control, leaving the 
step-size Au constant. See discussion in the section, 
Integration (Thrust) of reference [l] . 

A non-zero value of EROT causes the input ecliptic 
projection of the primer vector and its time derivative 
to be rotated about the z-axis through an angle equal 
to the difference in longitudes of the spacecraft between 
the last trajectory of the previous case (or zero if no 
previous case) and the first trajectory of the current 
case. This feature permits one to use the initial ad- 
joint variables from a 2 -dimensional trajectory as the 
initial-guess inputs for a 3 -dimensional trajectory using 
the ephemeris option. EROT must be set back to zero to 
avoid undesirable rotations on subsequent cases. 

Spinner indicator. Not used at present. 

Provides normal termination conditions for runs which 
require more machine time than is estimated. The 
value specifies the number of machine-time seconds 
(CPU and I/O) required to execute the summary tra- 
jectory after halting the iteration-sequence, [sec] 

Does not apply if subroutine REMTIM is dummied. 


Indicator for special print from MINMX3 iterator. 
Non-zero value invokes print. 


Jettison indicator j^ g for electric propulsion system 

prior to primary-target retro -maneuver, as used in 
expression (8) of reference Cl]. (Continued on next 


page). 
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JPP 

(cont) 


J PRINT 


JT 


KPART 


LAUNCH 


LOADX 


- 0 Propulsion system not jettisoned 
= 1 Propulsion system jettisoned prior to retro 
maneuver. 


Unit 11 printout-length indicator. A value of zero 
causes the iterator independent and dependent variables 
to be output only for each summary-trajectory; a value 
of one causes the same output additionally at each 
iteration of an iter at ion -sequence. 


Jettison indicator j for electric propulsion tankage 

prior to primary-tai'get retro -maneuver, as used in 
expression (8) of reference [l] . 

= 0 Tankage not jettisoned 

= 1 Tankage jettisoned prior to retro -maneuver. 


Option for automatically selecting improved independent 
parameter perturbations for generating the iterator's 
partial derivative matrix. The option is invoked by 
setting KPART = N (N > 0), where N is the maximum 
number of allowed steps, as discussed in the section, 
Perturbation Step Size Selector of reference [l] . 
KPART must be set back to zero if not desired on 
subsequent cases. 


Launch mode selector, pertaining to the optimization 
of the departure asymptote declination, invoked by 
LAUNCH = 1. Related to X10, Y10, X17, and 
Y17. 

Intermediate -target initial-guess feature. Should be 
used with NSET(5) = 1, and then set to zero on the 
subsequent case, A non-zero value of LOADX will 
invoke this feature, whereby the primer A and its 
derivative A will be loaded into the iterator inde- 
pendent-variable arraj^s at each intermediate -target 
(Continued on next page). 
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LOADX 

(cont) 

MAXHAM 

MEOOST 


MDAY 

MODE 

MONTH 


provided that the trigger of the independent variable 
is on* The sole purpose of this capability is merely 
to generate an initial-guess for a multiple -target 
inission, where the values loaded into^the iterator 
arrays represent continuous A and A at each target. 

Maximum number of times that the program will print 
the warning message BAD HAMILTONIAN on any given 
computer run. 

Launch vehicle selector. 

= 0 ATLAS (SLV3X)/CENTAUR 

1 TITAN III C 

2 TITAN III C (1207) 

3 TITAN III X/CENTAUR 

4 TITAN m X (1207) 

5 TITAN IIIX (1207)/CENTAUR 

6 SATURN IB/LM 

7 SATURN IB/CENTAUR 

8 SATURN IC/SIVB/ CENTAUR 

9 TITAN III X (1205)/CENTAUR 

10 TITAN III B (CORE)/ CENTAUR 

11 TITAN III D (1205)/ CENTAUR 

12 DELTA 

13 TITAN III D 

14 TITAN III D (1205)/ CENTAUR/TE364 (2250) 

15 TITAN III E/CENTAUR 

16 SHUTTLE /TRANSTAGE 

17 SHUTTLE /DELTA 

18 SHUTTLE/AGENA 

19 SHUTTLE/CENTAUR 

20 SHUTTLE/CENTAUR/BURNER II (2300) 

NEG Use input booster coefficients b , b , and b . 

12 3 

Day-of-month of reference date (e. g. , 26). Related to 
MYEAR, MONTH and HOUR. 

Power variation option selector. The value of MODE is 
equal to the option-number of the power -curve, discussed 
in the section, Electric Propulsion System of reference 
Cl] (which see). Possibly related to ASOL and GAMMAX. 
MODE = 1 has been eliminated. 

Month-of-year of reference date (e.g. , 8). Related to 
MYEAR, MDAY, and HOUR. 
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MOPT Ballistic option indicator* Using MOPT invokes option 2, 

discussed in the section, Ballistic Trajectory Option, 
of reference ClJ , as follows; 

= 0 No action (use input A , A , and v ). 

x 1 O O “o' 

= 1 Generate ballistic solution with flyby end 
conditions. 

- 2 Generate ballistic solution with orbiter end 

conditions. 

Related to REVS. 

MOPTX Array of five elements, the first three of which may be 

currently used* This array specifies the target-number, 
or planet-number, of the successive intermediate-targets, 
and a value of zero indicates absence of the intermediate - 
target. A zero-entry must not precede a non-zero entry. 
Planet selection is the same as for MOPT2. MOPTX(l) 
pertains to iterator parameters X41-X50 and Y41-Y50; 
MOPTX (2) pertains to X51-X60 and Y51-Y60; and 
MOPTX (3) pertains to X61-X70 and Y61-Y70. Times 
at the targets are X48, X58, and X68. Not to be used 
unless MOPT2 ^ 0. 

MOPT 2 Launch planet number and ephe me r is -option indicator. 

- 0 Analytical planetary ephemeris is not used. 

^ 0 Analytical planetary ephemeris is used and 
the specific launch planet is selected as 
follows: (Continued on next page). 


= 1 

Mercury 

= 14 

Icarus (1978)* 

2 

Venus 

15 

Eros 

3 

Earth 

16 

Geographos (1983)* 

4 

Mars 

17 

Encke (1977)* 

5 

Jupiter 

18 

Encke (1984)* 

6 

Saturn 

19 

Encke (1987)* 

7 

Uranus 

20 

Halley 

8 

Neptune 

21 

Betulia 

9 

Pluto 

22 

Toro (1983)* 

10 

Ceres 

23 

Pallas 

11 

Input Target** 

24 

Juno 

12 

D T Arrest (1982)* 

25 

Vesta 

13 

Encke (1980)* 

26 

Astraea 


* Year -value indicates apparition for which internal orbital elements are 
most accurate. 

**Input corresponding orbital elements (see CNI, CNIX). None are available 
for the launch planet. 
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MOPT2 

(eont) 


-27 Hebe 

28 Iris 

29 Flora 

30 Achilles 
31- Amor 

32 Hidalgo 

33 Alinda 

34 Grigg-Skjellerup (1977)* 

35 Kopff 

36 Grigg-Skjellerup (1982)* 

37 Ganymed 

38 Ivar 

39 Beira 

40 Kepler 

41 Giacobini-Z Inner (1985)* 

42 Borrelly (1987)* 

43 Tempel II (1988)* 

44 Tempel II (1983)* 

45 Tuttle -Giacob ini -Kresak 

46 Schaumasse 

47 Honda-Mrkos- Pajdusakova 

48 Giacob in t-Z inner (1979)* 

49 Icarus (1987)* 

50 Toro (1987)* 

51 Geographos (1987)* 


*Year-vaIue indicates apparition for which internal orbital elements are 
most accurate. 


MOPT3 Planet number of primary-target. Planet selection is the 

same as for MOPT2. If ephemeris is not used, MOPT3 
is used only for retro-stage mass computations. 

MOPT4 Array of ten elements, specifying up to ten post-swingby 

targets. Planet selection is the same as for MOPT2, and 
a value of zero indicates the absence of a post-swingby 
target. A negative value in MOPT4(l) selects multiple 
ballistic swingbys, rather than a set of single swingbys 
in which case also set MAXHAM = 0. Negative values 
(in absolute value) produce planet selection the same as 
for MOPT2. When MOPT4(l) <0, the remaining elements 
of MOPT4(i) may be positive or negative. See the section, 
Swingby Continuation Analysis for details and Sample Case 
C for an example -case. Should be used only for primary- 
target flyby missions. Related to T2, MSWING, NSWING 
and XSWING. 
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MPOW 


MPRINT 


MPUNCH 


Flag used in conjunction with the solar array degradation 
option. Value of zero results In the optimum orientation 
of the arrays relative to the sun line, A non-zero value 
forces the arrays to an orientation yielding the maximum 
power achievable at that instant. Related to T POWER. 

Indicator for printing the summary-trajectory (final 
trajectory of a case) as a function of time or for invoking 
extra printout. 

= 0 Small-size block print at thrust switch points 
only (SWITCH POINT SUMMARY page), 

= 1 Same as = 0, except expands to become a 
standard print -block of parameters for each 
computed point along the trajectory, including 
the trajectory extension controlled by the in- 
put variable TGO. 

= 2 Same as = 0, except each block contains extra 
lines consisting of target-relative coordinates 
and target magnitudes. 

- 3 Combination of =1 and =2. 


Punehed-card and trajectory-tape generation control. 


= 0 No special output. 

1 Punch final values of independent parameters. 

2 In addition, punch selected mission analysis 
parameters used for graphic documentation 
or other purposes. 


< 0 and> -100 Punch trajectory output used with the ASTEA 
program. The absolute value of MPUNCH 
determines the frequency of trajectory points 
output. 


^ -101 Trajectory tape output used with the ASTEA 
program. The absolute value less 100 de- 
termines the frequency of trajectory points 
output. Related to NTAPE. 
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MREAD 


MSWING 


MTMASS 


Card input option (iterator independent variables) 

= 0 No special cards input. 

1 The independent variables generated by 
a previous run by the MPUNCH option are 
input following the NAMELIST case, as 
discussed in the section. Program Output, 
of reference Ci]. 

Array of ten elements, used only when running multiple - 
target ballistic swingbys, such that MSWING (i) corres- 
ponds to MOPT4(i) and selects the type of swingby 
maneuver desired at the respective swingby target. 

Used only if MOPT4(l) < 0. The shooting method 
(MINMX3 iterator) is used, and values of -1, -2, or 
-3 correspond to a swingby passage distance initial 
guess of r = 03 (i.e., continuous heliocentric velocity). 
Each element MSWING(i) may have any of the following 
values : 

= - 1 Go* directly for unpowered swingby; if and 

only if it fails, go for powered swingby having 
flight time T2(i) = initial guess. 

= - 2 Go directly for powered swingby only, having 
T2(i) - flight time, of post- swingby leg. 

= - 3 Go directly for unpowered swingby; then, 

whether it succeeds or not, go for powered 
swingby having T2(i) ~ flight time. 

= - 4 Go directly for unpowered swingby, but 
using initial velocity guess loaded into 
XSWING(j,i), j =1,2,3, similar to 
MSWING (i) = - 1. 

= - 5 Same as = - 2, except use initial guess as 
in = - 4. 

*”Go for” means "attempt to obtain (solution)". 

Related to MOPT4, T2, XSWING, and NSWING. 

Mission-type selector pertaining to the primary target. 

= 0 Flyby mission 

1 Orb iter (high-thrust retro-maneuver without 
velocity loss). 

(Continued on next page). 
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MT MASS 
(cont) 


MUPDAT 

MYEAR 

NDIST 

NHUNG 

NORMAL 


= 2 Orbiter (high-thrust retro-maneuver with 
velocity loss). 

3 Specified arrival excess speed v . 

. If v — 0 , rendezvous mission 
00 n 

If v > 0, controlled flyby mission 
°°n 

No retro “maneuver in either case. 

4 Orbiter (Electric propulsion system performs 
spiral maneuver. Arrival excess speed v 

00 n 

must be specified as zero). 

Other parameters which may be related to MTMASS are 
DMRETR, CTRET, RPER, RAP, THRET, SPIRET, JPP, 
and JT. 

Flag indicating whether iterator independent variables 
at end of one case are to be used as first guesses for 
next case. 

= 0 Use input initial guesses. 

1 Use independent variables of previous case 
as first guesses. 

Year of reference date (e„g. , 1982). Related to MONTH, 
MDAY, and HOUR. 

Identification number of celestial body to be used as the 
reference for the communication distance and angle 
measurement printed in the Extremum Point Summary 
Table. Identification code is the same as for MOPT2. 

Maximum number of propuls ion-corner-proximity 
occurrences allowed in a given iteration-sequence. 
Related to GAP. 

Automatic adjoint -variable scaling. 

= 0 No action. 

1 All A and A are scaled such that A. 
becomes unity. 
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NPERF 


Identification number of end condition that is to be used as 
the performance index when employing the direct parameter 
optimization feature (Improve Mode). The identification 
code is the same as the i in the Yi end condition array. 

NPRINT Print selection flag. Permits selection of amount of 

printout desired on each case. 

= 0 Print only the case summary. 

1 Print switching point summary of final 
trajectory. 

2 Print MINPUT and case setup. 

4 Print trajectory summary on each iteration. 

8 Print partial derivative matrix each iteration. 

Combinations of options obtained by summing options desired. 
If NPRINT > 15, printout consistent with NPRINT = 0 is 
obtained. If the sign of NPRINT is reversed to negative, 
the iterator independent and dependent variables addi- 
tionally are printed for every trajectory which HILTOP 
generates (including neighboring trajectories). 

NSET Iteration-sequence control array. 

NSET(l) Not used for input. 

NSET (2) Not used for input. 

NSET (3) Maximum number of iterations permitted 
in attempting to satisfy (point) constraints 
in satisfy mode. If zero, no upper limit 
imposed. 

NSET (4) Flag indicating whether (point) constraints 

are to be satisfied prior to entering improve 
mode. 

= 0 Satisfy constraints first. 

1 Proceed immediately to improve mode. 

NSET (5) Maximum number of iterations permitted 

after entering improve mode. Setting 
NSET (5) - 1 causes iterator to be bypassed 
and computes single trajectory to obtain 
printout. 
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NSWING 


NSWPAR 


NTAPE 

OMI 

OMK 

POWFLX 


Swingby continuation analysis option indicator. When 
zero or positive, NSWING is the number of full revolutions 
which the post-swingby heliocentric trajectory must have, 
and subroutine A CONIC (not MINMX3) is used to obtain the 
solution. When negative, NSWING has the same definition 
as MSWTNG (which see), and the MINMX3 iterator is used; 
NSWING must be used when MOPT4(l) > 0, and may be 
used when MOPT4(l) < 0. If MOPT4(l) < 0 and MSWING(i) 
= 0, then MSWING(i) will be set to the value of NSWING. 
Related to MSWING, MOPT4 and T2. 

Iterator independent-variable perturbation- increment 
control. 

= 0 No action, 

1 Allows the iterator to vary a given in- 
dependent-variable perturbation Ax 
whenever a neighboring trajectory is 
detected which has a different number 
of thrust switch points than the associated 
nominal trajectory. Ax is varied until 
the same number of switch points is 
achieved. 

Specifies the unit-number for the trajectory tape. 

Pertains to when MPUNCH 5 - 101. 

Ascending node angle (with respect to vernal equinox) 
of primary-target orbit. Input only when MOPT3 = 11. 
Related to CNI, ECI, SAI, SOI, TPI, EMUODD, and 
R A DODD, [deg] 

Array of five elements, the first three of which may 
be currently used. Ascending node angles of inter- 
mediate-target orbits. Input OMIX(i) only when 
MOPTX(i) = 11. Related to CNEX, ECIX, SAIX, 

SOIX, TPEX, EMUODX, and RADODX. [deg] 

Launch-vehicle -independent (i.e. , no launch vehicle) 
trajectory option in which the value of POWFIX is the 
spacecrafts reference power, [kw] 
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PSIGN 


RADODD 

RADODX 

RAP 

REVS 

RPER 

SAI 

SAIX 

SOI 

SOIX 


Flag defining the sense of the launch hyperbolic excess 
velocity relative to the initial primer vector. A value of 
+1. results in the assignment of the geocentric right 
ascension of the excess velocity equal to that of the initial 
primer vector. A value of -1. causes the geocentric 
right asce7ision of the excess velocity to be 180 degrees 
from that of the initial primer. 

Radius of primary target. Input only when MOPT3 = 11. 
Related to CNI, ECI, OMI, SAI, SOI, TPI, and 
EMUODD. [meters J 

Array of five elements pertaining to the radii of inter- 
mediate targets. These inputs are not used at present. 

Apoapse distance of capture orbit about primary target, 
[planet radii] 

Number of complete revolutions of the ballistic trajectory 
generated when the associated input MOPT is used. Must 
be a positive whole number. 

Periapse distance of capture orbit about primary target, 
[planet radii] 

Semi-major axis of primary-target oi’bit (must be 
positive). Input only when MOPT3 - 11. Related to 
CNI, ECI, OMI, SOI, TPI, EMUODD, and RADODD. 

[AU] 

Array of five elements, the first three of which may be 
currently used. Semi-major axes of intermediate - 
target orbits (must be positive). Input SAIX(i) only 
when MOPTX(l) = 11. Related to CNIX, ECIX, OMIX, 
SOIX, TPDC, EMUODX, and RADODX. [Au] 

Argument of perihelion of primary-target orbit. Input 
only when MOPT3 = 11. Related to CNI, ECI, OMI, SAI, 
TPI, EMUODD, and RADODD. [deg] 

Array of five elements, the first three of which may be 
currently used. Arguments of perihelion of intermediate- 
target orbits. Input SOIX(i) only when MOPTX(i) = 11. 
Related to CNIX, ECIX, OMIX, SAIX, TPIX, EMUODX, 
and RADODX. [deg] 
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SPIRET 

STATE 


STEP1 

STEP2 

TCOAST 


TDV 


TGO 


Retro-stage specific impulse (pertaining to the retro- 
maneuver at the primary target), [sec] 

Array of six elements containing the Cartesian position 
and velocity components of the primary target. Use 
only when MOPT2 = 0 and the trigger settings of Yl(2) 
through Y6(2) are 0 or 1. [AU, AU/tau] 

(tau = 58. 132440991 days) 

Thrust-phase computation step size, Au. 

Coast-phase computation step size, A/?. 

Array of twenty elements, consisting of the durations 
of the coast phases corresponding to the coast-phase 
start-times input in the associated array TOFF. 

[days] 

Time of occurrence of an impulsive deep space burn, 
in days from the start of the trajectory, which may be 
used only if the entire trajectory is ballistic (i.e. , electric 
propulsion is not permitted with this option, nor is a third 
intermediate target). Iterator independent variables 
X64, X65, and X66 must be turned on, as these are used 
as the Av vector components of the deep space burn in 
EMOS. Also, set MAXIIAM « 0. The following special 
feature is available regarding a first intermediate -target. 

If 1.D5 < TDV < 2.D5, then the burn occurs (TDV - 1.D5) 
days after passage of that target; if TDV > 2.D5, the 
burn occurs (TDV - 2.D5) days before passage of that 
target, [days] 

Ballistic trajectorj^-extension print option. When zero, 
no action. When positive, TGO = the number of days 
that the trajectory is to extend ballistically beyond the 
primary-target when no swingby- continuation is re- 
quested, and ballistically beyond the (last) post-swingby 
target when swingby-continuation is requested (in addition 
to the post-swingby trajectory segment itself). Any nega- 
tive value will invoke printout of only the post-swingby 
trajectory segment or segments when swingby-continuation 
is requested. Applies also to trajectories with multiple 
swingbys . [days ] 
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THRET 

TOFF 

TPI 

TP IX 

T POWER 

TSCALE 

T2 

XANG1 

XANG2 

XSWING 


Retro-stage thrust, used only when MTMASS = 2. [lbs J 

Array of twenty elements, consisting of the times, in 
days from the start of the trajectory, at which imposed 
coast phases are to begin. Times must be in ascending 
order. Related to TCOAST. [days] 

Time from reference date (MYEAR, etc.) to perihelion 
passage, for the primary target. Input only when 
MOPT3 = 11. Related to CNI, ECI, OMI, SAI, SOI, 
EMUODD, andRADODD. [days] 

Array of five elements, the first three of which maybe 
currently used. Times from reference date (MYEAR, 
etc.) to perihelion passages, for the intermediate targets . 
Input TPDC(i) only when MOPTX(i) = 11. Related to 
CNIX, ECEX, OMIX, SAIX, SOIX, EMUODX, and 
RADODX. [days] 

Solar-cell degradation characteristic-time; nuclear 
electric propulsion radioactive-decay characteristic - 
time. Related to MPOW. [days] 

Iterator dependent-variable tolerance -intei’val scaling 
factor; scales all tolerances multiplicatively by the 
amount TSCALE. 

Array of ten elements consisting of initial estimates of 
swingby- continuation trajectory-segment flight -times, 
i.e., T2(i) correspond to MOPT4(i). [days] 

Latitude of the launch site. Used only if LAUNCH is 
non-zero. Related to XANG2. [deg] 

Maximum parking orbit inclination, permitted by range 
safety considerations. Used only if LAUNCH is non- 
zero. Related to XANG1. [deg] 

Array of velocity vectors consisting of initial velocity 
guesses of a given post-swingbj' trajectory segment. 

Used only when either NSWING or MSWING has a value 
of -4 or -5. See especially the description of 
MSWING = -4. Velocity consists of exactly the same 
values as found in the VX, V2, V3 locations of the 
(Continued on next page). 
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XSWING trajectory block print (first block) . Related to 

(cont) MS WING, NSWING, MOPT4, and T2. [AU/tau] 

XO Array of seven elements, the first six of which contain 

the Cartesian position and velocity components of the 
launch planet. The seventh element is not used for 
input. Used only when MOPT2 = 0. [AU, AU/tau] 


The following describes the iterator independent and dependent variable 
arrays of the boundary value problem. Input pertaining to the individual inde- 
pendent parameters is contained in the arrays XI through X70. The independent- 
parameter arrays have five elements for each variable, as follows (where i = 

X, 2, 3, .... 70): 

Xi(l) Input value of parameter. Must be input regardless of 

trigger setting. If trigger is on (i.e., Xi(2) =1), input 
value is used as initial guess of independent parameter 
and is varied at each subsequent iteration. If trigger is 
off, the parameter is not used as an independent para- 
meter and is not changed. 

Xi(2) Trigger indicating whether parameter is to be an in- 

dependent parameter in boundary value problem. 

Xi(2) = 0 Not an Independent parameter. 

(Trigger is "off"). 

1 Use as Independent parameter. 

(Trigger is "on 1 '). 

Xi(3) Maximum change to parameter permitted in a single 

iteration. Should be a positive quantity. Used only 
if trigger is on. Units are same as that of the parameter. 

Xi(4) Perturbation increment used to compute partial de- 

rivatives by finite differences. Used only if trigger 
is on. Units are same as that of the parameter. 

Xi(5) Weighting factor. Should be a positive quantity. A 

value of 1. Is generally recommended. The larger 
the weighting factor, the more the parameter is in- 
hibited from varying. Used only if trigger is on. 
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The independent variables are as follows: 


XI 

v> 

X2 

V 2) 

X3 

\< 3 > 


. 

X4 

A o (l) 

X5 

y> 

X6 

A o (3) 

X7 

*,0 

xs 

X 

r 

X9 


X10 

6 


Initial primer vector. 


Initial primer derivative. 


Initial mass -ratio ad joint -variable. 

Propuls ion -time adjoint -variable. 

Not used. 

Geocentric declination of launch hyperbolic 
excess velocity, [deg] 


There is no conversion from input to internal units for any of the adjoint variables 


Xll 
XI 2 

XX 3 


Reference thrust acceleration, g. [m/sec ] 

Electric propulsion system jet exhaust speed, 
c. [m/sec] 

Launch hyperbolic excess speed, v 
[m/sec] 


X14 


Hypei’bolic excess speed at primary target, v 
[m/secl n 


XI 5 


Initial time, t Q , measured from the reference 
date (MYEAR, etc.), [days] 


X16 Time at the primary target, t measured 

from the reference date (MYEAR, etc.), 
[days ] 
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XI 7 


XI 8 
XI 9 
X20 


x 

o 



z 

o 


Launch parking orbit inclination, i. Used 
only if LAUNCH ^ 0. Optimized internally 
by the program if both XI 7 and Y17 triggers 
are off. [deg] 


Initial spacecraft heliocentric velocity. Not 
required unless one of the three triggers is on. 
[AU/tau] (tau = 58.132440991 days) 


X21 


Constant thrust cone-angle, 0. Non- zero value 
invokes the constant-0 constraint. 0 < 0 ^ 180°. 
Zero -value implies that 0 is optimized along the 
trajectory (variable 0). [deg] 


X22 through X29 are currently not used (although some locations following X21 
are reserved for additional constant thrust cone-angles). 


X30 


X 

s 


Degradation -t ime adj oint-var table . 


X31 through X40 are currently not used. X41 through X50 pertain to the first 
intermediate target, X51 through X60 pertain to the second intermediate target, 
and X61 through X70 pertain to the third intermediate target. The correspond- 
ing intermediate- target parameters are ignored if the intermediate target is 
absent. Subscripts 1, 2, and 3 pertain to the first, second, and third inter- 
mediate targets, respectively. 


X41 

V) 

X42 

(2) 

X43 

A x (3) 

X44 

Ajfl) 

X45 

\< 2 ) 

X46 



* 

Primer vector (at start of trajectory segment) 


l Primer derivative (at start of trajectory segment) 
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X47 

X48 

X49 

X50 


Encounter speed at first intermediate target, 

v . . [m/sec] 

05 1 

Time at the first intermediate target, t^ , 
measured from the reference date (MYEAR, 
etc,). [days] 

Sample -mass factor, k „ . for sample - 

samp 1 

retrieval at first intermediate target. 

Drop-mass factor, k. „ , for instrument- 

drop 1 

package dropoff at first intermediate target. 


The independent variables X51 through X60 and X61 through X70 are 
identical to X41 through X50 except that they pertain to the second and third 
intermediate targets, respectively. A third intermediate target may not be 
present when simulating ballistic missions having a deep space burn (See TDV), 
in which case X64, X65, and X6S are used as follows: 

X64 Ax ’ 

X65 Ay V Deep-space velocity-increment. [AU/tau] 

X66 Az 

Inputs pertaining to the individual dependent parameters are contained in 
the arrays Y1 through Y70. The dependent -parameter arrays have three ele- 
ments for each variable, as follows (where i-1, 2, 3, ..., 70): 


Yi(l) Desired value of the dependent parameter. 

Yi(2) Trigger. If off (i.e., equal to zero), the para- 

meter is ignored and is not considered a dependent 
parameter. Then the other two inputs pertaining 
to the parameter need not be input. If trigger is 
on (i.e., not equal to zero), the parameter is 
considered to be a dependent parameter or con- 
straint. Certain of the parameters may have 
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Yi(2) up to three non-zero trigger settings. These 

(cont) will be discussed individually below. 

Yi(3) Tolerance of desired value (full interval width). 

It should be noted that the transversality conditions, which comprise some of the 
parameters, are developed under the assumption that all constraints are of the 
point constraint type. Therefore, the satisfy-mode is sufficient in solving any 
optimization problems for which a complete set of transversality conditions is 
available. 

The dependent-parameter arrays are as given below. T(x) represents 
"the transversality condition associated with x" and the function T(x) will have 
different values depending upon the constraints imposed on the problem. See 
NOMENCLATURE of reference [l] for definition of symbols and subscripts. 


Trigger 1 


Trigger 2 



Y1 

Ax [AU] 
n 

T 

f 

a [AU] 

Solar distance* [AU] 

T(tJ) 

Y2 

Ay [AU] 
n 

f 

| 

e 

T(0j* 

t 


T ( e t ) 

Y3 

Az CauI 
n 

1 

f 

i [deg] 



T <V 

Y4 

Ax [AU/tau] 
n 

1 

t 

T(ft) 

T (x n ) ■ 


T(0 

Y5 

Ay [AU/tau] 
n 

1 

t 

T (60) 

T V 

optimal 

flyby 

v 

o=0 

Y6 

Az [AU/tau] 
n 

T 

T(f) 

T(z ) 
n 


T(A) 


* Applicable only for two-dimensional motion in the xy plane. Also requires that 
MOPT2 - 0. 


Under Trigger 1 above, the first set of conditions applies to ordinary targeting 
conditions for position and velocity, and also to extra-ecliptic conditions to be 
satisfied when TOUT = 1; the second set of conditions applies to extra-ecliptic 
missions when TOUT = 2. T{60), and T(f) are symbols for the transversality 
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NERVA 



conditions yielding optimum final node angle, argument of perihelion, and true 
anomaly, respectively. 



Trigger 1 

Trigger 2 

Trigger 3 

Y7 

V 

n 

X 

Vn 

m , [kg] 
net 

Y8 

T(T) 

T [days ] 


Y9 

Currently not used. 


Y10 

T(6) 

<5 [deg] 

Used only if LAUNCH ^ 0. 

Yll 

T(g) 

g [m/sec^] 

P f [kw] 
ref 

Y12 

T(c) 

c [m/sec] 


Y13 

T(v ) 

v [m/sec] 
coo 


Y14 

T(v ) 
x ®n' 

v [m/sec] 
®n 

extra-ecliptic inclination [deg] 

Y15 

T(t ) 
o 

t [days] 


Y16 

T(t ) 
n 

t [days] 
n 

t - t [days ] * 
n o 


*Time transversality with flight time fixed is assigned to Y15 under trigger 1, 

Y17 T(i) i [deg] , Used only if LAUNCH £ 0. 

where i = parking 
orbit inclination. 


Y18 

T (x o ) 

x [AU/tau] 
0 

Y19 

T (y 0 ) 

y [AU/tau] 
0 

Y20 

T(z ) 
0 

z [AU/tau] 
0 

Y21 

T (0 ) 

0 [deg] for 0 = 
constant with time. 

Y22 through 

Y29 are 

currently not used. 
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Trigger 1 Trigger 2 

Y30 T(s) s [days] (Degradation time) 

Y31 through Y40 are currently not used, Y41 through Y50 pertain 
to the first intermediate target. 


Y41 

Ax, [au] 


Y42 

Ay 1 [AU] 


Y43 

Az [au] 


Y44 

Ax x [AU/tau] r T(x ) ' 


Y45 

[AU/tau] r T(y 1 ) 

optimal 
> flyby 

Y46 

Az x [AU/tau] | T(z ) 


Y47 

v „ [m/sec] 

00 1 

Y48 

T (t^) t i [days 

.] 

Y49 

m samp 1 


Y50 

m dropl Ckg] 



Y51 through Y60 and Y61 through Y70 are identical to Y41 through 
Y50 except that they pertain to the second and third intermediate targets, 
respectively. 
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G. 


DEFAULT VALUES OF INPUT PARAMETERS 


The following is a complete, alphabetical list of the default values of 
program input quantities having non-zero default values, except for the iterator 
arrays o AH other inputs are zeroed. The default values of the iterator arrays 
Xi(l), Xi(2), Yi(l), and Yi(2), for i = 1, 2, 3, . 70, are zero, and the default 
values of Xi(3) through Xi(5) and Yt(3) for the same range of i are listed in 
the listing of program inputs of Sample Case C. 


ALPHAA 

15. 

ALPHAT 

15. 

AN 

1.5 

AR 

1. 

BI 

.76 

CTANK 

.03 

CTRET 

1/9 

DI 

13. 

GAMMAX 

1. 

GAP 

.0001 

HOUR 

12. 

IRK 

1 

IRL 

1 

ITF 

3 

MAXHAM 

5 

MDAY 

1 

MODE 

4 

MONTH 

1 

MOPT3 

10 

MUPDAT 

1 

MYEAR 

1975 

NDIST 

3 

NHUNG 

25 

NPRINT 

7 


NSET(3) 

300 

NSET(5) 

300 

NSW PAR 

1 

NTAPE 

17 

POWFIX 

-1. 

PSIGN 

1. 

RADODD 

1. 

RAP 

38. 

RPER 

2. 

SAI 

1. 

SPERET 

300. 

STATE(l) 

1. 

STATE(5) 

1. 

STEP1 

.03125 

STEP2 

.125 

TDV 

-1. 

TGO 

-1. 

THRET 

400. 

TOFF 

20* -I. 

T POWER 

10. **30 

TSCALE 

1. 

T2(i) 

50*i 

X0(1) 

1. 

X0{5) 

1.0000015 
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rv. 


SAMPLE PROBLEMS AND RESULTS 


Three sample cases have been selected for presentation in the following 
page s 0 Each case was selected to display the use of one or more new features 
that have been added to the program. The first case is an extra-ecliptic mission 
that exhibits the effects of high launch asymptote declinations on launch vehicle per- 
formance and the low thrust trajectory. This case also displays the housekeeping 
power option that was recently added to the program. The second case is a comet 
rendezvous mission which includes the effects of solar array degradation due to 
radiation effects. The third and final case displays the HILTOP’s powerful capa- 
bility for ballistic mission design and optimization. The specific mission chosen 
is a cometary flyby past Giaeobini-Z inner followed by a deep space burn 10 days 
after passage, a return to and swingby of Earth (unpowered), a second swingby of 
Earth (powered), and finally encountering the comet Boi-relly nearly 1023 days 
from launch* The tremendous flexibility for creating imaginative, multi-target 
mission profiles is demonstrated in this example. 


preceding page 


blank not filmed 
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A. 


EXTRA ECLIPTIC MISSION 


The objective of this mission is to place maximum payload in a circular 
orbit of 1.001 AU radius inclined 45 degrees to the ecliptic. The mission is 
2-J years (912.5 days) in duration and departs from Earth parking orbit on 
April 21, 1979. The solution given in this example is contained in the class of 
six-burn solutions which, for the mission duration assumed, tends to restrain 
the trajectory from deviating far from the nominal 1 AU solar distance. The 
specific case chosen uses the Titan III E/Centaur launch vehicle. The reference 
power (the power delivered to the power conditioners at 1 AU) is 20.35 kw 
and the specific impulse of the thrusters is 3000 seconds. The launch excess 
speed is optimized. The extra-ecliptic end conditions are invoked by setting 
IOUT = 1 and defining the desired values for AR, AE and AAI, the radius, 
eccentricity and inclination, respective^. 

This case exhibits the use of several optional features of the program. 

A total of 0.65 kw of power developed by the solar arrays is reserved for 
housekeeping (non-propul sive) uses. This option is triggered with the input 
DPOW which is the ratio of housekeeping power to reference power. The power 
delivered to the power conditioners at distances below 1 AU is not permitted 
to exceed that delivered at 1 AU. This constraint is invoked by setting MODE 
equal to 5 and GAMMAX (the maximum permissible value of the power factor 
y) equal to 1. The effects of launch asymptote declination are included in the 
launch vehicle performance model by setting LAUNCH equal to 1. The equatorial 
inclination of the launch parking orbit is limited to a maximum of 36 degrees 
through the input parameter XANG2. Since the geocentric declination of the 
launch asymptote for extra ecliptic missions is usually much greater than this 
inclination limit, the solution will include a non-coplanar injection maneuver from 
the launch parking orbit. The declination of the launch asymptote is optimized. 
Finally, the option of inputting the coefficients of the power profile is illustrated. 
The inputs for this case are listed on the next page. 
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It should be ruled that the choice of final orbit radius of loOOl AU rather 
than l o 0 AU was made to alleviate numerical difficulties arising as a result of 
the corner in the power curve at 1 AU„ Neighboring trajectories terminating 
on opposite sides of the corner point tend to possess different partial derivatives 
(i 0 e» , they will behave differently when subjected to the same perturbation). Con- 
sequently, if the final desired distance is exactly the point of the discontinuity, 
one might expect convergence retardation when the end conditions are nearly 
satisfied* 


m INPUT X1{2)=1.D0,X2(2)=1.D0,X3{2)=1.D0,X4(2)=1.DG,X5(2)=1.D0 
XG( 2>*l.D0,X7*l.D0,X10(2)«l.D0,Xll(2)«l.D0,X12*2.941995O4 
X13( 2) *1. DO, XI 5 = 3. 10 3 7001, X 16 =9 .4 35370 2 

Yl(2)=l.DO,Y2(2)“l.DO,Y3(2)=l.nO,Y4(2)=l.DO,Y5( 2)=1.00, Y6(2)=1.00 
Y10(2)=1.DO,Y11=20.35DO,3.00,Y13(2>=1.00 
LAUNCH = 1,M BOOST* 15, M TM ASS = 3,M00 E = 5 
.MO PT2 = 3,MO PT3 =0,MY EAR=19 79, MO NTH = 3 , MO AY = 21 

BI=.6300,DI=0.00, CTANK*. 03500, 0AMHAX=1 . 0 0,0 P0W*3 . 1041031941030-2 
1 OU T = l, AAI =4 5.D0,AR = 1. 0 01D0, XANG1 = 28. 50 0, XANG2=36 .DO 
AS0L = 1. 4 38 2D0, 0 . D 0, - . 2 2 3 50 0, 0 . D 0 , - . 2 14 70 0 
Xl*-1. 8492690168360-02, X2*-3 . 444545S696 870-01, X3— 5. 29 23946 80826D 00 
X4 = 7.2838600516620-01, X5= 0.4 8995 82 7726 30-01, X6--3. 27764 1 8049 73D-01 
X10*-3. 8091507295560 01,X11= 2 . 79 05 74 06 12 560-04, X13= 5 . 32 86 06 5 3464 5D 03 
&END 


*^< 83 ? 
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05 


Y55 *» 0.0 « °*° 

Y56 - 0.0 » 0-0 

Y57 * 0-0 • 0.0 

Ysa « o.o * o*° 

Y59 = 0.0 . 0 .0 

Y60 = 0.0 ► 0.0 

Y61 = 0 .0 .0.0 

Y62 = 0.0 » 0.0 

Y63 = 0.0 .0.0 

Y64 a 0.0 ,0.0 

Y 6 5 =0.0 .0.0 

Y 66 =0.0 « 0.0 

Y 67 = 0.0 ,0.0 

Y 69 = 0.0 * °*° 

Y 69 =0.0 • 

Y 70 =0.0 ♦ 0»° 

tOUT PUT 

A A I = 4S. 0000000 0 OOO 0 00 

0.0 » an= 

-0.2235000000000000 , 0.0 

0.0 » 0*0 

0.0 . 0.0 » 

ctret® o.iinimiumu ,di= o.o 
O.o , ECI X= 0 .0 

0.0 ,£1=0.0 

0.0 ,0.0 * 

0 . 9 9999 99999999999D- 04 , HCUR= 12 .00000000000000 




Q? 


,AE= 0.0 

1 . 500000000000000 


9. 99? 999999999999 D-05 , 

9. 99999999999 9 99 9 D-0 5 . 

9. 999 999999999999 D-0 5 » 

9. 999 9 9 99 999 9 9 999 D-0 5 , 

9. 999 9999999999990-05, 

9. 999 9 999999999 9 90-05 , 

9. 999 9 999999999990-05 » 

9. 999 9 999999999 990-0 5, 

9.999 999999 99999 90-05. \J,£ 

9. 999999999999999D— 05, 

9. 9999999999999990-05 , 

9. 999 9 999999999990-0 5 . 

9. 9999999999999990-05 , 

9. 99? 9 99999999999 D-0 5 , 

9 . 999 9999999999 99 D — 0 5 • 

9.99? 99999999 99 990-0 5 , 

^AU^iAAs 15.00000000000000 , AL S HAT= 15.00000000000000 

,AR= J.001300000030000 . ASOL= 1. 638 2C 00 000 00 000 . 0.0 




, ACT STU = 


0.0 

1 .000 0000 000 0 000 0 , R AOOD 0= 

, 0.0 

2 . OOOO 0000 0000000 « S A 1 = 1.00000000000 00 00 


, 0.0 
, 0 . 0 

1 . COO 000 0 000 0 00 0 0 , 

,STEP1= 
, 0.0 

. 0.0 
. 0.0 
, 0.0 
,THREt= 


0.0 

PSIGN= 

0.0 

RP£R = 

0.0 
0.0 

5TATE= 

0.0 

0.0 

0.0 

0.0 

0.0 

-1 .000000000000000 

-1 .0 000 00 COO OOO 0 00 

-1 . 000000 oocoooooo 
-1 .000000000000000 
-1 .ooooooooooooooo 
0.0 

1 .ooooooooooooooo 

2 50 .0000 coo 0000 00 

5 0 0. 0000 CO 0 0000 00 
0.0 
0.0 
0 . 0 
0.0 
0.0 
0.0 
0.0 
0.0 


,-0.2 1 A 70000 OOO 0 0000 
, CNI = 0.3 
, 0.0 

*O^REr R= 0.0 

, .0 .0 

,EMUOOD= 0.3 

, 0. 0 

, 0 M I = 0 .0 

, 0.0 

OOOOOOOOOOOOOOO ,T AOODX= 

, 0.0 


,ei = o.o 


,B2 = 


0.0 


0.0 


0 .0 


0 CAP = 


0.0 

, RAP® 


, 0.0 

, 0. 0 

0.0 

0.31 250 00000 0 000 000-01 

, 0.0 
, 0.0 
. 0.0 
, 0.0 

400. OOO 00 00 00 000 0 

, -1 .ooooooooooooooo 
, -1 .ooooooooooooooo 
, -1 .ooooooooooooooo 
, -1 .ooooooooooooooo 

, 0.0 

,T2= 50 .00000000000000 

300.0 0000 00 00 00 00 


,S4lX= 0.0 

,SOI= 0.0 

. 0 .0 

,o.o ,o,o 

, STEP 2= 0.12500000000000 00 

. 0.0 
, 0.0 
, 0.0 
, 0.0 

,TOFr= -t .000003300000000 
-1 .00000 03 00000000 

-1 .ooooooooooooooo 
-1 .00000 0330300000 
TP 1 = 0.0 
0.0 

, 100. OOO 00 OOO 00 090 

350.0000300300000 


,SI= 0 .6300000000000000 

, CN I X = 0 .0 - 

C5TR= 0.0 ,CTANK= 0 , 34999999999999990-01 

, D?Qb'= 0,3 1941031 94 10299 90-01 «£C J «* 

,0.0 ,0.0 • 

,EMJODX= 0,0 • 

, GAMMAX* 1 .ooooooooooooooo 

,OMIX= 0.0 » 

,P0WF1X= -1.000000000000000 
, 0.0 • 
30.00000000000000 „REVS= 0.0 

, 0.0 ♦ 

„S01X= 0.0 

» SPI R£T = 300.00 0 00 00000 000 

, 1,00 00 000000 OO OOO 

,rcOAST= 0.0 « 

, 0.0 • 

, 0.0 • 

, 0.0 • 

,TDV= -1 .ooooooooooooooo 

-1.0030 OOO OOO 00 OOO d -1.0 00000000000000 

-1 .ooooooooooooooo • -l .ooooooooooooooo , 

, - 1 ,0 0 00000 0300000 0 , -1 .cooocoooooooooo , 

, -1*000000000003000 * —1.000000000003000 , 

,TPIX= 0.0 • 

,TPOWER= 0. 1 00000 000000000 OD 31.TSCALE® 


,TGO= 


,XANG1« 20.50000000000000 ,XAnG2= 


MPR 1NT= 


, 0.0 
, 0.0 
. 0.0 
, 0. 0 

♦ 0.0 
, 0.0 
, 0 . 0 

, t BAL = 

0. 1 TF = 3,ITPRNT= 

1 , LOAOX= 0 . MAX HAM® 

O.MOPTX® 0, 

0, 0 , 0 , 

0, MPUNCH= 0 ,MREAD= 

0. O, O. 

25. NORMAL= 


36.00000000000000 

. 0.0 


1 50 .00000000 00000 
400.0000000003030 , 

, XSW ING= 0. 0 


200.0000000000000 
450.0000000000000 


3.NHUNG* 

0. 


0.0 
. 0.0 
. 0.0 
, 0.0 
, 0.0 
. . 0.0 
. 0.0 

0,INTPR= 0,IOJT= 

0. JPP= 3 ♦ JPRINTs 

5.M800ST® 15.MDAY® 

0, 0 , 

0, 0, 

0 .MS H I N3 = 0 

3 , Mt MASS* 

N>£RF* O.NPRINT= 


0 , 


300 . NSW 1NG= 


0 . NSWPAR= 


1 . NT APE = 


0.0 
0.0 
0.0 
0 .0 
0 .0 

1 .000001500000000 
1,1 RK= 1 

0, JT = 

21 ,MOl>E = 

0 , MQPT2= 

0, 

0, 

3 , HUPOAT = 

7, NSET= 
1 7 


, 0.0 

• 0.0 

, 0.0 

& 0.0 

. 0.0 

, x 0= 1 .003 0000000000 00 • 

,0.0 a 

,HLs 1,1 ROT = 0.1SPIN= 

0,KPARt= o, launch® 

5 , MONTH® 3»M0PT= 

3, M0PT3= o,MQPTA= O 

0, 0. MPOW= 0, . 

0. o, o, 

1 .MYEARo 19 79.NO I ST = 

0 , 0 , 300 . 


tEND 
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NOMINAL. TRAJECTORY 1 (TOTAL 


INHIBI TOR IS 


5.82080-1 1 


. 1 .PRIMl (-1 .64926900-02 J 2.P 
6. PDOT 3 <-3.277 641 80-0 1 ) I 0.0 

l, DELTA X< 1.367340-02) 2.0 

6.DELT ZD{ 2. 579040-02) 10. T 

THRUST SWITCHING TIMES ( DAYS ) 
439.175 ON 567.999 OFF 


POWER 

20.3499999990 


IN IT IAL 

3123.1976935923 


INDEPENDENT PARAMETERS 

2.PRI M2 ( — 3. 44 454 5 90— 01) 3 .PR I M 3( -5 . 29239470 00) 

10. DECLN(-3. 80915070 01) 11 .ACCEL< 2.790 57 4 10-04) 

DEPENDENT PARAMETERS 

2.0ELTA Y< -2. 04257D-02 ) 3, DELTA Z( 2.52203D-02) 

10,T.DECLN( 1.75024D-09) U, POWER ( 2.03500D 01) 

S) 0.0 ON 80.374 a— 103.276 ON 

OFF 618.895 ON 749.493 OF = 809.322 ON 

ELECTRIC PROPU.5ION PARAMETERS 
EFFICIENCY PROP TIME J 

0. 6300000000 71 2.1252770929 1 1 .3661826370 

MASS COMPONENT BREAKDOWN 
R0PUL5 ION PROPELLANT TANKAGE 

0.2499999709 1813. 392062 9241 63.4683722023 


4. PDOT 1 C 7 .28386 0 ID-01 ) 
13.VINFM 5 . 3286 06 50 03 1 

4, CELT XD< 1.623970-02) 

13. T, VI NFM -3. 7X29 60 -091 

218.909 OFF 247.259 
912.530 ON 


PROPULS ION 
620.2499999709 


PROP TIME HAT 10 
0.78041 :;2626 


STRUCTUR E 

O . 0 


S»PDOT2( 6.4 09 95 8 3D— 01 » 


S.DELT YOI-2.40B070-02) 


AVE ACCEL 
0.00043091 16 


PAYLOAO 

626,0972534950 


— NOMINAL TRAJECTORY 2 (TOTAL 4) 


INHIBI TOR IS 1 .8190D-12 


1 .PRIMl (-9. 005 854 20-03 I 2.P 
6 .P00T3<-3.3398207D— 01 ) 10.0 

l.OELTA X(- 1. 06725D-05 ) 2.D 

6.DELT ZD< 2.135280-031 10.T 

THRUST SWITCHING TIMES (DAYS) 

439.396 ON 566,714 OFF 


POWER 

20 . 3 066 S7 325I 


I N I T 1 AL 

3262.5709661616 


INDEPENDENT PARAMETERS 

2.PRI M2 (-3.51759130-01 > 3 .PR I M31 -5 . 2 00 B632D 001 

10.0ECLN(-3 .7843 9590 01) 1 1 ,ACCEL{ 2.6655743D-04) 

DEPENDENT PARAMETERS 

2. DELTA Y( -9.998500-04 > 3, DELTA Z ( - 2 . 5 1 1 l I 0-03 I 

lO.T.DECLNt 2.709670-03) 11, POWER ( 2.030670 01) 

S) 0.0 ON 80. 395 0= ~ 103.651 ON 

OFF 616.029 ON 747.829 OF- 834.878 ON 

ELECTRIC PROPULSION PARAMETERS 
EFFI C CENC Y PROP T I ME J 

0.6300000000 715.6762058950 9.9909155676 

MASS COMPONENT BREAKDOWN 
ROPULSION PROPELLANT TANKAGE 

8.9289536062 1819.7281211922 63.5904842417 


PROPULSION 
61 8.9289536062 


4 , PDOT 1 ( 6.991 09700-01 I 
13.VINFM 5.1101. 37 7D 03) 

4,DElT XO( 1 .1 3.9830-03) 

13. T.VINFU -1.6 90610-0 3) 

219.58* OFF 247.918 

912.500 ON 

PROP TIME RAT 10 
0.7843026914 


STRUCTURE 

0.0 


5 , POOT2 ( 6.28887160-01 ) 


S.DELT YD ( — 7 « 224 38D— 0 4 ) 


AVE ACCEL 
0.0004008460 


PAYLOAD 

760.2234071214 


NOMINAL TRAJECTORY 3 (TOTAL 7) 


INHIBITOR IS 5.68430-14 


INDEPENDENT PARAMETERS 

1 .PRIMl ( 2.2736 0150-0 3 ) 2.PRIM2(-3.74 36282D-01I 3,PRIM3(-5.192 236 7D 0 0) 

6, PDOT 3 (-3.22266600-01 ) 1 0 . OECLN (-3 . 785484 OD 01) 11 ,ACCEL( 2.66600440-04) 

DEPENDENT PARAMETERS 

1, DELTA X( 9.717230-04) 2, DELTA Y ( -4 . 332260- 04 ) 3 , DELTA Z ( - 1 . 0 1 0 0 1 0-0 3) 

6.DELT ZO( 3.845110-04) 10, T ,OECLN< - 3 . T4 S 2 HD - 05 ) II . POWER ( 2.035000 01) 


THRUST SWITCHING TIMES (DAYS) 
439.925 ON 568.200 OFF 


0.0 

617.759 


POWER 

20.3499746867 


I N I T IAL 

3269.1 258074603 


EFFI CIENCY 
0.6300000000 


PROPULSION 

620.2492284722 


ON 8 0.29 9 OF r 10 3.546 ON 

ON 749.422 OF- 806.781 ON 

ELECTRIC PROPU-SION PARAMETERS 
PROP TIME J 

715. 011350 6695 9.37048199 92 

MASS COMPONENT BREAKDOWN 
PROPELLANT TANKAGE 

102 1 .61 64 01 3205 63 . 75657684 62 


4, PDOT 1 I 7 .09035730-01 ) 
13.VINFK 5.09859030 03) 

4.0ELT XDI-5. 133860-04) 
13.T.VINFK 3.24 076D-05I 

220 .043 OFF 248.369 

912.530 ON 

PROP TIME kAT 10 
! 0.7835740829 


5 » POOT 2 ( 6.36032500-01 ) 


STRUCTURE. 

0.0 


5.0ELT YD( 4.326990-05) 


AVE ACCEL 
0.0004006508 


PAYLOAD 

763.503520821* 


THIS CASE IS CONVERGED, 


8 TRAJECTORIES WITHOUT PARTIAL 0ERIVATIVE5 ANO 3 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FOR THIS CASE, 



case i 


SWITCH POINT SUMMARY 


PAGE 1 


TIME 

SEMI-MAJOR AXIS 

ECCENTRICITY 

INCLINATION 

NOOE 

ARG POS 

RMAG 

travel 

ftl 

R2 

R3 

VI 

V2 

V3 

MASS RATIO 

thrust acc 

LI 

L2 

L3 

L4 

L5 

L6 

L7 

MAM 

LG 

LC 

LPHl 

CONE 

CLOCK 

HMAG 

POWER FNCT 

SWITCH FNCT 

PSI 

THETA 

PHI 

LATITUOE 

LONGITUDE 

FLT PTH ANG.E 

VMAG 

PROP TIME 

EARTH 



START OF TRAJECTORY. THRUST ON 




o .o 

8.989777750-01 

1.211808740-01 

9.228427590 30 

3.093485450 01 

1 .BOOOOOOOD 02 

1 .004982350 00 

0.0 

—8 . 629263 290—0 1 

-5.151190790-01 

0.0 

4.964306953-01 

-7.601992670-01 

-1 .50185351 D-01 

1 .OOOOOOOOD 00 

A. 4654*9540-02 

1 .433022090-03 

—3 .7 127347 0D-0 1 

—5 .1922233 2D 00 

7 .085591 *63-31 

6.35 2863 4 8D- 0 1 

-3.264391610-01 

1 .OOOOOOOOD 00 

9.491 114440-02 

0.0 

0.0 

0 .0 

8.605271660 31 

1.71 366591 D 01 

9.41 156981 D-01 

9.930951630-01 

4.193075260 00 

-7.708403020 01 

8.064814240 01 

8.791846720 01 

0.0 

— 1 .4916514 6D 02 

-1.612086370 00 

9.368618580-01 

0.0 


SWITCH THRUST OFF 


8.03383162D 01 
3.411006150-01 
3 .977283 74D-0 1 
—3*874708 670 00 
-6.041975990 01 


9.240153250-01 
-8.023345320-01 
-4.694863440-02 
-7. 3016373 30-02 
4.478400230 01 


1 .0302221 70-01 
-1 .761244690-01 
-1 .119306190 00 
0.0 

6.94 6 999 1 0 0 01 


1 .196780310 01 
9.423323710-01 
—4 .701 773373-02 
9.740 7950 7D 01 
-l .142093890 01 


4.065933450 01 
5.256348160-01 
-6.320616820-01 
1.611051620 01 
-6.696794930 01 


2 . 5273204 8D 02 
—4 . 56244954 D- 02 
5. 36 71 1 3 51 D 00 
9.561423 62D- 0| 
-5.51 2069660 00 


6.894437360-01 
9.371125330-01 
1 .253030110 00 
1 .000000000 00 
1.079983000 00 


8. 236967300 01 
4.798246700-02 
9.491115210-02 
-3.996802890-15 
8.033831620 01 


a 


l| 

If 

O Q+ 

w4 

p2 ^ 

Q O 


SWITCH THRUST ON 


1.0 35874 76D 02 
6.63125730D-01 
4 .0732751 3D— 01 
— 3 . 874798 8 7D 00 
7.321791 1 ID 01 


9.240153250-01 
-5.070083970-01 
— 3.99412178 D- 01 
-7.301637330-02 
—3.521 54 79 5D 00 


1 .0302221 7D-01 
-1 .731113720-01 
1 .042972960 00 
0 .0 

7.3250S352D 01 


1.196780310 01 
6.287177063-01 
1.311 270363-01 
9.047700793 01 
-l .171610550 01 


4 • 0 5593345D 01 
9.29B30855D-01 
— 1 .050 1 2549D 00 
1.74 B343240 02 
-3 • 7400 5722D 01 


2.3166T783D 02 
6.258344780- 02 
5.227605910 00 
9.S6142362D-01 
-3.910803900 00 


8.52502666D-01 
9.3711 25330-01 
1 .25303011D 00 
1 .000000000 00 
1 .1 24188860 00 


1.11 324607D 02 
4.798246700-02 
9.491 1 1 521D— 02 
-1.33226763D— IS 
8.033831620 01 


SWITCH THRUST OFF 


2. 200471 31D 02 
-3. 173520540-01 
-2.433325720-01 
-B. 084407140 00 
6.398730570 01 


9.260122330-01 
8. 40718 02 1D-0 1 
-6.420144370-01 
-2.048093540-01 
1.084350530 02 


9.308915210-02 
2 » 376 265 S 2 0—0 1 
1 .197578380 00 
0.0 

9.797198290 01 


1 .575481 410 01 
-9.845399570-01 
-9.530832783-01 
9.927457730 01 
1.481196260 01 


4.107520600 01 
-2.997231 950-01 
4.973095690-01 
1 .6209 7586D 02 
1.106803Q8D 02 


7.031145360 01 
1, 199687430-01 
— 4 .746233 36 D 00 
9.5 85 30572 D- 0 1 
5.33 8755S5D 00 


9.295082010-01 
8.45914448D-01 
1 .611863020 00 
1 .OOOOOOOOD 00 
1 .0 3571 63 1 D 00 


2.60280593D 02 
S. 315545950-02 
9.491 116750-02 
-1 .754 152380-14 
1 . 96797972D 02 


SWITCH THRUST ON 


2 .4841 4939D 02 
-7 .321702 29D— 01 
— 5.74 16241 7D —0 1 
-8.084407140 00 
-5.849641790 01 


9.268122330-01 
5. 8 547 02 B 50 — 0 1 
— 3*784 3 384 2D— 0 1 
-2.048093540-01 
9.923817760 01 


9 .30 8915210-02 
2 . 6023S064D-0 1 
-l . 1 969597 4D 00 
0 .0 

9.481219310 01 


1 .575481*13 31 
-6.831248373-01 
— 4 .74 0 87 03 90- € 1 
8.594777300 31 
1 .551 429030 31 


4.10 75 20600 01 
—7.1 3729664D-01 
S.8574I621D-01 
3.685644790 01 
1.413528820 02 


9.989935080 01 
-2.51642*050-02 
-4.833448640 00 
9.5B53057 2D- 01 
4.51038882D 00 


9.7291 8798D-0 l 
8. 459144460-01 
1 .6U88302D 00 
1 .000000000 00 
9.882827610-01 


2.8 9 66 84 9 0 D 02 
5.315545950-02 
9.491116750-02 
2.886579860-15 
1 .96T97972D 02 



PACE 2 


TIME 

semi-major axis 

ftl 

R2 

LI 

L2 

LG 

LC 

PS1 

THETA 


ECCENTRIC TTY INCLINATION 
R3 VI 

L 3 L4 

LPHt CONE 

PHI LATITUOE 


NODE ARG POS 

V2 V 3 

LS L6 

CLOCK hmaS 

LONGITUDE FLT PTH ANGLE 


RMAG TRAVEL 

MASS RATIO THRUST ACC 

L7 HAN 

POWER FnCT SWITCH FNCT 

VMAG PROP TIME 


SWITCH THRUST OFF 


4.010070630 02 
2 .400741 36D— 01 
2 .355275990-01 
—I .559I2236D 01 
-6.605315360 01 


9. 926792430-01 
-8.601 594 050-0 1 
2. I 35501880-01 
-4.3869621 60— 0 1 
4.9030011 OQ 01 


4.021 5566 6D— 0 2 
-3 . 40779965D— 0 I 
-1 .692560920 00 
0.0 

7 .456609 1 7D 01 


2.212696930 01 
9.467544783-01 
4 .52371 5290-01 
9.627615370 01 
-2.043844010 01 


3.931884210 01 
3.53787121D-01 
-1.130723800 00 
9.56640768D 00 
— 7.42594 3 6 6D 01 


2.47986770D 02 
-1 .31 0635600-01 
4 .56849831 D 00 
9.956271 730-01 
-2.O65043OOD 00 


9.758854130-01 
7.276817250-01 
2 .33764 1 39D 00 
1 .COOOOOOOD 00 
1.020906030 00 


4.361040980 02 
6.179206530-02 
9.491 117440-02 
-2.220446050-16 
3.49390095D 02 


SWITCH THRUST ON 


4. 399952 2* D 02 
7.724348930-01 
6.064626 C3D —0 1 
-1 .559122360 01 
7.682147830 01 


9.9287924 3D- 0 1 
-4.50032489D-01 
-7.88401 03 70-0 1 
-4.3069621 60-01 
-2.34991726D 01 


4.021556660-02 
-3 .405740610-01 
1 .405868520 00 
0 .0 

7 .793 16696D 01 


2.212696933 01 
5.4 72 87 0 2 33-01 
5.70391 350D-01 
1 .0061791 30 02 
-2.085521473 01 


3.931884210 01 
8.751 594840-01 
-1 .60-09807 9D 00 
1.62708336D 02 
-3.022573793 01 


2.890614300 02 
1.342932270-01 
4.292903070 00 
9.956271 730-01 
-9.691 429060-01 


9 .5664 81 04D-0 1 
7.27681 7250-01 
2.337641390 00 
1 .0 00 0 00 OOD 00 
1 .040894360 00 


4.771787580 02 
6.1 79 20 05 3D- 02 
9.49111 7440-02 
-1 .421085470-14 
3.493900950 02 


SWITCH THRUST OFF 

00 5.681922520 02 9.8 805095 30-0 1 2 .60 924695D-02 2.853850543 01 3.92753538D 01 6.34800766D 01 9 .9 1 545 5 74D-01 6..11386147D 02 

-1 .50 680 7 64 D— 0 1 8 . 8 36 26 36 2 D- 0 1 4 .23 865999D-0 1 -9.469304013-01 -2,402550 930-01 2.248546780-01 6 .27292 2 52D-0 1 7. 168105630-02 

2 .74751 786D-01 -8. 1 264 707 70-0 l 1.796782390 00 -1.2080237DD 30 1.5392I638D 00 -3.874477950 00 3.135158230 00 9.491117390-02 

-2.149807820 01 -7.147457240-01 0.0 7.B8119551D 01 1.326645970 02 9.936690950-01 1. COOOOOOOD 00 -1 . 9984 0 1 4 4D- 1 5 

8.021703660 01 8.963810900 01 8 .99385092D 01 2.530774630 01 9.96773020D 01 1.4013643QD 00 1.002476670 00 4.775071220 02 


SWITCH THRUST ON 


6.177877340 02 
-8.126352050-01 
—7 .596 51 4 9 2D— 0 1 
-2.149807820 01 
-B. 496792010 01 


9.880 509530-01 
3 . 9844009 1 D— 0 1 
6. 6875964 0D-01 
-7.147457240-01 
4.8656 843 4 D 01 


2.609246950-02 
4.475007360-01 
-1 . 71 464309D 00 
0 .0 

8 . 66 783 8360 01 


2.853650643 01 
—5.1 3813741D-01 
-1 . 1 I 335732D 00 
7 .332521 31 D D1 
2.630982253 01 


3.927535380 01 
-8.222761 69D-01 
1.737417660 00 
3.469896040 00 
1 • 53 861 0 OOD 02 


1.119155000 02 
-1 .6927439 30-0 1 
—3 .831350020 00 
9.936690950-01 
8 .167061 320-01 


1 .009647160 00 
6 .272922520-01 
3.1 3S1 58230 00 
9 .866692860-01 
9 .84274 595D- 01 


6.598215690 02 
7 .072549660-02 
9.491 1 1739D-02 
0.0 

4.775871220 02 


SWITCH THRUST OFF 


7.496209950 02 
4 .07334674D-02 
-6 .461999620-01 
-2.8358934*0 01 
-8 .1 30258760 01 


1 .008048750 00 
-8.540298070-01 
8. 6 99 0 0 64 3D— 0 1 
-1.095770480 00 
3.625483640 01 


1 .32 2 28 6 5 7D — 02 
-5 .16 4 97 04 1 D— 0 1 
-2.032552950 00 
0 .0 

8 .299592590 01 


3 . 669791 4 1 D 01 
9.313290373-01 
I .22 0 096 7 B D 00 
1 .0723421 50 02 
-3.1 13573320 01 


3.858409670 01 
2.343064 76D — 01 
-1.967002170 00 
3.4767547 2D 02 
—8.726931 26D 01 


2.399107740 02 
—2. 961 16285D-01 
3.546784520 00 
1.003928530 00 
-5.50 05384 70- 01 


9.908970520-01 
5.2 54455 04 D — 0 1 
4.329961250 00 
1 .000000000 00 
1 .005083490 00 


7.869761060 02 
8.557494710-02 
9.491 1 1 731D-02 
-1 .332267630-15 
6.094203840 02 



time 

R1 
LI . 
LG 
PSI 


SEMI-MAJOR 

AXIS ECCENTRICITY 

inclination 

R2 

R3 

VI 

L2 

L 3 

L* 

LC 

LPHI 

CONE 

THETA 

PHI 

LATITUDE 


NODE 

ARG P05 

RMAG 

V2 

V 3 

MASS RATIO 

L5 

L6 

L7 

CLOCK 

HMAG 

POWER FNCT 

LONGITUDE 

FLT PTH ANG_E 

VMAG 


Q 


O 

CO 


A ^ 
ci 

. & 

h-3. ,..-15' 
02“' 


£ 


SWITCH THRUST ON 


8 *06 9674 900 02 
7 *973548400—01 
7.83S6I7570-01 
-2.83S09344O Ot 
7 *946959820 01 


INPUT TARGET 

9 .125000000 02 
3.00751 3090-0 1 
1 .21941 148D 00 
—3 .4 70770 BID 01 
8.38 1241470 0 1 


I .008048750 00 
-2.7212591 10-01 
-1 * 4294343 2D 00 
-1 .09S77048D 00 
-6.063345320 01 


1.00106538D 00 
8,246264640-01 
-1.874460100 00 
-1.512366970 00 
6.962279920 01 


1 .32 22065 7D- 02 
-S.29181042D-O1 
1 .627366430 00 
0 .0 

8.6! 8216410 01 


9 .482929360-05 
4 .807877500-01 
2.728663000 00 
0 .0 

8.78491690D 01 


3.669791413 01 
4.738377590-01 
1.502294470 00 
1 .135914070 02 
-3.213291830 01 


4.499978360 01 
-8.543975790-01 
-1.264304910 DO 
7.320074723 01 
2 . 87062 7 3 BO 31 


3.858489670 01 
8.478079300-01 
-2.195337040 00 
1.956872790 02 
-1.39440461D 01 


3.576237590 01 
9.159705550-03 
1.820503340 00 
2.321395460 02 
6.996698860 01 


2.971215400 02 
2 .7368 DSS2D-0 1 
3.286792330 00 
1 .003928530 00 
1 .300956420-01 


4.278717760 01 
5.196383550-01 
— 2.369610S3O 00 
1.000 532 54D 00 
-1 .594932950-03 


9.949174200-01 
5.254455040-01 
4.329961250 00 

i .eoooooooo oo 

1.00905975D 00 


1 .000974630 00 
4.428474530-01 
5.782196730 00 
9.986459540-01 
9. 99S 583450-01 


END OF TRAJECT DRY, THRUST ON 


PAGE 3 

TRAVEL 
THRUST ACC 
HAM 

SWITCH FNCT 
PROP TIME 


8,441868720 02 
8.557494710-02 
9.491 117310-02 
-2.220446050-16 
6.094203840 02 


9.482619080 02 
1.01 398S43D-01 
9.491117590-02 
9.355 19994D-01 
7.149528940 02 



CASE X 


ITERATOR SUMMARY 


INDEPENDENT PARAMETERS 


1 »PRIM1 ( 1 .4330 221 D-03 > 

6 ,PDOT3 (-3. 264391 60-01 ) 
IX » ACCEL ( 2.66648050-04 1 

1 6 ♦ T 1 ME 2 ( 9.4353700D 02 1 

21 .THET1C 0.0 » 

26 »THET6< 0.0 > 

31 . PH 1 l { 0.0 > 

36, PH 1 6 ( 0.0 1 

41.PR1-AI 0.0 ) 

46 ,P03-A< 0.0 ) 

51 « PR 1 — 0 C 0.0 ) 

56 . P03— B ( 0.0 ) 

61 ,PR1-C< 0,0 1 

66,PD3-C( 0.0 » 


1 . OELT A X< 3.625590-05 
6, DECT ZD( 4.21362D-05 
11, POWER ( 2.03500D 01 


2.PR1 M2 I -3. 71 2734 70-01 
7.LMASSC l.OOOOOOOD 00 
12, V J ET ( 2.9419950D 04 

17,1 PARK t 0.0 
22.THET2 < 0.0 
27.T HET7 ( 0.0 

32. PHI2( 0.0 
37, PHt7< 0.0 
42,PR2-A( 0.0 

47 , V INF A l 0.0 
5 2 , PR2-B ( 0.0 
57 , V [ N FB ( 0.0 

62,PR2-C< 0.0 
67 , V INFC ( 0.0 


2 .DELTA Y( — 1 * 85804D- 05 


7 , 

< 

0.0 

> 0 , 

12 . 

< 

0.0 

) 13 ,T 

17 . 

< 

0.0 

) 18 , 

22 . 

( 

0.0 

) 23 . 

27 . 

< 

0.0 

) 28 » 

32 , 

( 

0.0 

> 33 , 

37 . 

( 

0.0 

1 38 . 

42 , 

( 

0.0 

) 43 , 

47 , 

( 

0.0 

) 48 , 

52 . 

( 

0.0 

) S 3 , 

57 , 

( 

0.0 

> 58 . 

62 . 

( 

0 . 0 

) 63 , 

67 , 

( 

0 . 0 

) 60 . 


} 3,PRI«3{ -5.1922233D 001 

>) a, ltaj< o.o I 

.) 1 3 » VI NF I ( 5.0 99 6 2 6 5 D 03! 

) 18 . VE . 31 I 0.0 i 

> 23 .THE T 3 ( 0.0 

> 28,THET3( 0.0 

> 33, P-U3( 0.0 

) 3B . PH 18( 0.0 

) *3,PR3-4( 0.0 

) 48 ,T I 4EA( 0.0 

I 53.PR3-3C 0.0 

J 5b.tim=3( 0.0 

> 63.PR3-C1 0.0 

J 68 «T I MEC( 0.0 

DEPENDENT PARAMETERS 

i) 3.0ELT4 Z(-4. 362440-05 

> 8 , (0.0 


4 , ROOT 1 ( 7.06559150-01) 
9* ( 0.0 ) 

14.VIMF2C 0.0 ) 

1 9 • VEuO 2 ( 0.0 ) 

24,THcT4( 0.0 ) 

29. T MET 9 ( 3.0 ) 

34, PH I 4 ( 0.0 ) 

39, PHl9( 0.0 > 

44,P0I-A( 0.0 > 

49.KSAMP{ 0.0 > 

54 . PD1 - b ( 0.0 ) 

59 , KSAMP l 0 . 0 ) 

64 , P D1 -C ( 0.0 ) 

69, KSAMP ( 0.0 ) 


4, DELT KD( -4 . 2001 8D-C5 


5.POOT21 6,35286350-01 
10. OECLNI-3. 785367*0 01 
1 5, TI ME1 ( 3.10370000 01 
20 . VEL03 ( 0.0 
2 5* T HET5 ( 0.0 
30 , LOEGR ( 0.0 

35. PH15( 0.0 
40.PHI1 0( 0.0 

45. P02- A ( 0.0 

50 . KDRQP ( 0.0 
55 . PD2 — B ( 0.0 

60.KDRdP( 0.0 
6 5.P02-C ( 0.0 

70,KUROP( 0.0 


5. DELT YO( 2.217760-0* 
10 . T * D£CLN( 5. 7 7598 0-07 

15, ( 0.0 

20. ( 0,0 

25. ( 0.0 

30. ( 0.0 

35, ' < 0.0 

*0, ( 0.0 

45, ( 0.0 

50. < 0.0 

55. ( 0.0 

60. (0.0 
65. ( 0.0 

70. ( 0.0 


THRUST SWITCHING TIMES (DAYS) 

4 39.995 ON 568.1 92 OFF- 


0.0 ON 
617.788 ON 


80.338 0=F 

749.621 OF ~ 


103.587 ON 
836.967 ON 


POWER 

20.3499990796 


INITIAL 

3268.5*60092059 


EFFICIENCY 

0.6300000000 


PROPULSION 
620,2499 71 9473 


ELECTRIC PROPJLSION PARAMETERS 
PROP TIME J 

714. 95 2893 03 45 9.3696171166 

MASS COMPONENT BREAKDOWN 
PROPELLANT TANKAGE 

1821 .0787330785 6 3.7377556577 


220.047 OFF 
912.530 ON 


P*OP T IME RAT IO 
0.7835100206 


STRUCTURE 

0.0 


*01.007 OFF 


AVE ACCEL 
0.0004006926 


PAYLOAD 

763.4795485224 


SWITCH-COUNT HISTORY ALL 12 

*22 THRUST COMPUTE STEPS. 32 COAST COMPUTE STEPS 
LAUNCH ASYMPTOTE OFFSET FROM PRIMER = -32.786 DEGREES. 



ORIGINAL PAGE 10 
OF POOR QUALITY 


CASE 1 


EXTREMUM POINTS 3 s SELECTED FUNCTIONS 


X 

time 

ECLIPTIC 

SOLAR 


COMMON IC AT I ON 



SW ITCH 



LONG 1 TUDE DISTANCE 


ANGLE 

DISTANCE 


FJNCTION 

0 

0.0 

0.0 


1.005 


75.7 


0 

.0 

ON 

4 

.1 9D 00 

4 

1 .065 

1.0 


1.004 


75.3 


0 

.30 

MAX 

4 

.1 90 00 

5 

8.379 

7.6 

* 

1.000 


72.6 


0 

.03 


4 

.1 50 00 

4 

0.428 

7.7 


1.000 


72.5 


o 

.03 


4 

.1 5D 00 

4 . 

60.338 

B2.2 


0.809 


52.3 


0 

.24 

OFF 

-4 

.0 0D-15 

6 

01 .933 

96.6 


0. 870 


51 .4 


0 

.28 

MIN 

-7 

• » 5D-01 

T 

99.677 

106.6 


0.650 

MIN 

51.3 


0 

.31 


-3 

.230-01 

5 

103.587 

111.8 


0. 853 


51 .3 


0 

.32 

ON 

-1 

.330-1 5 

4 

136.984 

156.9 

min 

0. 830 


54 .2 


0 

.45 


2 

.3 90 00 

4 

1 59.653 

166 . B 


0.840 


57.0 


0 

.55 

MAX 

2 

. 9 3D 00 

S 

165.388 

194 . 2 


0. 846 


57.6 


0 

.57 


2 

.900 00 

4 

220.047 

259.0 


0.930 


62.3 


0 

.78 

OFF 

“1 

.750-14 

S 

234.249 

27 5. 5 


0.952 


63.3 


0 

.81 

MIN 

“7 

-0 9D-01 

5 

248.415 

290.5 


0.573 


64.5 


0 

.82 

GN 

2 

.390-15 

4 

257. 189 

299. S 


0. 904 


65.3 

max 

0 

.82 


6 

.500-01 

4 

272.871 

31 4.9 

* 

1.000 


66.7 


0 

.82 


1 

.740 00 

5 

315.465 

353.8 


1.019 


70.4 


0 

.74 


3 

.470 00 

4 

3 t 6.969 

355. 2 

MAX 

1.019 


70.5 


0 

.74 


3 

.490 00 

5 

325.138 

362.3 


1.018 


70.9 


0 

• 72 

MAX 

3 

• 5 3 D 00 

4 

335.862 

371 . B 


1.016 

MAX 

71 .1 


0 

.70 


3 

.4 50 00 

4 

356.876 

390.7 


1 . 006 


70.3 

MIN 

0 

.69 


2 

.340 00 

3 

366. S44 

399.0 

* 

1.000 


69.5 


0 

.59 


2 

.390 00 

6 

401 .007 

434 . 9 


0.976 


65.2 


0 

.75 

OFF 

~2 

.2 20-16 

6 

420.432 

456.6 


0. 965 


62.8 


0 

.80 

M IN 

-1 

.240 00 

5 

439,995 

478.9 


0. 957 


61.1 


0 

.84 

ON 

-1 

.4 20-1 4 

4 

461.762 

503. 1 

min 

0. 953 


60.2 


0 

.07 


1 

.520 00 

4 

467.030 

508. 7 


0. 954 

MIN 

60.2 


0 

.97 


1 

.3 10 00 

5 

476.946 

51 0.9 


0.955 


60.3 

MAX 

0 

.97 


2 

,270 00 

4 

S03. 064 

544.6 


0. 962 


61 .2 


0 

.97 

MAX 

2 

.790 00 

6 

509.882 

551.0 


0.965 


61 .5 

MIN 

0 

.07 


2 

.750 00 

5 

518.520 

559.2 


0.969 


61.0 


0 

.97 


2 

.510 00 

4 

529.847 

570. 1 


0.974 


62.1 


0 

• B 7 


2 

» 2 6D 00 

4 

544 .170 

584.0 


0.900 

MAX 

62.3 


0 

*99 


1 

» 5 3 D 00 

4 

568.192 

600.0 


0.992 


61 .9 


0 

.94 

OFF 

“2 

.Q0D-1S 

5 

588.100 

630.6 

• 

1.000 


61 .3 


0 

.97 


-1 

.510 00 

6 

S93.002 

636. 1 


1.002 


61.2 


0 

.90 

MIN 

-1 

.720 00 

5 

604.245 

640.5 


1. 006 

MIN 

61 .1 


0 

.99 


-1 

.040 00 

5 

613.892 

658.9 


1. 009 


61 .2 

MAX 

1 

.00 


-2 

.920-01 

6 

617.708 

663.0 


1.010 


61 .3 


1 

.30 

ON 

0 

.3 

5 

633.894 

679.4 


1 .01 3 


62.1 


0 

.98 


1 

.1 OD 00 

5 

651 .984 

696.2 


1. 01 5 


63.6 


0 

.94 


2 

•04D 00 

4 

6S8. 478 

701.9 

MAX 

1.015 


64.3 


0 

.92 


2 

• 2 SD 00 

4 

668.726 

71 0.5 


1.014 


6S.4 


0 

.90 


2 

.560 00 

6 

603.476 

722.6 


1.013 


66.6 


0 

.93 

MAX 

2 

.7 1 0 00 

7 

692.955 

730.2 


1.01 1 


67.4 


0 

.91 


2 

• 5 5 D 00 

4 

704.354 

739.6 


1 .009 

MAX 

67.7 


0 

. 79 


2 

,410 00 

6 

711.876 

745.9 


1.007 


67.6 

MIN 

0 

.78 


2 

.ISO 00 

4 

743.566 

775.5 

* 

1 .000 


64 .5 


0 

.85 


4 

.290-01 

4 

749.621 

701.9 


0. 999 


63.6 


0 

.07 

OFF 

-1 

• 33 D-I 5 

6 

778.265 

81 5.S 


0.995 


59.6 


0 

.99 

MIN 

-2 

.000 00 

9 

796.940 

03 B , 7 

MIN 

0.995 


50.1 


1 

.03 


-7 

.51 D— 0 1 

4 

806.967 

050.3 


0.995 


57.6 


1 

.04 

ON 

-2 

. > 2D-16 

9 

609.1 71 

052.0 


0.995 

MIN 

57. B 


1 

.04 


1 

.570-01 

4 

010.095 

053.0 


0.995 


57.6 

MAX 

1 

.04 


2 

.220-01 

4 

857.024 

096.6 


1 . OOO 


62.1 


0 

.92 


2 

.290 00 

5 

663.1 33 

900.6 

• 

1.000 


62.8 


0 

.90 


2 

.330 00 

4 

666.463 

903-1 


1.000 


63.3 


0 

.99 

MAX 

2 

.340 00 

5 

6 60 .26 2 

904 . 4 


1.000 


63.5 


0 

.98 


2 

.340 00 

6 

895.587 

924.9 


1.001 


65.7 

MIN 

0 

,82 


1 

.750 00 

5 

096.651 

925.0 

MAX 

1.001 


65.7 


0 

.83 


1 

.71D 00 

4 

698.104 

926.9 


1.001 

MAX 

65.7 


0 

.83 


1 

.550 00 

S 

912.300 

939.1 


1.001 


65.1 


0 

.85 

ON 

9 

.3 5D-01 




THRUST ANGLES 



input 


ARRAY 


PS I 


Theta 


PHI 


POWER 


ANGLE 


-77. 1 


bo. 5 


87 .9 


20.2 

ON 

0.0 


-77. 1 


80.5 


07,9 


20.2 


0.0 


-77.2 


79.6 


87.7 


20.3 


0.0 

MIN 

-77.2 


79.5 


87 .7 

MAX 

20.3 


0.7 


-60.4 


44.8 


69.5 


20.3 


37.7 


4*4* 4 


*44444 


4*4 44 


0.0 


90.0 


4*44 4 


*44444 


4**4* 


0.0 


90.0 


73.2 


-3.5 


73 .3 


20.3 


43.4 


05.4 


-80.5 


09.2 


20 .3 


46. S 


07. 2 


-133.4 


91 .9 


20.3 


45.1 

MAX 

07.3 


-154.1 


92 .4 


20.3 


44.4 


64.0 


108.4 


98.0 


20.3 


30.2 


44*44 


*44444 


44*** 


0.0 


90.0 


-50.5 


99,2 


94 .8 


20.3 


18.8 


-66. 1 


97.2 


92 .7 


20.3 


14. S 


-74.7 


93.7 


91 .0 


20.3 


0.0 

MIN 

-70. 1 


84.3 


88.8 


1 9. 9 


0.0 


-78. 1 


64.0 


88 .8 

MIN 

19.8 


0.0 


-78. 0 


82.1 


SB .4 


19. a 


0.0 


-77. 7 


79.5 


07.0 


19.9 


0.0 


-76.5 


73.6 


06.2 


20.2 


0.0 


-75. 6 


70.3 


. 65.2 


20.3 


0.0 


-66. 1 


49. D 


74.6 


20.3 


17.8 


*444 4 


4*4*44 


* 4* ** 


0.0 


90.0 


76. 0 
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£5132a5i22I^_fi2_iULlAa_eai£ depart earth 


x 


Y 


Z 


X03T YDOT 


ZDOT 


RADIUS LAT. 1.0NG* 


PLANET -8 .6292634D-01 -5 . I 5 1 1 9080- 01 0.0 

S/C -8.6292634D-01 —5.1 5119 OSD- 01 0.0 


4. 96313420-01 
4.9S43363D-0I 


-B. 62411130-01 
-7, 80199270-01 


0.0 

-1.50135950-01 


1*00498240 00 
1. 0049824D 00 


0.0 -149.165 

0.0 -149.165 


CASE 1 t CONVERGED 1 


PERFORMANCE SUMMARY 


-o 

to 


OUT OF ECLIPTIC MISSION 


FINAL INCLINATION » 44.9998 OEG 


LAUNCH VEHICLE IS TITAN III E/CENTAUR 


(COEFFICIENTS = 167238.9500 3460.2038 1753.6965) 


LO a APR 21 ( 1979. 12.8880 HOURS GMT 

JULIAN DATE 43985.0370 


AD 


OCT 20. 1981, 

JULIAN 


0.8880 HOURS GMT 
DATE 44897.5370 


FLIGHT TIME a 912.5000 DAYS. 


ELECTRIC PROPULSION SYSTEM PARAMETERS 

ALPHA A (KG/KW) ALPHA T ( KG/KK ) TANKAGE FACTOR STRUCTURE FACTOR 

15.00 00 15 .0000 0.0 3S0 0.0 


EFFICIENCY COEFFICIENTS 
8 D ( KM/SEC > 

0.63000 0.0 


E 

0.0 


INITIAL 
32 66.54 6 0 


ELECTRIC PROPULSION SYSTEM MASS SUMMARY <KG) 


POWER PLANT PROPELLANT 

620.2500 1821.0737 


TANKAGE STRUCTURE 

63.7378 0.0 


NET MASS 
763.4795 


PH AU) (KW> P(HSKP) ( K W ) 

20.3500 0.6500 


ELECTRIC PROP UL5 1 DN SYSTEM PERFORMANCE SUMMARY 

P(TARG) CKW) THRU AU i IN) ACC 1 1 AU > (M/SEC**2) ISP ISEC) EFFIC CHAR OEG (DAYS) 

20 .3224 0.8 71 55 1 2.66648 0 0-04 30 0 0.0 00 0*63000 1 .000000CD 30 


MAX D 1ST (AU) MIN DIST ( AU 1 

1 .0X89599 0.8296067 


EXTREME TRAJECTORY ANO PERFORMANCE CONDITIONS 

MAX POWER ( KW) MAX TRRUST (N) BURN TIME (DAYS) 

20.349999 0,37155141 714.95289 


DEGRO TIME (DAYS) TRAV ANG (OEG) 
710,01587 948.26191 


DEP DECL (DEC) 
-37.8S37 


DEPARTURE AND ARRIVAL CONDITIONS 

PARK INC (DEG) DEP VI NF CM/SEC) C3 ( KM* *2/SEC* * 2 ) ARR VINF (M/SEC) 

36.0000 5099.62649 26.006211 0.0 

AT LIMIT 


C4 <KM**2/SEC**2) 

0.0 


B. COMET RENDEZVOUS MISSION 

The objective of this mission is to deliver maximum payload to the comet 
Tempel II, rendezvousing at perihelion of the comet f s path in the 1988 apparition. 
For the case shown, the arrival date is fixed at September 16, 1988; the launch 
date and launch excess speed are optimized to yield maximum payload for the 
nominally 4 -year class of solutions. The launch vehicle assumed is the Titan 
III E/Centaur, The electric propulsion system parameters are representative 
of projections of the SERT III spacecraft. Specifically, the reference power is 
8,671 kw, the specific impulse is 2900 seconds, the efficiency is 0.63376 
and the specific propulsion system mass (arrays plus power conditioning and 
thruster subsystem) is 30.285 kg/kw» 

This example illustrates the use of the solar array degradation option. 

This option is invoked by setting the characteristic degradation time T POWER 

10 

to a value less than 10 . In this case, it is set to 7121 days which means 

that if the arrays were situated at 1 AU and oriented normal to the sun line for 
this amount of time, the power developed by the array would degrade to l/e of 
its initial power output. The input MPOW = 1 forces the arrays to be oriented 
normal to the sun line throughout the mission. Setting the triggers X30(2) and 
Y30(2) to 1 results in the optimal adjustment of the trajectory to accomodate 
the degradation. The complete set of inputs for this case follows and the resultant 
program output begins on the next page. 


&M I NPUT Xl(2)«l.DO,X2(2) !a l.DO,X3(2)=l.nO / X4(2)=l.DO,X5(2)“l.DO 
X6(2) S 1.DO,X11(2)=1.DO,X13(2)=1.DO,X15(2)=1.DO,X7=1.DO 
Y1C2) =1.0 0, Y2( 2) =1.0 0, Y3( 2) =1.D0,Y4 (2) =1.00,75(2 ) =1.0 0, Y6( 2) =1.0 0 
Y13(2)=l. DO, Y15(2)=l. DO, Al PH AT = 15. 285DO, MBOOST^l 5, MTM AS S=3 
MO PT2 = 3,M0 PT3 =4 3,MYEAR=19 8 8, MO NTH = 9 ,MDAY=16,MOUR=17. 2 20 0 
X12=2.8439285D4,B{=.63376D0,DI=0.D0 / CTANK=. 1D0, Yll =8 . 6 71D 0, 3. DO, 1. 0- 3 
TPOWER=7121.DO,MPOW = 1,X3 0( 2) =1,00,73 0(2) =1.0 0 
Xl= 5.0729985882240 00, X2- 3 . 3 501 836 72 12 6D 00, X3=- 2 . 1 59 80 814 2 3 SOD 00 
X4 s -2. 06 931493331 2D 00, X5= 4 . 56 5 24 774 1 305D- 00, X6= 3 , 317322088 503D-01 
Xll= 1.64797965415 ID -04, XI 3= 6 . 718 039 79 849 3D 03, X15=-l, 5176 252844 70D 03 
X30=-4. 6574062226460-03 &END 
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Y 3 T 

a 

0.0 

* 0.0 


Y 38 

a 

0.0 

. 0 .0 


Y 39 

a 

0.0 

* 0.0 


Y 40 

3 

0.0 

, 0.0 


Y 41 

a 

0.0 

, 0.0 


Y 42 

a 

0.0 

, 0.0 


Y 43 

3 

0.0 

, 0.0 


Y 44 

a 

0.0 

, 0.0 


Y 45 

a 

0.0 

, 0.0 


Y 46 

a 

0.0 

, 0.0 


Y 4 7 

3 

0.0 

, 0.0 


Y 48 

3 

0 .0 

. 0.0 


Y 49 

M 

0.0 

, 0.0 


Y 50 

8 

0.0 

• 0.0 


YS 1 

3 

0 .0 

, 0.0 


Y 52 

a 

0.0 

* 0.0 


Y 53 

w 

0.0 

, 0.0 


Y 54 

3 

0.0 

, 0.0 


YS 5 

a 

0 .0 

,• 0.0 


V 56 

a 

o 

o 

. 0 .0 



3. OOOOOOOOOOOOOOOD 00 

3.0000000000000000 00 

3. 000 0 00 00000 0000 D 00 
3.0003000000000000 00 
3. 000 OOOOOOOOOCOOO 00 
3. OOOOOOOOOOOOOOOD 00 

5.0000000000000000 02 
5. OOOOOOOOOOOOOOOD 01 
3. OOOOOOOOOOOOOOOD 00 
3. 033 30000030 OOOOD 00 
9. 999 9 9 99 9 99 99 999D -0 5 
9. 999 9 999 99 99 9 9990-0 5 
9. 9999999999999990-05 
9. 999 9 999999999990-0 5 
9. 9999999999999990-05 
9. 999 9 999999999990-0 5 
9.9999999999999990-05 
9. 999 9999999999990-05 
9.9999999999999990-05 

9. 9999999999999990- 05 

9. 9999999999999990- 04 

9. 9999999999999990- 05 
9. 9 99 99 99999999990-05 
9. 999 9 999999999990-05 
9. 999 9 999 99 99 9 99 9 0-0 5 

9. 9999999999999990- 05 
9. 999 9999999999990-0 5 
9. 9 99 9 99 99 9 99 99 9 9D-0 5 

9. 9999999999999990- 05 
9. 999 999999 9999 99D-0 5 
9. 9999999999999990-05 
9. 999 9 999999999990-0 5 
9. 999 999999999999D— 0 5 
9. 999 9999 999999990-0 5 
9. 9999999999999990-05 
9 . 999 9 9999 9 999 99 9D -0 5 
9. 999 9 999999999990-05 
9. 9999999999999990-05 
9. 999 9999999 999990-05 
9. 9999099999999990-05 
9.9999999999999990-05 
9. 9999999999999990-05 
9. 999 9 999999999 990-0 5 
9. 999 9 99999999999D— 0 5 
9. 999 9999999999990-0 5 
9. 999 9 999999999990-0 5 
9.9999999999999990-05 
9 . 9999999999999990-05 
9.9999999999999990-05 
9.9999999999999990-05 
9.9999999999999990-05 
9. 999 9 999999999990-05 
9. 9999999999999990-05 
9. 9999999999999990-05 
9. 999 99999999 99990-0 5 
9.9999999999999990-05 
9. 999 999999999999D-0 5 
9 . 999 999999999999D-C5 
9. 9999999999999990-05 
9. 9999999999999990-05 
9. 9990999999999990-05 
9.9999999999999990-05 
9. 999 999 99999999 9 0—0 S 
9.9999999999999990-05 
9. 999 9 999999999990-05 
9. 999 9999999 999 950-05 


I • 00 OOO OOO 0 U 00 0 OOD— OH* 
1 .OOOOOOOOOOOOOOOD-08, 

i .oooooooooooooooo-oa, 

1 .00 OOCO 000000 OOOD-OB, 

1.00000000 00 000 00 0- oa, 

1. 0000000003000 OOD- 08. 

9. 9999999999999990- 05. 

9.9999999999999990- 07. 

1.0000000000 000 00 0- oa, 

1 . 0000000000000000 - 08 . 


i .ooobooo 000000 0013 00 . 
1 . 0^00000000000000 00 . 
1 .OOOCOOOOOOOOOOOO 00, 
1,0000000000000000 00. 
1 .OOOOOOOOOOOOOOCD 00. 
1.0000000000000000 00 . 
1 . OOOOOOOOOOOOOOOD CO. 

: .oooooooocooooooo oo. 
1.0000000000000000 00, 
1.0000000000000000 00 . 



<J3 


Y57 » 0.0 , Q .0 

Y58 *0,0 ,0.0 

YS9 = 0.0 , 0.0 

Y 60 =0.0 , 0.0 

Y6I =0.0 ,0.0 

Y62 = 0.0 . 0.0 

Y 63 = 0.0 . 0.0 

Y 64 = 0.0 , 0.0 

Y6S =0.0 .0.0 

Y66 = 0.0 . 0.0 

Y67 * 0.0 , 0.0 

Y 63 = 0.0 , 0.0 

Y 69 =0.0 ,0.0 

Y 70 = 0.0 , 0.0 

tOUTPUT 


9.9999999999999990- 05, 

9. 9999999999999990- 05. 

9. 9999999999999990- 05 . 
9. 999 999999999999D— 0 5 , 

9.9999999999999990- 05, 
9. 999999999999999D— 05 . 

9.9999999999999990- 05, 
9. 9999 9999999999 9D-0 5 , 

9. 9999999999999990- 05. 
9. 9999999999999990-05. 
9. 9999999999999990-05. 

9. 999 9 9999999999 90-0 S , 

9.999 9999999999990-0 5, 
9. 9 99 9 9 9999 99 9999 D-0 5 . 


AAI= 0.0 
0.0 
0.0 
0.0 
0.0 

CTRET= 0.1 1 1111 tllll Ml 

0.0 

0.0 

0.0 

0.99999999999999 990-0 A 

0.0 

PS1GN= 1 . 0000000000000 

0.0 

RPER= 2.00000000000000 
0.0 
0.0 

STATE* 1.0000000000000 
0.0 
0.0 
0.0 
0.0 
0.0 

-1 .000000000000000 
-1 .000000000000000 
-1 .000000000000000 
-1 .000000000000000 
-1 .oooooocoooooooo 
0.0 

1 .000000000000000 
250.000 0 COO 0 00 0 00 
5 00. 0000 CO 0000000 
0.0 


00 .RAOOOOS 1.000000000000000 ,R ADQD 

* 0.0 , 0.0 

0 ,SAl= 1.000000000000000 ,SAIX= 0.0 

. 0.0 , 0.0 

. 0.0 , 0.0 

00 , 0.0 , 0.0 


* AE* 0.0 

AN= 1 . 500000000000000 

0.0 

B3= 0.0 

0.0 

1 , Dl = 0.0 

EC 1 X= 0.0 
El = 0.0 
0.0 


, ALPH A A= 15.00000000 00 3000 , AI_PHAr= 1 5 .28500000 0000 0 0 

,AR= 1.000000000000000 . ASOL* 0.0 , 0.0 


■ ALT I TO® 


0.0 

» CMI = 0.0 

0.0 

,oviRer^= o.o 

, o ,o 

, EMUOOD = 0.0 

0.0 


.61 = 0.6337600000000000 
, CN I X= 0.0 
,CSTR= 0.0 

• DPO *= 0.0 

« 0.0 

,EMJODX= 0.0 


,B1= 0,0 .02= 

, 0.0 , 

»CTANK= 0.10000000000000000 00, 

» EC I = 

, 0.0 « 


HOUR* 17.22000000000000 .DMI= 0,0 

0.0 , 0.0 


,OMlX= 0.0 


, it A OODX= 0.0 


,GAHMAX= 1.000000000000000 


» 0.0 .POWFIXs -1 .000000000000000 

.0 , 0.0 , 

,RAP= 36.00000000000000 ,R£VS= 0.0 

, 0.0 , 

.SOI = 0.0 ,SOIX= 0.0 

,0.0 , SPIRET® 300.0000000000000 

, 0.0 . 1 .000000000000000 


STEP1= 0. 31 250 000000000 0 OD-Ol , S TEP2 = 3.12500000 0 00000 00 

0.0 , 0.0 , 0.0 

0.0 , 0.0 , 0.0 

0.0 . 0.0 , 0.0 

0.0 ,0.0 ,0.0 

THRET® AOO .0000000000000 , T QF F a -1 . 000003 OO 0 DO 0 03 0 


, T COAST = 0.0 , 

, 0.0 • 

, 0.0 • 

• 0.0 , 

,T0V= -1.000000000000000 ,TGO® 
-1.0000003000 00 000 , -l . 00 00 00 00 00 00 000 


0. ITF= 
O.LOADXa 
O.MOPTX® 
0. 


-1 .000000000000000 
-1 .000000000000000 
-1 .000000000000000 
-1 .030000000000000 
0.0 

T 2- 50.0 00 000 00000 00 0 

300.0000000000000 

XANG1® 0.0 
0.0 
o.o 
0.0 
0. 0 
0.0 
0. 0 
0.0 

IBAL* O.INTPR 

3 , 1 TPRNT a 0 

0. M AXHAM= 

0, o, 

0. 


. -I .000000000000000 
, -1 .000000000000000 
, -1.000000300000000 
» TP 1 = 0.0 
, 0.0 


. -1 .000000000000000 

, -1.0000 000 oooooooo 

, -1.030000000000000 
»TP I X= 0.0 


t -1,000000000009000 
« -1 .COCOOO 000000000 
• - 1 .000000000000000 

. 0.0 


,TPO*ER= 71 21.00000000000 0 .TSCAL.E® 


• 100.0000000000000 , 150. 000 00 00000003 » 200.0000000000000 

350. OOOODOOO 30D00 , 400.0000000000300 . 450.0000000000000 

« X ANG2 = 32.5 0000000000000 ,XSV*InG= 0.0 , 


. 0.0 
, 0.0 
. 0.0 
. 0.0 
, 0.0 
. 0.0 
. 0.0 

0. INTPR® 

0, JPP= 

5,MS00ST= 

O, 0, 

0 , 


» 0.0 
, 0.0 
. 0.0 
. 0.0 
. 0.0 
• 0.0 

• 1 .030001500000330 


, 0 . 0 
. 0.0 

• 0.0 
, 0.0 
, 0.0 

,X0= 1,000000000000000 

, 0.0 


o, iour= 

0 . JPRINT® 

15 » MDAY = 

0 , 

0 , 0 , 


, I RK= MILS 1, 

0,JT= O.KPARTa 

16. MODE® 4. MONTH® 

0.M0PT2® 3.MOPT3® 


1 , IROT = 0.1 

0, LAUNCH® 

N» 9, MOPT= 

4 3. MOPT 4 s 


0, MPUNCH= 0 

0 , 0 , 0 , 

3 , NMUNGa 2 S. NORMAL* 

0, 300 , NSW ING® 


O.MREAD® O.MSrflNG® Oi 

0, 3.HTMA5S® 

= O.NPERFa 0 , NPR I NT = 

O.NSWPAR® 1 , NTAPE® 


3 , MUPO AT = 

7 . Nser = 

1 7 


O.HPOrfa 1, 

o. o , o ■ 

1 »MYE AR= 1 966 .NO I ST® 

0. 0, 300. 


fcEND 



'JNAL page is 
GG ii QUALITY! 


c; 

fa 


O 
.rj 
■ -< 


CASE 1 


ITERATOR PARAMETERS 


NO. INDEX 

1 I 

2 2 

3 3 

A 4 

5 5 

6 6 

7 11 

0 13 

9 IS 

10 30 


VALUE 

3.07299856022400000 00 
3 .350 163672 1260000D 00 
■2 .1 596081 423690 000D 00 
■2.06931493331200000 00 
4.565247741 30500000 00 
3.31 7322 C 685030 0000-01 
1 .64 79 79 654 151 00 0 00-04 
6.71803979649299900 03 
1 .51762S2644700000D 03 
4 .657406222645999 00-03 


INDEPENDENT 
STEP LIMIT 

3. OOOOO COO OOOOOOOOD 00 
3.00000000000000030 00 
3 • OOOOOOOO DOOO OOOOD 03 
3.00000 00 OOOOO 00 OOD 03 
3 . 000000000000 OOOOD 00 
3 .00000 OCOOO OOOOO OD 3 3 
9. 99999 9999999999 30-04 
5.00 000 00 OOOOO 00 OOD 02 
B. 00 00 0000 0000 OOOOD 00 
1 . OOOOO OOOOOOOOOO 3D 01 


variables 

PERTURBATION 

l . oooooooooo ooooooD-oa 
1 . OOOOOOOOOOOOOOOOD-O 8 
1.0000000 00 000 00 OOD- 08 
1.000 000 000 OOOOO OOD- 06 

1 . 000 0 00 000 OOOOO OOD-0 B 
1 .00 00 000 00 OOOOO OOD- 08 

1 .00 0 00(3000 0 0 000 000-1 1 

9.99999999999999900- 05 

9.99999999999999900- 07 

9.99999999999999900- 05 


WEIGHT 

1 .000 OOOOO OOOOOOOOD 00 
1 .000 OOOOO OOOOOOOOD 00 
1 .00000000000000000 00 
1 .0000 00 OCOOOOOOOOO 00 
1.00000000000000000 00 
1 .000 OOOOO OOOOOOOOD 00 
1 .3000 00 OCOOOOOOCOD 00 
1 .0000 OOOCOOOOOOOOD 00 
1 .00000000000000000 00 
1 .00000000000000000 00 


-0 

-3 


NO. INDEX 

1 1 

2 2 

3 3 

* 4 

5 5 

6 6 

7 1 1 

e 13 

9 IS 

10 30 


DEPENDENT VARIABLES 


VALUE 

0.0 
0.0 
0.0 
0.0 
0.0 
O .0 

8 .67100000000000000 00 
0 .0 
0 .0 
0 .0 


TOLERANCE 

9. 99999 999 9999 999 3 D-0 5 
9. 9999999999999993 D-0 5 
9.99999999999999930-05 
9. 99 9 99 999 9999 99 9 3 D-0 5 
9.99 99 9 599 999 7 99930-0 5 
9. 99 999 999 9 9 99 999 00-05 
9. 999999999999 9993 0-0 4 
9. 99999 999 5999 99933-05 
9.99999999999999930-05 
9. 99 999 999999999 3 3 D-0 5 



INHIBITOR IS 


5.8208D-11 


NOMINAL TRAJECTORY 1 {TOTAL 1> 





INDEPE 

1 .PRIM1 I 

5.0729986D OO ) 

2.PRIM2I 3.35018370 

00 ) 

3 

6.PDDT3I 

3.31732210-01 > 

11 , ACCEL ( 1 .64797970- 

04) 

1 3 




DEPEN 

1 . OELT A 

X(-8. 888960-02 ) 

2, DELTA Y ( 3.21013D- 

01 ) 

3 

6, CELT Z 

D< 5.616260-03) 

11 , POWER ( 8. 67 1 0 OD 

00) 

1 3 


THRUST SWITCHING TIMES (DAYS) 


0.0 


POWER 

8.67100001 09 


IN IT IAL 

234S.0S74S92971 


EFFIC IENCY 
0.6337600000 


PROPULSION 

262.6012353301 


4.PODT1 <-2 .06931 4 90 00) 
15, T IMEK -1 .517625 30 03) 


4.DELT XD( -1 .1 406 30—0 1 ) 
15,T,riMEU 3.36I96D-06) 


ON 1517.625 ON 
ELECTRIC PR0PJ.5I3N PARAMETERS 
PROP TIME J 

1517.6252844700 1 . 091 5980627 

MASS COMPONENT BREAKDOWN 
PROPELLANT T ANK AGE 

443.01 C046 8725 44.301 0046873 


PROP T 1 M£ R IT 10 

1.0000000000 


STRUCTURE 

0.0 


5.PDOT2< 4.56524770 0 0) 
30. LOEGR 1—4 . 65740 6 2D— 03 ) 

5.0ELT YD< 2.328850-01) 
30. T .DEGROI 2. 325040-03 ) 


AVE ACCEL 
0.0001829860 

PAYLOAD 

1595.1451724072 


NOMINAL TRAJECTORY 2 (TOTAL 


6) IN-Udl TOR IS 


5.96050-08 


oo 


I.PRImM 5.31326930 00) 2,PRIM2< 3 
6.POOT3I 3. 18813270-01 ) tl.ACCELC 1 

1. DELTA X{ 2.528050-01) 2.0ELTA Y< 
6.DELT ZD< 4.368550 — 02 ) 11. POWER < 

THRUST SWITCHING TIMES (DAYS) 0 

POWER EFFICIENCY 

8.6595241412 0.633760000 

IN IT IAL PROPULSION 

2 3 6 9.6 2 4 9 9 2 9 4 5 2 2 6 2 .2 5 36 8 8 61 5 


INDEPENDENT PARAMETERS 

.3270277D 00) 3 .PRI M3 ( -2 . 1 7 1 7 1 520 00) 4, 

.62873550-041 13.VINFK 6.6779946D 03) 15, 

3EPEN0E NT PARAMETERS 

3.1 05510-02) 3. DELTA 2 ( -3 . 77 54 1 D-02 ) 4, 

8.6 5 9520 00) 1 3 ,T ,/IN= t ( -7.253690-0 2) IS, 

.0 ON 1519.205 ON 

ELECTRIC PROPULSION PARAMETERS 
PROP TIME J 

0 1519.2050826052 1.0241766936 

MASS component BREAKDOWN 
PROPELLANT T ANK AGE 

4 429.0728049393 42.9072804939 


PDOTH-2. 069818 
TIMEK-l .519205 

DELT XD( —4 • 1 76 4 
T.TIMElt 4.7243 


PROP TIME RAT 

1.00000000 

STRUCTURE 

0.0 


OO 00) 
ID 03) 

70-02) 

90-03) 


IO 

00 


5 • PDCJT2 ( 4 . 783981 OD 00) 
30.LDEGRI— 2.8739127 D— 03) 

5. DELT YOI-1. 357670-03) 
30 « T . OEGRD ( 3.867910-03) 


AVE ACCEL 
0.000179981* 

payload 

1635.3912188969 


NOMINAL TRAJECTORY 2 (TOTAL 16) 


I NH 1 3 1 TDR IS 3 . 90 63D-03 


ll 

© *■* 


I.PR1MK S.313387SD 00) 2.PRIM2( 3. 
6 • POOT 3 ( 3. 1B78 92 10-01 ) 1 1 .ACCEL ( 1 . 


1. DELTA X{— 7.41178D-04) 2. DELTA YC 

6, DELT ZD< 8.262910-03) 11. POWER I 

THRUST SWITCHING TIMES (DAYS) 0. 


POWER 

8.4557672000 


EFFIC IENCY 
0.6337600000 


INITIAL 

2340.9387961429 


PROPULS ION 
256.0829096522 


INDEPE 

NOE NT P 

32688420 00) 

3 

,PRI M3 ( 

60990060-04) 

1 3 

. VI NF1 ( 

DEPEN 

DENT PA 

8,365560-02) 

3 

, OE- r a 

8.455770 00) 

13 

,T, V (N= 


0 ON 1206.867 Or 

ELECTRIC PROPU-5IO 
PROP TIME 
1350 .02441 l 6041 

MASS CCMPDNEN 
PROPELLANT 
373 .4807477299 


ARAMETERS 
-2.17170720 00) 

6 « 72*75960 03) 
RAMETER5 
Z( -1 .094680-02) 

1 (-2.379390-02) 

F 1375.891 ON 
N PARAMETERS 
J 

0 .869121 5491 
I BREAKDOWN 
TANKAGE 
37.3480747730 


4, 

IS, 

4 , 
15, 


POOT 1 (-2.069532 2D 00) 
TIMEK-l .51904310 03) 

DELT XO( -5 . I 536 40-02 ) 
T,TI ME1 ( 1 .6 64 0OD-C3) 

1519.048 ON 

PROP T I ME HAT 10 
0.38873050 43 

STRUCTURE 

0.0 


5 , POOT2 ( 4.78407360 00) 
30, LOEGR (-3.1 76374 BD-03 ) 

5. DELT YOI 4.915210-02) 
30.T.DEGROI 1 ,443200-03) 


AVE ACCEL 
0.0001756068 

PAYLOAD 

1674.0270639878 



INHIBITOR IS 


NOMINAL TRAJECTORY A (TOTAL 19) 


I .52590-05 


■O 

X> 



l.PRIMH 5.31320200 00 1 
6.POOT3< 3. 18897250-01 ) 


l.OELTA X ( -5.60 128D— 0 3 I 
6.0ELT Z0< 8. 5 A 559 D — 03 ) 


INDEPENDENT PARAMETERS 

2 ,PRI M2 ( 3 .32 09Z87D 00 > 3. PR 1 M3 (-2, 171 5 7690 0 0> 

l 1. ACCEL ( l • 625 8 74 OD — 04 ) 13 ,VI NF1 I 6. 71 31 835D 03) 

dependent parameters 

2. DELTA Y I — 1 » 72 5 1 80 - OH ) 3. DELTA Z( 2.423120-04) 

11. POWER ( 8.565550 00 I 1 3 .T . Y I N*l I -2 . 650 650-02 ) 


THRUST SWITCHING TIMES (DAYS) 


0.0 


POWER 

0 . 5655498465 


IN IT I AL 

2348.0346249507 


EFFIC IENCY 
C. 6337600000 


PROPULSION 

259.4076771027 


ON 1236.139 OFF 1 372.865 ON 
ELECTRIC PROPJ_SION PARAMETERS 
PROP TIME J 

1 38 2. 37 70 86 98 19 0,095214 76 29 

MASS COMPONENT BREAKDOWN 
PROPELLANT TANKAGE 

380 .6594236109 38.0859423611 


4, PD3T 1 ( -2 . 068 90 2 00 00 ) 
15.TIMEM-1 .51910230 031 

4.D6LT XO(-l. 737730-03) 
15.T.TIMEK 4.810760-04) 
1519.102 ON 

PROP TIME R XT ID 
0 . 9099960737 

STRUCTURE 

0.0 


5.PDOT2C 4.78270080 00) 
30.L0EGR (-6. 1 9774630-03 ) 


S.OELT YD t — 1 • 227 570—02 1 
30 . T » DEGRO< -1 .4 68120-03) 


AVE ACCEL 
0.0001776305 


PAYLOAD 

1669.6815818760 


NOMINAL TRAJECTORY 5 (TOTAL 22) I NH I 31 T OR IS 


5.9605D-08 — 


t.PRIMH 5.21913620 00) 
6 . ROOT 3 ( 3 • 06 52 06 5D— 01 ) 


l.OELTA x( 1.3 351 OD- 03 ) 
6.DELT ZD( 3.88323D-03) 


INDEPENDENT PARAMETERS 

2«PRIM2( 3. 2734 061 D 00 ) 3 .PR I M3 ( - 2 . 1 659 385D 0 0) 

lit ACCEL { 1.63737730-04) 13.VINF1I 6.71453760 03) 

DEPENDENT PARAMETERS 

2.0ELTA Y I 2.145690-03) 3.0ELTA Z ( -2 . 394 79D-03I 
11. POWER < 8.62288D CO) ) 3 . T , V l N s 1 ( - 1 . 21 59 7D-02) 


THRUST SWITCHING TIMES (DAYS) 


0.0 


POWER 

8.6228773519 


INIT I AL 

2347.1 431552637 


EFFICIENCY 

0.6337600000 


PROPULSION 
261 .14 3840 6008 


ON 1232.409 OFF 1374,451 ON 
ELECTRIC PROPJ-SION PARAMETERS 
PROP TIME J 

1377 .0217396572 0.9071750874 

MASS COMPONENT BREAKDOWN 
PROPELLANT TANKAGE 

382 .7034426903 38,2703442390 


4 , POOT 1 ( — 2.028736 9D CO) 
15.T1MEK-1 .5190635D 03) 

4, CELT XD ( - 3 » 35654D— 03 ) 
15.T.TIME1 ( 4.415320-04) 

1519.064 ON 

PROP TIME RAT 10 
0.9064938549 

STRUCTURE 

0.0 


5,PD0T2( 4.68436390 00) 
30. LDEGR (-7.2061 8160-03) 


5.DELT YD ( -1 . 5 36 9 7 D— 0 3 > 
30 « T , OEGRD ( -2 . S31 45D— 03 ) 


AVE ACCEL 
0.0001789778 


PAYLOAD 

1665.0255274835 


NOMINAL TRAJECTORY 6 (TOTAL 27) 


INHIBI TOR IS 7.45060-09 


1 .PRIMK 5.1 599 1 53 D 00 
6 * PDOT 3 ( 2.89 8B390D-01 

1, DELTA X( 7. 452710-04 
6.DELT Z D( 1.701280-03 
THRUST SWITCHING TIMES ( 

POWER 

8.6591006747 
INIT IAL 

2347.4620903736 


INDEPENDENT PARAMETERS 
2 1 P R I M2 ( 3.2549 972D 00) 3‘,PRIM3< -2.1596 53* 

1 1 , ACCEL C 1 . 64 4 0 32 30- 04) 13 .VINF1 ( 6.714117* 

DEPENDENT PARAMETERS 

2. DELTA Y<-6. 392670-04) 3,Dc_TA Z(-l. 51191 
11. POWER ( 8.6591 OD 00) 13.T .YIN 1 I (-3.470 4< 


OAYS) 


0.0 


EFFICIENCY 
0 . 63 37 60 000 0 


PROPULS ION 
262.2408639344 


ON 1237.403 OFF 1369.671 

ELECTRIC PR OP J » SION PARAMETERS 
PROP TIME J 

1 386 .6839462631 


00) 

4. PDOT 1 (-2 .0 10462 OO 

00 

03) 

15 1 T I ME 1 ( - 1 .51895160 

03 

03) 

4.0ELT XO(-7.33O03D 

-04 

03) 

15.T.TIME1 ( — l .053600 

-05 

ON 

1518.952 ON 



0.9231469282 


MASS COMPONENT BREAKDOWN 


PROPELLANT 

387.0664540193 


TANKAGE 

38.7066454019 


PROP TIME RAT 10 
0.9129217404 


STRUCTURE 

0.0 


5 » PDOT2 ( 4.62858000 00) 
30, LDEGR (-5. 905 19820*03) 

S.DELT YD ( —3 .717060-03) 
30 » T . DEGRD ( — 1 .24976D-03) 


AVE ACCEL 
0.0001799027 

PAYLOAD 

1659.4481270183 



NOMINAL TRAJECTORY 7 < T OTA L 32) 


iNHiairoR is 


1 .4901D-08 


l.PRIMK 5 • 1336 404 D 00) 
6 » PDOT 3 < 2. 8387 8630-01 ) 


1 » DELT A X( 
6, DELT ZD( 


2. 359270-04 > 
7.414300-04) 


thrust Switching times (days) 


INDEPENDENT PARAMETERS 

2.PRIM2I 3.26293130 00) 3 ,PRl M3 ( - 2 . 1 55 77 1 40 00) 

11. ACCEL ( 1.64&04170-04) 13.V1NEK 6.71374320 03) 

DEPENDENT PARAMETERS 

2. DELTA Y( - 1 . 07608D- 04 ) 3,0ElTA Z ( - 5 . 52 1 00D- 04 ) 
11. POWER ( 8.67053D 00) 1 3 ■ T . V I Mf 1 ( - 6 . 54 91 1 0-0 4) 


0.0 


POWER 

8.6705313436 


IN IT l AL 

2347.6914743591 


EFFIC IENCY 
0.633760 0000 


PROPULS ION 
262.5870417409 


ON 1239.035 OF- 1368.048 ON 

ELECTRIC PROPULSION PARAMETERS 
PROP TIME J 

1389. 758596 7S 65 0.92e0672l61 

MASS COMPONENT BREAKDOWN 
PROPELLANT TANKAGE 

386 .42 966794 4 9 38 . 34 296 679 45 


4. POOT1 (-2.01417280 00) 
15, TIMEK-l .513771 OD 03) 

4 . DELT XDl-3.291 4 20-04 > 
IS, T, U ME l (-2.2338 10-C4 ) 
1518.771 ON 

PROP TIME RAT IO 
0.91505471 69 

STRUCTURE 

0.0 


5. POOT2 < 4.604B40SD 00) 
30 .LOEGR C-S .1 41 29970-03) 

5. DEL T YD( -3. 090670-03 ) 
30 , T . D£GRD( —4 . 997740-04 ) 


AVE ACCEL 
0.0001801835 


PAYLOAD 

1657.8317978767 


00 

o 


NOMINAL TRAJECTORY 8 (TOTAL 36) 


IN-4181 TOR IS 3.7253D-09 


l.PRIMK 5 . 10 68 403D 00 
6 * POOT 3 ( 2.9845884D-01 

1 .DELTA XI 2.397290-04 
6 .DELT Z 0 ( 4. 0151 70 -04 

THRUST SWITCHING TIMES I 

POWER 

8.6719973228 

INITIAL 

2347.0578738041 


) 2,PRIM2( 3. 

> 11. ACC EL ( 1 . 

> 2. DELTA r<- 

) 11. POWER < 

DAYS) 0. 

EFFIC IENCY 
0.633760 0000 

PROPULS ION 
262.6314389196 


INDEPENDENT P 
29632600 00) 3,PRIM3< 

64 6764 5D-04 ) 1 3 . V 1 N = 1 < 

DEPENOENT PA 
1.272850-04) 3 .DELTA 

6.672C0D 00) 13.T.VIN- 

0 ON 1238.824 OF 
ELECTRIC PROP JLS I 3 
PROP TIME 
1 393 .3738931363 

MASS C0MPQN5N 
PROPELLANT 
389 .71 50553729 


arameters 

-2. 1 550970 D 00) 
6.71 47767D 03) 
RAMETERS 
Z(-8. 130750-05) 

1 (-1.331910-04) 

F 1363.853 ON 
N PARAMETERS 
J 

0.932450501 1 
T BREAKDOWN 
TANKAGE 
38.9715055373 


4, 

15. 

4, 

15. 


P03T 1 ( -2 .035820 80 00) 
TIMEH-1 .51840320 03) 

DELT X0(-1 .093890-04) 
T.TI MEl ( -1 .4870 ID-04) 
1518.4JJ on 

PROP TIME RAT 10 
0. 91 765740 35 

STRUCTURE 

0.0 


5.POOT2C 4.58725360 60) 
30. LDEGR ( -4 .82 285010-03 I 

5. DELT YD(-1. 961720-03) 
30 * T « DEGRD ( — 1 .801870-04 I 


AVE ACCEL 
0.0001803267 

PA YL 0 AO 

1655.739 8739742 


NOMINAL TRAJECTORY 9 (TOTAL 39) 


I NH 131 TOR IS 


1.16420-10 — 


l.PRIMK S.0760147D 00) 
6 « PDOT 3 ( 3 . 269 8 985D — 0 1 ) 


1. DELTA X( 4.840580-04) 
6. DELT ZD( 1.523960-04 ) 


THRUST SWITCHING TIMES (DAYS) 


INDEPENDENT PARAMETERS 

2, PR I M2 ( 3 . 3463 43 OD 00 ) 3 .PR! M3( -2 . 1 58 2 7 550 0 0) 

11. ACCEL ( 1 .64774060- 04 ) I3.VINFK 6.7174625D 03) 

DEPENDENT PARAMETERS 

2. DELTA Y(-9. 147750-04) 3, DELTA Z ( - 1 .9 39 1 4 D-05) 

11. POWER ( 0.671C5D 00) 1 3 ,T . V ( N= 1 ( - 1 . 43 676D-0 4 » 


4, PD0T1 (-2 .06716460 00) 
15, T IME1I-1 . 51788330 03) 


0.0 


POWER 

8 .6 7104993 30 
INITIAL 

2345.4111595641 
THIS CASE IS CONVERGED. 


EFFIC IENCY 
0.6337600000 


PROPULSION 

262.6027472223 


ON 1 237. 601 OFF 1356.688 ON 

ELECTRIC PR GPU. SION PARAMETERS 
PROP TIME j 

1398.7961204634 0.9390740281 

MASS COMPONENT BREAKDOWN 
PROPELLANT TANKAGE 

391. 54 90 83 8753 39.154 90 8 38 75 


4. DELT X0( 
IS. T.TI ME1( 
1 517,883 


2. 631540-04) 
9.6415 8D-06) 
ON 


PROP TI ME RAT 10 
0.921 54392 60 


STRUCTURE 

0.0 


5. POOT2 ( 4.56691700 001 
30, LOEGR ( -4.66131 070-03) 

5, DELT Y0( -7. 398 4 6 D— 04 ) 
30 , T . DEGRD ( -3 .7 04 22D-06 ) 


AVE ACCEL 
0.0001 005308 


P A YL DAD 

1652. 1 044200790 


40 TRAJECTORIES WITHOUT PARTIAL DERIVATIVES AND 


9 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FOR THIS CASE 



X31NAL PAGE IS 

O /^kTV a 


S 


.O 


cas e 1 

TIME 

Rl 

LI 

LG 

PSI 

s 




SWITCH 

POINT SUMMARY 

SEMI-MAJOR 

AXIS ECCENTRICITY 

I NCLINATION 

NODE 

R2 

R3 

VI 

V2 

L2 

L3 

L 4 

L5 

LC 

L PHI 

CONE 

CLOCK 

THETA 

PHI 

LATITUDE 

LONGITUDE 

LS 

DENSITY 

OPOWR 

DPOWO 


ARG PQ5 RHAG 

V3 HASS ratio 

L6 L7 

hmag power fnct 

FLT PTH ANGLE VMAG 

degrad chi 


page 1 

travel 
thrust acc 

HAM 

switch fnct 

PROP TIME 
CHI REF 


start of trajectory, thrust on 


EARTH 

0.0 

4.963849430-01 
5 .073909600 00 
0 .0 

-1 .594890530 01 

0.0 


1 .8805164 70 00 
-8,864045150-01 
3, 3 4 88 387 5D 00 
0.0 

9.409532320 01 
-4.655926720-03 


4.S9968606D-01 

0 .0 

-2.1588149ZD 00 

0.0 

9.39374288D 01 
9.687549330-01 


3.610570230 00 
1 .033722943 30 
-2.068522770 00 
7.631060770 01 
0.0 

-1.245431433 30 


1.19 246 5 440 02 
6. 020455560-01 
4.565730810 00 
6.672519420 01 
-6.07534559D 01 
6.530838.3 60-01 


1 .800000000 02 
-7 .54 731 520 D-02 
3.310356990-01 
1.21 7642600 00 
-9.682721 730-01 
1 .000030000 00 


1 .0 1599043D 00 
1 .000000000 00 
1 .000000000 00 
9.798629180-01 
1 .198640130 00 
0.0 


0.0 

2 .722821470-02 
-1 . 37300066 D— 02 
5.238368410 00 
0.0 
0.0 



switch thrust off 

4.906303780-01 1.088257373 31 1.25632040D 02 6.67674540D 01 2.73670624D 00 2.529460510 02 
4 ,747e4359D-0I 3.965664360-01 -4.460203160-01 - 1 . 20 1 1 54 030- 02 8 .71 3523 71 D- 01 6.032873110-03 
7.031355300-01 1.929923320-31 1.901253400-01 -9.358536950-01 1.770179140 00 -1.373012670-02 
°*° 8.092104460 01 1.301452630 02 1.42395914D 00 1.907430540-01 -1.270511470-14 
8.823744160 01 9.990652720 30 —1.679797550 02 — 2.93507184D 01 5 .9694 50 0 6D- 01 1 .236831020 03 
1 .335191740-01 -1.27 8136903-01 1 .309861 140 00 9.726 068000-01 9.000000000 01 9 . 0000 OOCOD ■ 01 


SWITCH THRUST UN 

1 .35715295D 03 2.670456720 00 4.906303780-01 1.086257373 31 1.25632040D 02 9.638366720 01 2,067645740 00 2.825622640 02 
•1 .506145910 00 -1.36241 769D 00 3 .879 444 S2O-0 1 7 . 1 06 320 6 7D- 3 1 -2.856 1 182 30-0 1 -7.90 536 0 060-02 8 .7 1 35237 1 D- 0 1 1 . 004 ZS66SD-02 
7.141749260-01 -8.623381710-01 -1.171796060 00 1.305279230-31 3.725147820-01 -8.547469840-01 1.77017914D OQ -1.373012670-02 
■1.952116370 01 -1 . e 846 7069D-0 1 0.0 7.29375879D 01 5.706362440 01 1.423953140 00 3.149087240-01 -8.239936510-17 
■4.449965370 01 9.855055800 01 9.608756000 01 1.081428250 01 -1.378683790 02 -2.6561244BD 01 7 .69949037D-01 1.236831020 03 
1.97788487D 02 -2 . 2762 £84 8D-04 2.339094340-01 -2.649232773-01 1.170896S2D 00 9 .7260 680 8D- 01 0.0 0.0 


1.23683102D 03 2.670496720 00 

—2 .636 1 1 2 1 1 D 00 -5 . 6 1 29652 8D- 0 1 
3.753070330-01 -1.41086288D 00 
-1.95211637D 01 -1.884670690-01 
2 . 14624934D 01 8.81 06 068 3D 01 

1.977884870 02 -2 . 2762584 8D-04 


END OF TRAJECTORY. THRUST ON 


TEMPEL IK 19881 

1 .517856370 03 
8.867712800-01 
6.32371 1 340-01 
-2.05246983D 01 
—5 . 1 1364 2 66 D 01 
2.S9774436D 02 


3.036730040 00 
-1 . 0603228 4D 00 
9.3981671 2D-0 1 
— 2.6134328 50— 01 
9.837360070 01 
2.21 14851 50-08 


5 .44431 64 60-01 
-5.703626660-02 
-2 .184943770 00 
0 .0 

9.524279680 01 
5.22493074D-01 


1 .243166500 31 
7.862621480-01 
-4.080751530-01 
7.957 856500 01 
-2.362855150 00 
-5.790674790-01 


1.191181310 02 
6.695707070-01 
8.725228400-01 
2.022546020 01 
-5. 00935201 D 01 
9. 030641530-01 


1 .910410760 02 
-2.232464320-01 
5.181773140-01 
1.481720 080 00 
2.633763800-04 
9.641773290-01 


1.383438100 00 
8.331846790-01 
1 .864739160 00 
6.064626290-01 
1 .05658S1 OO 00 
0 .0 


3.708293540 02 
2.022630780-02 
-1.373014410-02 
8. 164061 83D— O 1 
1.397534440 03 
0.0 



CASE 


ITERATOR SUMMARY 


INDEPENDENT PARAMETERS 


OO 

CO 


1 »PR I Ml ( 

5.0739096D 00 ) 

2 *t 

6.PDOT3C 

3,31 03 570 D — 01 ) 

7.1 

11 . ACCEL ( 

1 .64765100-04 ) 

12. 

16.T IME2 ( 

0.0 ) 

17, 

21 .THETl ( 

0.0 ) 

22. 

26 , T PET 6 I 

0.0 ) 

27. 

31. PH 1 1 ( 

0.0 ) 

32, 

36, PH I 6 ( 

0.0 ) 

37, 

41 , PR 1— A ( 

0.0 > 

42, 

46 , PD3-AI 

0.0 ) 

47, 

51 , PR 1 — B ( 

0.0 ) 

52, 

56 ,PD3-B( 

0.0 ) 

57, 

61 » PR1 — C ( 

0.0 ) 

62, 

66 , PD3-C< 

0 .0 ) 

67, 

1 .DELTA 

X< 4.6S24SD-06 ) 

2, 

6,0 ELT 20( 1.012920-06) 

7, 

11 . POWER 

< 6.67100D 00) 

12, 

16, 

(0.0 ) 

17, 

21 , 

( 0.0 ) 

22, 

26 . 

( 0.0 ) 

27, 

31 , 

( 0.0 ) 

32, 

36. 

( 0.0 > 

37, 

41 • 

(0.0 > 

42, 

46 , 

( 0.0 ) 

47, 

51 , 

( 0.0 ) 

52, 

56 , 

(0.0 ) 

57. 

61 , 

(0.0 ) 

62. 

66, 

( 0.0 ) 

67, 


,V JET< 


PH 1 2 ( 
PH I 7 < 


3 ■ 34B8387D 00 
1 .OOOOOOOD 00 
2.84 39285D 04 
0.0 
0.0 
0.0 
0 . 0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


3.PRI M3f -2. 1SB31 A9D 00 
e» ltau< o.o 

I3.VINFH 6.71 77 4 I OD 03 

10.VE-OK 0.0 

23 ,THET3t 0.0 

2B , TH 2 T 9 ( 0.0 

33 . PH I3( 0.0 

38 , »H I9( 0.0 

A3 ,PR3- A( 0.0 

40 »T 1 HE A< 0.0 

53 ,PR3-3( 0.0 

53 *TI ME3( 0.0 

6 3 , PR3 -C ( 0.0 

68 .T I V|£C{ 0.0 


DEPENDENT PARAMETERS 


2. DELTA Y ( — 8. 982 9 2D — 06 ) 


0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0. 0 
0.0 
0.0 
0.0 
0. 0 
0.0 


3.0ELTA 2 ( -6 . 49 573 D- 07 
8, ( 0.0 
13.T.I/INF1 <-l .06427D-06 


18 . 

23 . 
28, 
33 . 
38 . 
43 , 
46 . 
53 . 
56 , 
63 , 
68 • 


0.0 

0.0 

0.0 

0.0 


( C .0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


4.PDQT1 (-2 .068522 80 00 
9, ( 0.0 

1 4 , V I NF 2 t 0.0 
1 9, V EL 32 ( 0.0 
24 , THET 4 ( 0 .0 
29 • THET 9 ( 0 .0 
34. PHI4( 0.0 
39. PH191 0.0 
44,PD1-A( 0.0 

49 . K S AMP ( 0.0 

54,PDl-a< 0.0 

59.KSAMP( 0.0 
64,P01-C( 0.0 
69 « K SAMP ( 0.0 


4, CELT XO( 2.398510-06 


9, 

14, 

19, 

24. 
29. 
34, 
39. 
44, 
49, 
54 , 
59, 
64. 
69, 


0 .0 
0 .0 
0 .0 
0.0 
0 .0 
0 .0 
0 .0 
0 .0 
0 .3 
0.0 
0 .0 
0 .0 


( 0.3 


THRUST SWITCHING TIMES (DAYS) 


0.0 


ON 


1 236. 631 OFF 


1357,153 ON 


poke r 

e. 671 0000072 


INITIAL 

2345 .2405989059 


EFFICIENCY 

0.6337600000 


PROPULSION 

262.6012352180 


ELECTRIC PROPULSION PARAMETERS 
PROP TIME J 

1397.5344446940 0.9382774144 

MASS COMPONENT BREAKDOWN 
PROPELLANT tankage 

391 .2220623059 39,1222062306 


1517.856 ON 


PROP TIME R AT 10 
0.9207290424 


STRUCTURE 

0.0 


5 

10 

1 5 
20 . 
25. 
30. 
35. 
40. 
45. 
50 , 
55, 
60 . 
65, 
70. 


P3QT2 ( 4.S65730BD 00 
DECLNI 0.0 

TIME1 (-1 . 51 78564D 03 
VEL03I 0.0 
TMET5 ( 0.0 

LDEGR ( -4 .6559267D-0 3 
PMI5( 0.0 
PHIUH 0.0 
PD2-A ( 0.0 
, KOROP ( 0.0 

i PD2-BI 
KDRQP( 

P02-C ( 


0.0 
0.0 
0.0 
KDHOP( 0.0 


5 . DEL T YO< -6.471390-06 

10. < 0.0 

1 5, T, T I ME 1 { 7.734 4 90-00 

20, < 0.0 

25. (0.0 

30 . T, OEGRD ( 2.21149D-08 


35, 

( 

0.0 

40 , 

( 

0.0 

45, 

( 

0.0 

50, 

( 

0.0 

55, 

( 

0.0 

60, 

( 

0.0 

6 5, 

( 

0.0 

70, 

< 

0.0 


AYE ACCEL 
0,0001805291 


PAYLOAD 

1652.2950951515 


SWITCH-COUNT HISTORY 


2. 2. 4. 4,4. 4. 4. 4. 4/ 


204 THRU5T COMPUTE STEPS. 


7 COAST COMPUTE STEPS 



■-i-/AL PAGE IS 


CASe 1 EXTREMUM POINTS 3= SELECTED FUNCTIONS 


1 

time 

ECLIPTIC 

SOLAR 


COMMUNICATION 


S* ITCH 



THRUST ANGLES 



input 


ARRAY 



LONGITUDE 

DISTANCE 


ANGLE 

DISTANCE 


FUNCT1 ON 


PS I 

THETA 


PHI 


POWER 


ANGLE 

0 

0.0 

0.0 


1.016 


66.1 


0.0 

ON 

5.240 00 


-1 5.9 

94.1 


93.9 


6.5 

ON 

0.0 

4 

1.634 

2.1 


1.016 


87.7 


0.01 


5.250 00 


-IS. 8 

93.4 


93.3 

MAX 

8.5 


0.0 

4 

2.644 

3.1 

MIN 

1 .016 


60.4 


0.01 


S.25D 00 


-15.7 

93.1 


93. C 


8.S 


0.0 

4 

13.058 

16.2 


1.024 


96. 7 


o.os 

MAX 

5.30D 00 


-14.4 

96.0 


88.9 


8.4 


0.0 

5 

03.425 

00.3 


1 . 364 

MAX 

166.0 


0.30 


4.450 00 


-0.4 

69.7 


69.7 


S.5 


0.0 

4 

152.250 

114.2 


1.667 


H 5.4 


1.22 


3.83D 00 


19. 1 

62.6 

MIN 

64 .5 


3.3 


0.0 

5 

173.169 

121.2 


2.016 


10 1 .0 


1.58 

MIN 

3.30D 00 


24.6 

62.0 


64 .8 


2.9 


0.0 

4 

226.549 

135.4 


2. 372 


67.2 


2.57 


3.94D 00 


36.5 

MIN 61.3 


67 .3 


2.2 


0.0 

S 

339.72 3 

155.6 


2. 992 

MIN 

0.9 


4.01 


4.560 00 


51.5 

63.0 


73.6 


1.4 


0.0 

T 

367.048 

159.4 


3. 1 I 5 


16.9 

MAX 

4.07 


4.70D 00 


53.6 

63.7 


74 .8 


1*3 


0.0 

T 

339.71 2 

170.8 


3.679 


171 .0 

MIN 

2.70 


5.1 30 00 


61.0 

68.7 


79.9 


1,0 


0.0 

5 

546.123 

179.4 


3.694 

MAX 

170.3 


2.71 


5.1 3D 00 


61.1 

60.9 


eo .o 


0.9 


0.0 

4 

555.605 

100.3 


3.714 


167.0 


2.75 

MAX 

5. 1 3D 00 


61.3 

69.2 


80.2 


0.9 


0.0 

4 

714.156 

194.5 


3.91 3 


24.1 


4.82 


4.850 00 

MAX 

62.7 

73.5 


e2 .5 


0.8 


0.0 

5 

748.370 

197.5 


3.922 


2.0 

MAX 

4.93 


4.59D 00 


62.7 

74.3 


82 .9 


0.8 


0.0 

5 

740.744 

197. S 


3.922 

MIN 

2 .8 


4.93 


4.5 90 00 


62.7 

74.3 


82.9 


0.0 


0.0 

4 

757.403 

190.3 

MAX 

3.922 


6.6 


4.93 


4.550 00 


62.6 

74.5 


83 .0 


o.a 


0.0 

5 

759.240 

196.4 


3,922 


7.8 


4.92 


4.540 oo 


62.6 

74.6 


83.0 

MIN 

0.0 


0.0 

5 

946.76 I 

215.3 


3.742 

MAX 

172.0 


2.76 


3.230 00 


59.7 

76.9 


04.4 


0.9 


0.0 

5 

953.421 

216.0 


3.729 


168.6 

MIN 

2.75 


3.1 70 00 


59.5 

79.0 


04 .5 


0.9 


0.0 

4 

1131 .643 

236.1 


3. 207 


17.6 

MAX 

4.1 6 


1 .1 3D 00 


47.6 

83.6 


85. 7 


1.2 


0.0 

4 

I 1 57 .45 2 

239.7 


3. 103 

MIN 

6.7 


4.1 0 


8.S5D-01 


43.7 

04. 5 


86.1 


1.3 


0.0 

4 

1236.031 

252.0 


2.737 


SO. 4 


3.26 

OFF 

-I .2 7D-14 


21.5 

08.1 


88 .2 


1 .7 


0.0 

5 

1296.317 

265.5 


2.420 


89.3 


2.22 

MIN 

-2.37D-01 


***** 

****** 


***** 


0.0 


90.0 

7 

1357.153 

2 02.9 


2.068 


134.8 


1.24 

CN 

-0.240-17 


-4 4.5 

98.6 


96.1 


2.7 


0.0 

5 

1394.950 

297.1 


1 . B4S 

MAX 

150.0 


0.87 


3.370-01 


-54.7 

102.7 


97 .3 


3.3 


0.0 

6 

1432.034 

315.0 


1.6 36 


133.2 

MIN 

0.76 


6. 55D-01 


-58.3 

105.2 


97.8 


4.0 


0.0 

5 

1441 .045 

319.9 


1.591 


125.7 


0.77 


7.1 70-01 


-59. 1 

105.4 

MAX 

97.8 


4.2 


0.0 

4 

1444.313 

321 . 3 


1.579 


123.8 


0.77 


7.330-01 


-59. 1 
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4.9630494D-O1 —8.064 8451 D— 01 0.0 
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CASE X 


( CONVERGED) 


performance summary 


EARTH 


TO TEMPEL 11(1 980 ) WITH FIXED ARRIVAL EXCESS SPEED 


ARRIVAL AT TEMPEL IM198BI PERIHELION 


LAUNCH VEHICLE IS TITAN 


III E/CENTAUR 


(COEFFICIENTS * 167238. 950 0 3480.8030 1753. 6965) 


LO 


JUL 21* 19e4. 

JUL IAN 


20.6671 HOURS GmT 
DATE 45903.3611 


AD a SEP 16, 1988, 

JULIAN 


17.2200 HOURS GMT 
DATE 47421. 217S 


flight time s isi7.856* days. 


ELECTRIC PROPULSION SYSTEM PARAMETERS 

ALPHA A (KG/KW) ALPHA T (KG/KW) TANKAGE FACTO? STRUCTURE FACTOR 

15.00 CO 15.285 0 0.1000 0.0 


EFFICIENCY COEFFICIENTS 
B O (KM/SEC) E 

0.63376 0.0 0.0 


INITIAL 

2345.2406 


ELECTRIC PROPULSION SYSTEM MASS SUMMARY (KG) 

POWER PLANT PROPELLANT TANKAGE STRUCTURE 

262.6012 391.2221 39.1222 0.0 


NET MASS 
1652. 29S1 


oo 


ELECTRIC PROPULSION SYSTEM PERFORMANCE SUMMARY 


P(1 AU) (KW) P(HSKP) (KW) 

8.6710 0.0 


P( TARG) (KW) 

5.2586 


HR( 1 AU) ( N) 
0.386461 


ACC(1 AU) (M/SECW+2) 
1 .6478510-04 


ISP (SEC) EFFIC 

2900.000 0.63376 


CHAR DEG (DAYS) 
7.12100000 03 
CONSTRAINED MAX 


MAX DIST (AU) MIN DIST ( AU I 

3.9221510 1 .0155371 


EXTREME TRAJECTORY AMD PERFORMANCE CONDITIONS 

MAX POWER (KW) MAX THTUST CN) BURN TIME (DAYS) DEG RD TIME (DAYS) 

8.498784 0.37876516 1397.53444 259.77444 


TRAV ANG (DEG) 
370.82935 


DEPARTURE AND ARRIVAL CONDITIONS 

OEP DECL (DEG) PARK INC (DEG) DEP VINF (M/SEC) C3 ( KM**2/SEC**2 ) ARR VlNF (M/SEC) C4 ( KM* «2/SEC**2 I 

-5.7634 28*5000 6717. 7 4104 45.126045 0.20776 0.000000 



C. MULTIPLE BALLISTIC SWINGBY MISSION 

This computer run demonstrates how the program may be effectively used 
in the investigation of all-ballistic missions (with electric propulsion absent),, 

In general, trajectories are forced to be all-ballistic by setting the pro- 
pulsion-time adjoint variable X8 to a "large" negative number (X 8 = - 1.D3} 
with respect to the mass ratio adjoint variable X7 = 1.D0. These settings (to- 
gether with the reference thrust acceleration Xll = l.D-4, jet exhaust speed 
X12 - 1.D0, and non-zero primer vector components X2 = . IDO at the launch 
planet and X42 = .IDO at the first intermediate-target) guarantee that the 
electric propulsion thrust switch-function will be maintained negative, yielding 
ballistic flight at all times. The electric propulsion specific masses ALPHAA 
and ALP HAT are also set to zero, so that the net spacecraft mass equals the 
launch vehicle payload. 

The particular ballistic mission simulation demonstrated here involves 
a March 6, 1985 launch of a 1,635 kg payload by an Atlas/Centaur launch 
vehicle (using the launch vehicle selection default -value of MBOOST = 0; actually, 
changing MBOOST alters the payload mass computation but not the or the 
trajectory). 

The primary target is the Earth, specified by MOPT3 = 3, and there is 
one intermediate target, the comet Giacobini-Zinner, specified by MOPTX - 41. 
The launch occurs fifteen days before (X15 = - 15. DO) the reference date 
(MYEAR = 1985, MONTH = 3, MDAY = 21), and the spacecraft passes Giacobini- 
Zinner 174.22 days (X48) after the reference date and arrives back at Earth 
353 days (XI 6) after the reference date. 

Ten days after passing Giacobini-Zinner a deep -space burn (of 188. 9 
m/sec) makes possible the re-targeting to Earth. This is accomplished by 
inputting TDV = 1. 00010D5 together with X64(2), X65(2), and X66(2) = 1.D0. 


OR' 
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The deep space burn velocity-increment (X64, X65, X66) initial guess is the 
zero vector, by default The deep- space burn could just have well occurred 
before arriving at Giacobini-Z inner (e.g. , TDV = 2.00005D5). The iteration 
sequence consists of a 6x6 hunt, with the independent variables being the initial 
heliocentric velocity at Earth (X18, X19, X20) and the deep-space burn velocity 
increment (X64, X65, X66), and the dependent variables being the position target- 
ing at Giacobini-Z inner (Y41, Y42, Y43) and at Earth (Yl, Y2, Y3). The setting 
MAXHAM = 0 is made to avoid an unwarranted BAD HAMILTONIAN warning 
message due to the presence of the deep- space burn. 

The detailed printout of the iteration sequence is omitted by setting 
NPRINT = 3 (compared to the default value NPRINT = 7), and three extra lines 
are added to each print -block by setting MPRINT = 2. 

Finally, the multiple ballistic swingby option is invoked by means of 
MOPT4 = -3, 42, which directs the spacecraft to swingby the primary target 
(Earth) in such a manner as to re-target to Earth (MOPT4(l) = -3), and then, 
at the second Earth passage, to again swingby in such a manner as to target to 
the comet Borrelly (MOPT4(2) = 42). Furthermore, the first Earth swingby is 
specified as unpowered (MSWING(l) = -1) having T2(l) - 480 days as the initial 
guess of the Earth -to- Earth transfer time, and the sec raid Earth swingby is 
specified to be powered (MSWING(2) = -5) with a specified transfer time from 
Earth to Borrelly ( 12 ( 2 ) - 127 days) and using the input Earth-departure helio- 
centric velocity initial-guess XSWING(1,2) - ,62D0, . 91D0, O.DO. Both swing- 
by iterations converged, and the powered-swingby incremental speed turned out 
to be -13.7 m/sec, the minus sign denoting a braking burn. The default value 
of TGO - - 1. DO resulted in the printout of several pages describing the tra- 
jectory segments following the time of primary target passage. 

It should be noted that the trajectory simulation would have terminated 
at the primary target (first Earth encounter) if MOPT4(l) 'were zero, and would 
have terminated at the second Eai’th encounter had MOPT4(2) been zero. The 
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mission, as simulated, thus consists initially of a comet-flyby, followed by a 
second comet-flyby making use of a double Earth-swingby to reach the second 
comet. 


The inputs for this computer run are listed below, and the resulting pro- 
gram output is displayed on the following pages,. 


&M INPUT X7=1.D0,X2=.1D0,X42=. 1D0,X8— 1.D3, AlPHAA»0.D0,ALPHAT«0.D0 

X64(2) J *l*DO / X6 5(2)=l.DO/X66(2)=l,r>O / X18(2)=l«DO/Xin(2)=l.DO / X20(2) a l.DO 

Yl(2)=l*DO / Y2(2)=l.nO / Y3(2> a l.DO / Y41(2>=l,DO / Y42(2)=l.DO / Y43(2)=l.nO 

MYEAR=1985 / MONTH=3,MDAY=21,MOPT2=3,MOPT3=3 / MOPTX=41, Xll=l.D-4 / X12=l.D0 

TDV-1.00010D5, X15s-15.no/X16-353. DO, X48-174.22D0 

X18=-8.3276D-2, X19=-l. 00 3 40 0 , X2 0 =- 2 . 0143D-3 

T2-480. DO/127. D0/XSWINfi(l/2)*.62O0/.91DO/0. DO- 

MO PT4 =- 3/ 4 2/MSW I NG=-1, -5/MPRI NT«2,MAXHAM=Q, NPR I NT“3 A END 
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er* o.o 


, EMUDOD* 

0.0 


• 

0.0 

9 

0.0 

0.99999999999999990 

-04. 

HOUR* 12 . 00000000000000 

* OM I = 0 

0.0 


, 

0. 0 

, 

0. 0 

PS I GN = 

1 .00000000000000 

0 .RADDDD* 1.000000000000000 

0.0 


• 

0.0 

, 

0,0 

RPER* 

2. OOOOOOOCOOOOOOO 

, 5A I* 1. 00000000000 OOOO ,SA 

0.0 


* 

0.0 

t 

o 

O 

0.0 



0.0 

t 

0.0 

STATE* 

1 . 000000900000000 . 0.0 


. 0. 3 

0.0 



STEP 1 * 0. 31 250 0000000000 OD- 01 .STEP2 

0.0 



0.0 

, 

0.0 

0.0 



0.0 

, 

0. 3 

0.0 



0.0 

, 

0.0 

0.0 



0.0 

• 

0.0 

-1 .000000000000000 


THRETa 400 .0000000000000 ,TOFF= 

-1 .000000000000000 


-1 .0000000 00000000 


-l .0000000 

-1 .000000000000000 


-1 .000000000000000 


-i .oooooeo 

-1 .OOOOOOOOOOOQOOO 


-1 .000000000000000 


-l .000 0003 

-1 .OOOOOOCOOOOOOOD 


-t .ooocooooooocooo 


TPI * 0.0 

o 

• 

o 



0.0 


0. c 

1 .000000000000000 


T2= 480.0000000000000 


, 127.00 

250 - 

0000000000000 


300.0000000000000 


350.0 OOOO 

500. 

oooo eoo oooooo 


XANGI* 0.0 


, X ANG2 

o.o 



0.0 


0. 62000000 

0.0 



0.0 


0.0 

0.0 



0.0 


0.0 

0.0 



0.0 


0.0 

0.0 



0. 0 


0.0 

0-0 



0.0 


0.0 

0.0 



0.0 


0.0 

0.0 



IBAL* O.INTPR* 


0. 


O.ITF* 


3 . ITPRNT* 0. 

jpp= 


O.LOADX* 


0 ,M AXHAMa 

0 

,mbqost= 


O.MOPTX* 


41 . 0. 


0, 


42. 

0 . 

0 . 0 

* 

3 

MPR I NT 

* 2.MPUNCH* O.MREAD* 


0 .MSN 


0, 

0, 

0 . 0 

• 

3 


3.NHUNG* 


25, NORMAL* 

0 

•NPERF* 


0, 

300, 

NSW ING= 0. NSWPAR* 


9999999999999990 - 05 , 
999 99999999 99 9 9 0-0 5 . 
999 9 9 9999 99 999 9 0-05 » 
9999999999999990 - 05 . 
9999999999999990 - 05 , 
9999999999999990 - 05 , 
999 9 999999999990-0 5 , 
9999999999999990 - 05 . 
9999999999999990 - 05 . 
9 99 9 99999999 9990 - 05 . 
9999999999999990 - 05 , 
9999999999999990 - 05 , 
999 999999999 9990 - 0 5 . 
9999999999999990 - 05 , 


AA = 0.0 « ALPHAT -■ 0.0 ,ALTITU* 

QOQOOOOOOOOODO . A50L* 0.0 • 0.0 • 

.51* 0.7600000000000000 , B1 * 0.0 ,B2* 

.CNIX* 31.57374000000000 • 0.0 , 

,CSTR= 0.0 .CTANtCs 0.29999999999999990-01 , 

R£T R* 0.0 .D?0W* 0.0 ,ECl* 

, 0.0 , 0.0 , 

0.0 .EHJDOX* 0.0 , 

,0.0 . G AM MA X* 1.00000 00 00000000 .GAP* 

.0 > OM 1 X* 195.20 73800000000 . 

, 0.0 .P0WF1X* -1.000000000000000 , 

.RAD0DX* 1.000000000000000 , 0.0 , 

»RAP = 38.00000000000000 .REVS* 0.0 , 

IX* 3.516014000000000 ,0.0 , 

.SOI* 0.0 .SQIX* 172.4922300000000 . 

• 0.0 « S PI RET* 300.00 00000000000 . 

. O.o , 1.000000000000000 . 

a 0.1250000000000000 .TCOAST* 0.0 » 

,0.0 .0.0 t 

,0.0 ,0.0 , 

,0.0 ,0.0 , 

. O.o ,TDV= 100010.0000000000 ,TG0» 

-1.000000000000000 . -1.000000300000000 , -1 .OOOOOOOCOOOOOOO • 

03000000 , -1.000000000000000 • -1.000000000000000 , 

, -1 . 000000000300000 , -1.000000000000000 • 

, -1.000000000000000 , -1 .coooooooooooooo • 

»T P 1 X* 1 60 . 1958333333332 . 0.0 , 

,TPOwER= 0.100000000000 00 00D 31.TSCALE* 

, 150.0000000000000 , 200.0000000000000 . 

. 400.0000000000300 • 450.0000000000000 • 

■ ASKING* 0 0 . 

,0.9100000000000000 ,0.0 , 

,0.0 .0.0 , 

•0.0 .0.0 « 

•0.0 .0.0 , 

,0.0 .0.0 , 

, O.o ,X0* 1.000000000000000 ♦ 

, 1.000001500000000 . 0.0 , 

I our a o .1 RK= 1 . iftL* 1 , IROT* 0. {SPIN* 

O.JPRINT* O.JT* O.KPART* 0. LAUNCH* 

O.MDAY* 21, MODE* 4, MONTH* 3.MOPT* 

0, 0 , MOPT 2= 3, HOPT3* 3.M0PT4* -3. 

« 0. 0, 0. O.MPOM* 0. 

I N3 * —1. —5. 0. 0, 0. 

♦ MT.vlASS* O.MUPOAT* 1 , M YE Aft* 1965.NOIST* 

O.NPR1NT* 3,N5Er* 0. 0, 300. 

1 . NT APE* 17 


00000000 

30300000 


09000000000 
D0900000 
* 0.0 
000 30000 


tEND 



CASE 1 


iterator parameters 


NO. INDEX 

1 18 

2 19 

3 20 

* 64 

5 65 

6 66 


INDEPENDENT VARIABLES 


VALUE 


-8. 32 76 00 0 00 00 00 00 OD-O 2 
-1 .00340000000000000 00 
-2 .01 43 00 000 0000 00 00—03 
0.0 
0 .0 
0.0 


STEP LIMIT 

1.00000000000000030 00 
1 .00000000000000003 00 

1.00000000000000000 00 

3. 00000 OOOQQOOOQOOD 03 

3.00000000000000000 00 

3.00000000000000000 00 


PERTURBATION 

9.99999999999999900-05 

9.99999999999999900-05 

9.99999999999999900-05 

1 . 00000000000000000-08 

1 .00000000000000000-08 

1 .00000000000000000-03 


WEIGHT 

1 .30000000000000000 00 
1 .0 00 OOOOO 0 0 0000 OOO 00 

1 .00000000000000000 00 
1 .OOOOOOOOOOOCOOOOD 00 
.1 .00000000000000000 OO 

1.00000000000000000 00 


DEPENDENT 

NO. INDEX VALUE 

1 1 0.0 

2 2 0.0 

3 3 0 .0 

♦ 41 0.0 

9 42 0.0 

8 43 0.0 


VARIABLES 

TOLERANCE 

9.9999 9 999 999 9 99 9 0 0-0 5 
9. 99 9 99 999 99 99 999 0 D-3 5 
9.99999999999999930-05 
9. 99 99 9 999 9999 999 3 0-0 5 
9. 99999 999 9999 9993 D-0 5 
9.99999 999 9999 9990 0-0 5 


CO 


THIS CASE IS CONVERGED. 


8 TRAJECTORIES WITHOUT 


PARTIAL DERIVATIVES ANO 


3 trajectories „,th partial oer, e.i.ves required for this cas 


11 

& *-0 

3ss 



ORIGINAL PAGE IS 
OF POOR QUALITY 


CASE 1 Su'tTCH POINT SUMMARY 


TIME 

SEMI-MAJOR AXIS 

ECCENTRICITY 

I NCL I NAT I ON 

NODE 

ARG POS 

RMAG 


R I 

R2 

R 3 

VI 

V2 

V3 

MASS 

RATIO 

LI 

L2 

L 3 

L4 

L5 

L6 

L7 


LG 

LC 

L PH I 

CONE 

CLOCK 

HMAG 

POWER 

FNCT 

psi 

THETA 

PHI 

LATITUDE 

LONGITUDE 

FLT PTH ANGLE 

VMAG 


R 1 REL 

R2 REL 

R 3 REL 

VI REL 

V2 REL 

V 3 REL 

RMAG 

REL 

S/C NUC MAG 

5/C TOT MAG 

GEO NUC MAG 

3 E 0 TOT MAG 

ANG ( V ,R ) 

ANG l V ■ X Y ) 



R1 REL ECL 

R2 REL ECL 

R 3 REL ECL 

VI REL ECL 

V2 REL ECL 

V 3 REL ECL 

RMAG 

ECL 


EARTH 




START OF TRAJECTORY, THRU5T OFF 




0.0 

9,966431370-01 

9.87735807D 

-02 

1 .012293300-51 

3 .458709 420 

02 

1 .800000000 02 

9 .922389820- 

01 

-9.622221 04D-0 1 

2.422123450-01 

0 .0 


-1 .481 254593-01 

-9.951537530- 

01 

-1.768923210-03 

1 .000 0009 OD 

00 

0.0 

1 .00090000 D- 0 1 

0 .0 


0.0 

0.0 


0 .0 

1 .000000000 

00 

0.0 

0.0 

0 .0 


9.67331 7390 01 

2.824917620 

02 

9.9343B306D-01 

1 .009919740 

00 

-9 .816716!: 00 — 02 

— 7.S8709206D 0 1 

7 .587094 l 8D 

01 

0.0 

1.658709420 

02 

-5 « 66291 3 34 D 00 

1 .0061 1 8920 

00 

-3.495828e9D 08 

-l .429000780 Oe 

3 .840250440 

07 

4.0 84 34 7400 0 1 

-1 .501755770 

01 

-1 .453723040 01 

3.796095100 

08 

1 .91253369D 01 

1.893708290 01 

1 .91 253 369D 

01 

1.893708293 01 

2.710054390 

01 

-1 .84 7243 120 01 



-3.7970S924D 07 

3.398198320 08 

-1 .648760240 

06 

-2,41 70451 90 01 

-3.39301 08 2D 

01 

2.298566570 00 

3.796095160 

08 


Gt A 

COS-ZlN|I9e5> 


END 

OF 

TRAJECTORY SEGMENT 1, THRUST 

OFF 





1 

.892200000 02 

9.966431370-01 

9 .877358070-02 

1 

• 012 295 3 OD- 3 1 

3»4 5870942D 

02 

2 .92569026D 

01 

1 

.031 030110 

00 

9 

.960727 600-0 1 

2. 6927573 1D-0I 

8 . 9096022 l D— 04 

-1 

.651 041053-01 

9. 5271901 5D- 

01 

1.551126430- 

0 3 

1 

.COOOCOOOD 

00 

8 

.905341390-02 

2.1 95455710-02 

2 .91 5005070-04 

7 

.030975570-02 

-5.991 054940- 

02 

-1 .249993590- 

0» 

1 

.000000000 

00 

0 

.0 

o 

o 

0 -0 

9 

.010628970 01 

7.627306 1 ID 

01 

9.934383060- 

01 

9 

.603559 02D- 

0 I 

9 

.1181 37630-02 

I .5 1 31 984 30 0 1 

1 .51 3225260 01 

4 

.947351943-02 

1 . 51 2730630 

01 

5.2362 7209D 

00 

9 

• 66920 588D- 

01 

-I 

.634932600 01 

-2. 3874501 OD 0 1 

-1 .40731 1070 01 

5 

.033225320-01 

-1.392004540 

01 

l .51 302097D 

01 

3 

.217683210 

0 1 

-t 

.698502370 01 

-1 .826868580 01 

1 .51 10601 30 01 

1 

.339544310 01 

8.832212 59D 

01 

4.745291 770 

01 




-8 

.552795700 00 

-3. 1 01 268090 0 1 

6.417390670-01 

l 

.558600063 00 

-3. 5204S428D 

00 

2.024222310 

o:i 

3 

.217683210 

01 


GIACOfi-Z 1 N ( 19851 




START 

OF TRAJECT 

DRY SEGMENT 2, THR 

UST 

OFF 




1 

.892200000 

02 

9-966431370- 

01 

9 

.877358070 

-02 

1 .0122953 

00-01 

3.458709420 

02 

2.925690260 01 

1 

.03X830110 

00 

9 

.960727800-01 

2.6927973 10- 

0 1 

8 

.9096022 10' 

-0 4 

-1.6510440 

50-01 

9.5271901 50 

-01 

1 .591126430-03 

1 

.000000000 

00 

6 

.905341390 

-0 2 

5.19545571 D- 

02 

2 

,91 5005070 

-04 

7.0309755 

70-02 

-5.991054940 

-02 

-1 .249993590-04 

1 

•OOOOOOOOD 

00 

0 

.0 


0.0 


0 

.0 


9.010 6289 

70 0 1 

7.5273061 10 

01 

9.934303 06 D- 0 1 

9 

.603559020- 

01 

9 

.1 181 3763D-02 

1.5131984 30 

0 1 

1 

.51 3225260 

01 

4 .947351 9 

49-02 

1 .512780630 

01 

5.296272090 00 

9 

.669205880- 

01 

-1 

.208744450 

07 

6.81 feSl 8530 

07 

1 

.332667580 

OS 

-9.8266590 

40 00 

-2.636773230 

00 

4.649802060-02 

6 

.924641930 

07 

1 

.573990770 

01 

I .334148520 

01 

-1 

.000000000 

30 

-1 . 0000000 

OD 30 

1 . 64705936D 

02 

2.615128950-01 




1 

• 2492 3929D 

06 

6.923702200 

07 

1 

.332867500 

05 

-1 .015991 2 

60 01 

-7.470970320 

-01 

4.649802060-02 

6 

.924641 930 

07 


DEEP SPACE 

BURN 

188.9 METERS/SECOND AT 199. 

22 

DAYS 







EARTH 








END OF TRA 

4- 

m 

n 

TORY, THRUST OFF 



3.680000 

OOD 

02 

9.986165720- 

01 

1 

.034847080 

-0 1 

1 

.0792551 

4 0- 

0 1 

3.4B629200D 

02 1 .799982650 02 

9.929569700- 

01 

-9 .7346 1 0 

71D- 

01 

1 .9579863 OD— 

01 

5 

.664450510 

-0 8 

-9 

.5434335 

30- 

32 

-1 .001842980 

00 -1 .085530520-03 

1 .OOOOOOOOD 

00 

4.423972 

260-02 

5.726471970- 

01 

9 

-006427640 

-04 

-4 

.7164241 

73- 

01 

9.017742720- 

02 -1 .065531000-05 

1 .OOOOOOOOD 

00 

0.0 



0.0 


0 

.0 


9 

.6129704 

00 

01 

2.027036970 

02 9. 939428060-01 

1 .00899794D 

00 

-I .648836 

690-02 

-8.304507040 

01 

8 

.304507070 

01 

3 

.26851 l 2 

73- 

36 

1 . 6962746S0 

02 -S. 93102027O 00 

1 .006379950 

00 

9.370437 

650 

02 

—4 .4948231 9D 

03 

8 

.473961320 

00 

-3 

.5510542 

20 

00 

- 1 .776539290- 

01 -5.616036930-02 

4.591465820 

03 

-4.526972 

090 

00 

-7.580248210 

00 

-1 

.OOOOOOOOD 

30 

-1 

.0000033 

OD 

30 

1.769965160 

02 -9.049320670-01 



-3.246200 

580 

01 

4.5913432 50 

03 

6 

.473961320 

00 

3 

.5163791 

90 

00 

-S. 259509930- 

01 -5.616336930-02 

4 .5 9146 58 2D 

03 


PAGE 1 

travel 
thrust ACC 

HAM 

SWITCH FNCT 
PROP TIME 
VM AG REL 

vmag ecu 


0.0 

1 .703031 06D-O2 

— 2 . 479404690-02 

— 8 . 650362 A 1 D OA 

0.0 

A . 5 88 079 63 D 01 
A .5 B8 0796 3D 01 


2.092569C3D 02 
1 .619452160-02 
-2.479404690-02 
-9.153409050 04 
0.0 

2.060510690 01 
2 • 06051 0690 01 


2.092569C3O 02 
1 .619452160-02 
-2.47940469D-02 
-9.15340905D 04 
0.0 

1 .018745010 01 
1.01 8745 0 1 0 01 


3.627S651B0 02 
1 . 701477440-02 
-2. S2 0507910-02 
—8 .855679360 04 
0.0 

3.55593883D 00 
3.555936830 00 



CASE 1 


ITERATOR SUMMARY 


INDEPENDENT PARAMETERS 


© o 

*3 & 

§§ 


C? ro 

£ £ 


Co 


<& 

w 


1 .PR I Ml C 

0.0 1 

2,PRIM2( 

1 .00000000-01 ) 

3 ,PR I M3< 

0.0 ) 

4 »PDOTl ( 

0.0 ) 

5,PDOT2< 

0.0 


6 . ROOT 3 ( 

0.0 ) 

7 . L MASS ( 

1.00000000 00) 

a ,• Li a j( 

-1.00000000 03) 

9, t 

0.0 > 

10, DECLN< 

0,0 


11 .ACCEL ( 

1.00000000-04) 

12. V JET ( 

1.00000000 00) 

13 .VINFl < 

0.0 ) 

14,VINF2< 

0.0 ) 

1 5, T I ME 1 < 

-1,50000000 01 


16.T IME2{ 

3.5300000D 02 ) 

17.IPARK C 

0.0 ) 

1 8 .VE-31 < 

-1 .401 2547D-01) 

19. VE-32( 

-9.95153750-01) 

20 , VEC03 < 

-1 .76892320-03 


21 » THET I t 

0.0 > 

22.TMET2 < 

0.0 > 

23 ,THET3( 

0.0 ) 

24 . T MET 4( 

0.0 > 

25, THETS ( 

0.0 


26 . T MET 6 C 

0.0 > 

27.THET7 ( 

0.0 ) 

2B.THZT3< 

0.0 ) 

29 . THET 9 ( 

0.0 ) 

30 • LDEGR ( 

0.0 


31. PHI It 

0.0 > 

32. PHI2< 

0.0 ) 

33, =>-U3< 

0.0 ) 

34. PH I 4 < 

0 .0 ) 

35, »H|5( 

0.0 


36. PNI6( 

0.0 } 

37. PH I 7 < 

0.0 ) 

3B , 3 H I 3 ( 

0.0 ) 

39. PH I 9 ( 

0 .0 ) 

40, PHI 1 0< 

0.0 


41 ,PR1-A< 

6. 90 534140-02) 

42. PR2-A < 

5.19545570- 02) 

43.PR3-M 

2.9150051 D— 04 ) 

44 .PDl-At 

7.03097560-02) 

45, PO 2 — A < 

-5.9910549D-02 


46 . PD3— A ( 

-1 .24999360-04) 

47, V [NFA ( 

0.0 ) 

48 ,T I ME A ( 

1 .7422000D 02) 

49 ,KSAMP< 

0.0 ) 

50 , KOROP ( 

0.0 


51 » PR 1 — 8t 

0.0 I 

S2.PR2-B ( 

0.0 ) 

S3 »PR 3 -3 1 

0.0 ) 

54, PO) -B( 

0.0 1 

55 , PD 2-B ( 

0.0 


56,P03-B< 

0.0 ) 

57 . V [NFS ( 

0.0 ) 

53 ,T I ME3< 

0.0 ) 

59 » < SAMP< 

0.0 ) 

60 ,KDROP< 

0.0 


61 .PBl-C < 

0.0 ) 

62.PR2-C ( 

0.0 ) 

63 ,PR 3 -C < 

0.0 ) 

64 , PD1-C< 

5 .57343290-03) 

65, PD2-C < 

3.030941 2D— 03 


66 .PD3-C( 

1.46433320-04 > 

67, V INFC ( 

0.0 ) 

60 * T 1 M E C ( 

0.0 ) 

69 » KSAMP < 

0,0 > 

70 , KDRGP < 

0.0 





DEPENOENT PARAMETERS 






1.0 ELT A 

X(-2. 169930-07) 

2. DELTA 

Y< 3. 0691 00-05) 

3 .DELTA 

Z< S. 664450-08) 

4, 

(0.0 ) 

5, 

( 0.0 


6. 

t 0.0 ) 

7. 

(0-0 ) 

a , 

<0.0 ) 

9, 

(0.0 ) 

' 10. 

( 0.0 


11 . 

<0.0 ) 

12. 

<0.0 ) 

13 , 

< 0.0 ) 

14, 

( 0.0 ) 

IS, 

( 0.0 


16 . 

(0.0 ) 

17, 

<0.0 ) 

la . 

<0.0 ) 

19, 

< 0.0 ) 

20, 

( 0.0 


21 . 

(0.0 ) 

22. 

(0.0 ) 

23. 

<0.0 ) 

24, 

<0.0 ) 

25, 

< 0.0 


26. 

( 0.0 > 

27, 

<0.0 > 

2B, 

<0.0 ) 

29, 

<0.0 ) 

30. 

( 0.0 


31 , 

( 0.0 » 

32. 

( 0.0 ) 

33, 

<0.0 ) 

34, 

< 0.0 ) 

35, 

( .0.0 


36, 

( 0.0 ) 

37, 

(0.0 > 

38, 

< 0.0 ) 

39, 

<0.0 ) 

40. 

( 0.0 


41, DEL X 

At-5. 717150-08) 

42, DEL Y 

A<-2. 073050-07) 

43. DEL Z 

A < 4 .28973D-09) 

44 , 

<0.0 > 

45, 

( 0.0 


46. 

(0.0 ) 

47, 

(0.0 > 

48 , 

o 

o 

49, 

<0.0 ) 

50, 

( 0.0 


51 . 

( 0.0 ) 

52, 

(0.0 ) 

53 . 

<0.0 ) 

54, 

(0.0 ) 

55, 

< 0.0 


56, 

(0.0 > 

57, 

(0.0 ) 

58 , 

< 0.0 ) 

59, 

<0.0 ) 

60. 

( 0.0 


61 , 

(0.0 ) 

6 2, 

<0.0 ) 

63 , 

< 0.0 ) 

64, 

<0.0 ) 

65. 

( 0.0 


66 . 

<0.0 ) 

67. 

<0.0 ) 

68 . 

<0.0 ) 

69, 

< 0.0 ) 

70, 

( 0.0 



THRUST SWITCHING TIMES (DAYS) 0.0 OFF 189.220 VISIT 199.220 ON 368.000 OFF 


POWER 

18175.6839354431 


ELECTRIC PR OPUL 5 ID N PARAMETERS 
EFFICIENCY PROP TIME J 

0.0000000045 0.0 0.0 


PROP TIME RAT IO 

0.0 


AVE ACCEL 

0.0001000000 


INITIAL 

PROPULSION 

MASS COMPONENT 
PROPELLANT 

BREAKDOWN 

TANKAGE 

1634,7360602603 

0.0 

0 .0 

0.0 


STRUCTURE PAYLOAD 

0.0 1634.7360602603 


SWITCH-COUNT HISTORY ALL 5 


0 THRUST COMPUTE STEPS 


52 COAST COMPUTE STEPS 



CASE l 


extremum points o- selected functions 


1 

TIME 

ECLIPTIC 

SOLAR 



LONG I TUDE 

DISTANCE 

0 

0.0 

0.0 


0.992 

4 

02.559 

93. 1 

M IN 

0.898 

5 

179.513 

200.2 


1.016 

4 

109.220 

20 9.3 


1. 032 

0 

109.220 

209.3 


1.032 

4 ' 

199.220 

21 6.3 


1 . 047 

4 

267.702 

275.6 

MAX 

1.102 

4 

368.000 

362.0 


0.993 


COMMUNICATION 


ANGLE 

D15TANCE 


3.4 

0.0 

OFF 

54 .6 

0.23 


77.6 

MAX 0.47 


79.9 

0.46 

OFF 

79.9 

0.46 

OFF 

82.5 

0.45 

OFF 

106.2 

0.29 


101.6 

0.00 

OFF 


Srf I TCH 


FUNCTION 

PS I 

■B.55D 

04 

4*444 

•8.1 ? D 

04 

4 4 4*4 

■9.030 

04 

4 4 44 * 

9,1 50 

04 

**** * 

9.1 50 

04 

4 4 4* 4 

9.2 7D 

04 

4444 4 

9.72D 

04 

4*44* 

B.3SD 

04 

4**4* 


thrust angles 


THETA 

PHI 

* * * 4 * * 

* *4 4* 

444*44 

4 4* *4 

44**44 

4 4* 44 

44444* 

44444 

444444 

4***4 

4*4444 

4 44 4 4 

4444*4 

***** 

44**44 

***** 


INPUT 


ARRAY 

POWER 


angle 

0.0 

ON 

90.0 

0.0 


90,0 

0.0 


90.0 

0. 0 


90.0 

0.0 


90.0 

0.0 


90.0 

0.0 


90.0 

0.0 

ON 

90.0 


CASE 


MIS5I3N SCHEDULE 


CO 


£S4£121.*.2£fle_£a_yyU.ANJMI£ DEPART EARTH 


X 


V 


2 


XDDT roor ZDOT RADIUS LAT. 


PLANET 

S/C 


-9 .622221 OD-O 1 £.42212340-01 0.0 

-9 .62 2221 OD-O 1 2.4 2212 340-01 0.0 


-2.63453550-01 
-1 .401 254 70-01 


-9,73554920-01 
-9 .95 1 S3 7 50 — 01 


0.0 

— 1 . 76B9232D— 03 


9.92238960-01 

9.92230900-01 


0.0 

0.0 


5E£IE*iSB_llA_12S5__li22SSQflSQQ_21_0 Jl M J ,3 Jt 

£±4£j!L22i22Se_fi2_JyLiAN_Q£3;£ PASS GI ACDO-ZI N < 1985 > AT 20.605 rCW/SEC 


X 


V 


z 


XDO T 


YDOT 


ZOOT 


RADIUS LAT. 


PLANET 

S/C 


9 .96 072850-0 1 
9 .96072790-0 I 


2.6927994D-01 
2. 6927973 D- 01 


8. 9095593 D— 0 4 
8. 90960220-04 


— 2.1743251 D— 0 1 
-1 .651044 1 D-01 


1 .0 7091490 00 
9.52719010-01 


-6.733 5209D-01 
1.56112640-03 


1.03183020 00 
1. 03183010 OO 


0.049 

0.049 


TWO-BODY TRANSFER ANGLE BETWEEN 


EARTH 


AND GIACOB-ZINI1985) IS 209.2569 DEGREES. 


— MAset! — a*_a2S£^_i*2flfl£aflflse_ai_SjjjjLiA 

■22£6£t22*£fl2e_££_jykIAN_eAX£ pass earth at 3.556 kw/sec 


X 


Y 


Z 


XOOT YDDT 


ZDOT RADIUS LAT. 


PLANET -9.73460850-01 
S/C -9.73461070-01 


1.9 S 76794 D— 01 
1 .9S79863D-01 


0.0 

S. 66445050-08 


-2.1 349335 D-01 
-9.543433 5D-02 


-9.6418468D-01 
— 1 .00164300 00 


0.0 

-1 .83553050-03 


9. 92950710-01 
9 * 92 9 569 BD- 01 


0.0 
0.0 00 


TWO-0CDY TRANSFER ANGLE BETWEEN EARTH AND EARTH 


IS 2.7503 DEGREES. 


LONG. 

165.871 
165. B71 


LONG* 

15.128 

15.128 


LONG. 

168.629 

168.627 



CASE 1 


C CONVERGED) 


PERFORMANCE SUMMARY 


A 


TO EARTH 


WITH VISITATION OF GI AC03-ZI N( 1 985 ) 


LAUNCH VEHICLE IS ATL AS (SLV3X J/CENTAUR 


< COEFFICI ENTS « 77360.1300 3652.7910 1653.7160) 


LD a MAR 6, 1985, 12.0000 HOURS GMT AD *» MAR 9, 1986, 12,0000 HOURS GMT FLIGHT TIME » 368.0000 DAYS. 

JULIAN DATE 46131.0000 JULIAN DATE 46499.0000 


ALPHA A (KG/KW) 

0.0 


ELECTRIC PROPULSION SYSTEM PARAMETERS 


ALPHA T (KG/KW > 

0.0 


TANKAGE FACTOR 
0.0300 


STRUCTURE FACTOR 

0.0 


EFFICIENCY COEFFICIENTS 
B D (KM/SEC) E 

0.76300 13.00000 0.0 


INITIAL 

1634.7361 


ELECTRIC PROPULSION SYSTEM MASS SUMMARY <KG) 


POWER PLANT 

0.0 


PROPELLANT 

0.0 


TANKAGE 

0.0 


STRUCTURE 

0.0 


NET MASS 
1634.7361 


P(1 AU) <KW) 
1 0175.6839 


PIHSKP) (KW) 

0.0 


ELECTRIC PROPULSION SYSTEM PERFORMANCE SUMMARY 


Pf TARGI (KW) 
18339.2277 


THRU AU) ( N 1 
0.163474 


ACCIl AU) <M/SEC**2) 
1.0000000-04 


ISP (SEC) 

0.1 02 


EFFIC 

0.00000 


CHAR DEG (OAVS) 
1.00000000 30 


MAX D 1ST <AU > 
1.1019561 


MIN DIST (AUI 
0.89820 1 1 


EXTREME TRAJECTORY AND PERFORMANCE CONDITIONS 


MAX POWER < KW) 

0.0 


MAX THRUST (N) 

0.0 


BURN TIME (DAYS) 

0.0 


OEGRD TIME (OAYS) 

0.0 


TRAV ANG ( DEG ) 
362.75652 


,o o 
^ w 

% £ 

Is 

§g 

£ o 

*3 fc*J 


DEP OECL (DEG) 
-5.1246 


PARK INC (DEG) 
26.5000 


OEPARTURE AND ARRIVAL CONDITIONS 


DEP VINF ( M/SEO 
3407.29049 


C3 (KH**2/$£C**2> 
1 1.609629 


ARR VINF (M/SEC) 
355S. 93883 


C4 ( KM**2/SEC**2 I 
12.644701 



SWINGBY CONTINUATION ANALYSIS 


THIS CASE IS CONVERGED. 

11 TRAJECTORIES WITHOUT PARTIAL DERIVATIVES AND 4 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FOR THIS CASE. 


EARTH SWING 3Y CONTINUATION TO EARTH 


PASS DIST 

IRADI I 

) SPEED I M/SEC 

> INCLIN (DEG) 

N03E (DEG) ARS 

PER 

< DEG > 


LEG TIME (OAYS) 

MISSION TIME (DAYS) 

ARR VINF (M/SEC) 

60 . 

0441 


3037 .49 

162.1681 

185.7729 

108 

.232 



527.95 

095.95 

3557.35 

ARRIVAL 

VOO 


1. 100590570-01 

— 1 . 765 030020-02 

-1 . 88553O52D-0 3 


MAG 

= 

1 

.193072350-01 

(ECLIPTIC REFERENCE 

SYSTEM) 

DEPARTURE 

voo 

s 

1. 140582540-01 

-3.527078530-02 

5. 0058 599 2 0-07 


MAG 


1 

. 1 938724 2D— 0 1 

( ECL IPTIC REFERENCE 

SYSTEM) 

ARR IV AL 

VOO 

= 

9.800750470-0 1 

- 1 .479077740-01 

-1 . 5793401 40-02 


MAG 


1 

,000000000 00 

(ECLIPTIC REFERENCE 

SYSTEM) 

DEPARTURE 

VOO 

- 

9.553630390-01 

-2.954317790-01 

4. 1 9296052D-05 


MAG 

= 

1 

.000000000 00 

(ECLIPTIC REFERENCE 

SYSTEM) 



HELIOCENTRIC APPROACH ANGLE * 177.0. 
SWINGBY INCLINATION W.R.T* ECLIPTIC = 

DEPART ANGLE = 174.2, 
-6.0 DEGREES. 

SEND ANGLE * 

8.7 OEGREES » 


<LQ 

POWERED SWINGBY INCREMENTAL SPEED ® 

0.0 METERS/S6C0ND. 

3ENO ANGLE - 

8.7 DEGREES. 

<P_An£TOCENTR1CI 


CD 


POWERED SWINGBY ANALYSIS ONLY. FOR FIXEO SW INGBY LEG FLIGHT TIME 127,0 OAYS . 

THIS CASE IS CONVERGED. 

5 TRAJECTORIES WITHOUT PARTIAL DERIVATIVES AND 2 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FOR THIS CASE. 


ARRIVAL 

VOO = 

1 .050964 I9D— 01 

-5.673956340-02 -4.655296630-07 

MAG = 

1. 19434649D-01 

I ECL 

IPTIC 

REFERENCE 

SYSTEM) 

DEPARTURE 

VOO S 

8. 326140160-02 

0 . 240 4 065 3D — 02 -1,798177080-0 2 

mag = 

1.185168140-01 

( ECL 

IPTIC 

REFERENCE 

SYSTEM) 


ARRIVAL 

VOO a 

0. 799491610-01 

-4.750676620-01 

—3 . 89 77 772 7D-0 5 

MAG a 

1 .OOOOOOOOD 

00 

( ECL IPTIC 

REFERENCE 

SYSTEM) 

DEPARTURE 

VOO a 

7.02520184 D— 0 X 

6. 952 943030-01 

-1 .517233740-0 1 

MAG a 

1.000000000 

00 

< ECL IPTIC 

REFERENCE 

SYSTEM) 


HELIOCENTRIC APPROACH ANGLE = 5.7. DEPART ANGLE = 78.3, SEND ANGLE = 73.3 DEGREES. 

SWINGEY INCLINATION W.R.T. ECLIPTIC = 9.1 DEGREES. 

POWERED SW INGBY INCREMENTAL SPEED =■ -13.7 MET ERS/SEC ONO. 3EnD ANGLE * 73.3 DEGREES. (PL AnETOCENTR I C I 


EARTH SWING9Y CONTINUATION TO BORR ELL Y I 1 98 7 » 

PASS DIST (RADII ) SPEEO (M/SECI INCLIN ( DEG I NODE I DEG ) AR 3 PER < 06G > 

3.3636 7057.96 15.9036 1 5.0632 263.592 


ARR VINF (M/SEC) 
X 7349.7* 


LEG TIME (DAYS) MISSION TIME ( DAYS I 
127.00 1022.95 



DETAILED PRINT OF P0ST-SW1NGBY TRAJECTORY SEGMENT TO 


EARTH 


FOR SOLUTION HAVING €0.04 PASSAGE DISTANCE 


TIME SEMI-MAJOR AXIS 

Rl R2 

LI L2 

LG LC 

PSI THETA 

Rl REL R2 REL 

S/C NUC MAG S/C TOT MAG 

Rl REL ECL R2 REL ECL 


ECCENTRICITY I NCL I NAT I DN 

R 3 VI 

L 3 LA 

L PH I CONE 

PMl LATITUDE 

R 3 REL VI REL 

GEO NUC MAG GEO TOT MAG 

R 3 REL ECL VI REL ECL 


NODE ARG PDS 

V2 V3 

L5 L6 

CLOCK HMAG 

LONGITUDE FLT PTH ANGLE 

V2 REL V 3 REl 

ANG (V)R| ANGIV.XY) 

V2 REL ECL V 3 R£L ECL 


PAGE 2 

RMAG TRAVEL 

MASS RATIO THRUST ACC 

L7 HAM 

POWER FNCT SWITCH FNCT 

VMAG PROP TIME 

RMAG REL VMAG REL 

RMAG ECL VMAG ECL 


EARTH START OF TRAJECTORY SEGMENT 3. THRUST OFF 

3,660000000 02 1.03625958D 00 I .09 300 439D-01 2 . &6S71 t 31 0- 5 5 1.799999910 02 3.43627474D 02 9 .929569780-01 3.6Z7565I8D 02 

-9 .734610 710-01 1 .9 579863 00-0 1 5. 664 45 OS ID -08 -9 .943 51 0390-02 -1.0X9455460 00 S. 00 58 S9 92 0-07 I.OOOOCOOOD 00 1.701477440-02 

4 . 423972260-02 5. 72647 197D-0 1 9.086427640-04 -4 . 7 1 6 424 1 7D- 0 1 9.017742720-02 -1 .06 5531 000- 05 1 .000000000 00 -2.550375910-02 

0-9 0.0 0.0 9.612970400 01 2.62703697D 02 1.011859460 00 1.008997940 00 -3.566150740 04 

9.06712656D-02 -Q. 304506190 01 8 .304507070 01 3.268511 2 70-0 6 1.686274650 02 -5,901671430 00 1 .024293310 00 0.0 

9.37043785D 02 -4.494823190 03 0.473961320 00-3.537655370 00 3 .60 1 293 00D- 0 1 1.490991200-05 4.59146502D 03 3.555939050 00 

-4.526972090 00 -7.58O24021D 00 -1 .00000000D 30 - 1 . 000 00 DO CO 30 1.741873710 02 2.40 2389 420-04 

-3.246200500 01 4.5913432SO 03 8.473961320 00 3.397215590 00 — 1 . 0505374 OD 00 1.493991200-05 4.59146502D 03 3.555939050 00 


ENO OF TRAJECTORY, THRUST OFF 


earth 

8.959453970 02 
8.383806430-01 
-5.76592725D— 01 
0.0 

—5 . 1 5304 272D-02 
3 .7 6962 €320 01 
-1 .31 5072610 01 
2.729361300 01 


1 ,036259500 00 
-5.668903470-0 1 
-7. 098020000-01 
0.0 

-9.S0224599D 01 
-7.02304520D 00 
-1.764379180 01 
-2.693306730 01 


1 .093004390-01 

1 .407953980-07 
-8 .220367550-04 

0.0 

9.502245790 01 
2.10628S070 01 

-1 .oooooeooo 30 

2 .106285070 01 


2.865711310-05 
6.490 0891 13-01 
8.094332940-01 
7.942040390 01 
7.970926233-06 
3 . 539 75 1 3 5 D 30 
-1 .000000003 30 
3.130288093 00 


1 .7999999 ID 02 
7.680909400-01 
—5 . 4 9973269D — 0 1 
2.605688390 02 
-3.406 54 9 57D 01 
3.534217540-01 
5.701724 44D 00 
-1.609983170 00 


1 .459345 13D 02 
-4.655296600-07 
-5.163421370-04 
1-01 1869460 00 
6.131090520 00 
-1 .38657021 0-05 
-2.233251 070-04 
-1 .385576210-05 


1 .01 205077D 00 
1 .000000000 00 
1 .000000000 00 
9.B4792437D-01 
1 .0 05S7260D 00 
4.3749001 2D 01 

4.374900120 01 


6.800635580 02 
1 .660659590-02 
-2. 550375910-02 
-3.50O57O58D 04 
0.0 

3.557351060 00 
3 * 557 35 1 0 6D 00 



CASE 1 


EXTREMUM POINTS 3F SELECTEO FUNCTIONS 


<£> 

OO 


I TIME ECLIPTIC SOLAR COMMUNICATION 

LONGITUDE DISTANCE ANGLE DISTANCE 


0 

0.0 

0.0 


0.992 

3.4 


0.0 

4 

82.559 

93. 1 

MIN 

0.B9B 

54.6 


0.23 

5 

179.513 

200.2 


1.016 

77.6 

MAX 

0.4 7 

4 

189.220 

20 9.3 


1.032 

79.9 


0.46 

0 

189.220 

209.3 


1.032 

79.9 


0.46 

4 

199.220 

21 8.3 


1 . 047 

82.5 


0.45 

4 

267.702 

275.6 

MAX 

1.102 

106.2 


0.29 

4 

360.000 

362.8 


0.993 

101.8 


0.00 

0 

365.000 

362 . a 


0.993 

101.0 


0.00 

4 

36a . 00 5 

362.8 


0.993 

84.2 

MIN 

0.00 

4 

434.075 

436 . 2 

MIN 

0.923 

54.6 


0.17 

5 

521.214 

530.2 


1.032 

63.4 

MAX 

0.34 

4 

626.726 

616.2 

MAX 

1.150 

169.8 


0.1 6 

9 

632.049 

620.2 


1.149 

179.5 

MIN 

0.16 

7 

632.345 

620. 4 


1.149 

MAX 180.0 


0.16 

5 

742.568 

712. I 


1.012 

82.6 

MA X 

0.35 

4 

SI 9 . 377 

796.2 

MIN 

0.923 

57.9 


0.20 

3 

895.946 

eao. i 


1 . 012 

92 ,6 

MIN 

o 

o 

o 

4 

695.946 

8 80 . 1 


1.012 

149.5 


0.00 


._AUSUSI_19a_12§I LA2Z2A§ISSe_M_a.*MAl.A 

2A4Z&26 J .94§fi_fl2_JtiU.Afl_eAlg 

X Y Z 

PLANET S • 30 30 04 60-0 1 - £.6 6090 170-01 0.0 

S/C 8 .38 380 6 4D — 0 1 -5.6609 0 350-01 1 .4079540 D-0 7 



SWITCH 


THRUST ANGLES 


INPUT 


ARRAY 


FUNCTION 

PS 1 

THETA 

PHI 

POWER 


ANGLE 

OFF 

-8.65D 

04 

***** 

****** 

***** 

0.0 

ON 

90.0 


-8.1 90 

04 

***** 

****** 

***** 

0.0 


90.0 


-9.03D 

04 

***** 

****** 

***** 

0.0 


90.0 

OFF 

-9.1 SD 

04 

***** 

****** 

***** 

0.0 


90.0 

OFF 

-9.1 50 

04 

***** 

****** 

* **** 

0.0 


90.0 

OFF 

-9.270 

04 

***** 

****** 

*** ** 

0.0 


90.0 


-9.72D 

04 

***** 

****** 

*** ** 

0.0 


90.0 

OFF 

-8 .3 50 

04 

***** 

****** 

* **** 

o 

o 

ON 

90.0 

OFF 

-3 .5 7 D 

04 

***** 

****** 

***** 

0.0 


90.0 


-3.570 

04 

**** * 

****** 

***** 

0.0 


90.0 


-3.520 

04 

***** 

****** 

***** 

0.0 


90.0 


-3.500 

04 

**** * 

****** 

***** 

0.0 


90.0 


-3.590 

04 

***** 

****** 

***** 

0.0 


90.0 


-3.690 

04 

***** 

****** 

***** 

0.0 


90.0 


-3.590 

04 

** ** * 

****** 

***** 

0.0 


90.0 


-3.530 

04 

***** 

****** 

***** 

0,0 


90.0 


-3. 5 2D 

04 

***** 

****** 

***** 

0.0 


90.0 


-3.530 

04 

***** 

****** 

***** 

0.0 


90,0 

OFF 

-3.390 

04 

***** 

****** 

* ***• 

0.0 

ON 

90.0 


PASS EARTH AT 3.557 KX/5EC 





XDOT 

YOOT ZDOT 

RADIUS 


LAT. 

LONG. 

5.439124 90-01 

8.24Q3050D-01 0.0 

1.01205050 

00 

0.0 

-34. 065 

6.490039 1 0-0 1 

7. 6 00 9094 0 -01 —4 .65529663—0 7 

1.01 2050 0D 

00 

0.000 

-34.065 


TWO-0OOY TRANSFER ANGLE BETWEEN EARTH AND EARTH IS 160.0636 DEGREES. 




DETAILED PRINT OF POST-SWINGBY TRAJECTORY SEGMENT TO BQRRELLY < 1 967 1 


FOR SOLUTION HAVING 3.36 PASSAGE DISTANCE 


TIME SEMI-MAJOR AXIS 

R 1 R2 

LI L2 

LG LC 

PS I THETA 

ft 1 REL R2 REL 

S/C NUC HAG S/C TOT MAG 

HI REL ECL R2 REL ECL 


ECCENTRICITY I NCLI NATION 

R3 VI 

L3 LA 

LPHI CONE 

PHI LATITUDE 

R3 REL VI REL 

GEO NUC MAG SEO TOT MAG 

ft 3 REL ECL VI REL ECL 


NODE ARG POS 

V2 V3 

LS L6 

CLOCK HMAG 

LONGITUDE FUT PTH ANGLE 

V2 REL V3 REL 

ANG < V * R> ANGt V • XT ) 

V2 REL ECL V3 REL ECL 


PAGE 3 

RMAG TRAVEL 

MASS RATIO THRUST ACC 

L7 HAM 

POWER FNCT SWITCH FNCT 

VMAG PROP TIME 

RMAG REL’ VMAG REL 

RMAG ECL VMAG ECL 


EARTH 

<£> 

ZD 6.959*58970 02 1.316758960 OO 

6 . 38 3806 4 3D— 0 I -S. 668903* 70-0 1 
-5.765927 25D— 0 1 -7. 09802000D-0 1 
0.0 0.0 
-9 • B20 Z5299D-0 X -9.50231976D Ol 
“I .693679870 08 -6.70264746D 07 
1.73856513D 01 1.698047380 01 

-1.199069720 08 -6.182465150 07 


START OF TRAJECTORY SEGMENT 4, THRUST 

2.3162*9270-01 9.3*1092223-31 1.45934993D 02 

1 .40 7 953 9 8D -07 6.271 7389*0-01 9.3723*5690-01 

— 8 « 220 3675 SD- 04 S . 094 332 9 4 0-0 1 -5 . 4 9973269D-0 1 
0.0 7.942840593 01 2.605668390 02 

9 .50 2 24 5790 01 7.970 92 52 30-3 6 -3.40654957D 01 

7.213193150 07 2.633309523 31 5.91323267D 00 

1 .73856509D 01 1.69S0473S3 01 1.6652B869D 01 

1 .421 4457 7D 08 2.603 831230 01 3.05278047D 00 


OFF 

1 .7 999 951 1 D 02 1.012050770 00 8.600635580 02 

-1 .798I7708D-02 1 .000000000 00 1 .66 0659590-02 

-5.163421 370-04 1.000 00 00 OD 00 6. 854707930-02 

1.116295080 00 9.64792*370-01 -3.580570580 04 

5 .90 6 4 3S52D- 01 1.103061650 00 0.0 

5 . 16989 1170 00 1.959108950 C8 2.745458260 01 

1.089650330 01 

-8.246352700 03 1.959108950 08 2.7454S826D 01 


END OF TRAJECTORY, THRUST OFF 


BORRELLYI 1987) 

1 .022945900 03 
3.802012360-01 
8.117853470-0 1 
0.0 

3.768403030-01 
4 • 865 1 5 B 1 1 D 00 
-X .701931 CIO 01 
8.27254589D 00 


1.316758960 OO 
1.3041941 70 00 
-4.1964893 60-01 
0.0 

-1.010069430 02 
-9.596949150 00 
-2.000980980 01 
2.547955290 00 


2 • 3 1 6 24 92 7D— 0 1 
-2.103789930-02 
4.265490640-03 
0.0 

1 .0 1 086 7 0 1 0 0 2 
-2 .09 22 7858 D 00 
1 .40530 1 6 7D 01 
-6 .724732440 00 


9.34 1092J23-01 
-7.345 3552 23-01 
2.103731530-01 
1 .030228420 02 
-8.893301810-01 
4.700472120 00 
1 .427062080 01 
4.722179333 00 


1 • 4 593 49 93D 02 
4, 1 601 5951D-01 
4.354921170-01 
2.62S22979D 02 
7.374743690 01 
— 1 .1 35577 88D 01 
7.428070770 01 
3.255345610 00 


2.878102260 02 
1 . 03930260D-03 
5.252351040-03 
1 .116293080 00 
1 .327035960 Ol 
-1.224603870 OL 
—4 .48969591 O 01 
-1 .637428960 01 


l .358646430 OO 
1 .0 00 00 0 OOD 00 
1 .000000000 00 
6.4681 9844D-01 
8.441640200-01 
1 .096125530 01 

1 .096 1 25S3D 01 


9.878742730 02 
X • 09373 5000—02 
6.854707930-02 
-3.895213250 04 
0.0 

1.734974390 01 
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